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ABSTRACT 


Topics in vacuum decay 


Ali Masoumi 


If a theory has more than one classically stable vacuum, quantum tunneling and thermal jumps 
make the transition between the vacua possible. The transition happens through a first order 
phase transition started by nucleation of a bubble of the new vacuum. The outward pressure of the 
truer vacuum makes the bubble expand and consequently eat away more of the old phase. In the 
presence of gravity this phenomenon gets more complicated and meanwhile more interesting. It 
can potentially have important cosmological consequences. Some aspects of this decay are studied 
in this thesis. Solutions with different symmetry than the generically used 0(4) symmetry are 
studied and their actions calculated. Vacuum decay in a spatial vector field is studied and novel 
features like kinky domain walls are presented. The question of stability of vacua in a landscape 
of potentials is studied and the possible instability in large dimension of fields is shown. Finally a 
compactification of the Einstein-Maxwell theory is studied which can be a good lab to understand 
the decay rates in compactification models of arbitrary dimensions. 
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An artistic depiction of the multiverse. There is a whole universe inside each bubble. 
Left, the potential for a particle in one dimensions. By vacuum here we simply mean 
xi,X 2 and X 3 . Right, the field potential for a single scalar field. Here by vacuum we 
mean a spatially homogeneous configuration where the field is resting at (j)i , (p 2 or (j)^. 


An example of a potential with a metastable vacuum at A and a stable vacuum at B. 
An infinite well sectioned into two pieces, left a false vacuum of width b and potential 

Vi separated by a wall of width w and potential V 2 from the rest of the well. 

Chance of finding the particle at the left of the barrier as a function of time. Left is 
for the case vi = 4, U 2 = 40, U 3 = 0, Right for the case vi = —8, V 2 = 40 and U 3 = 0. 
Time is measured in units of 2mb^h~^ and the vertical axis has logarithmic scale. In 

both cases, a = 20, w = 0.2 and 6 = 1 . 

Probability of finding the particle on the left side of the barrier as a function of time. 
The potential is described by vi = 4, U 2 = 40, U 3 = 0, 6 = 1 and w = 0.2. Time is 
measured in units of 2mb‘^h~^ and the vertical axis has a logarithmic scale. The red 
filled circles, blue squares and black empty circles correspond to a = 40, 20 and 5. . . 


The anharmonic oscillator potentials defined in Eq. (2.5) for different values of A. 

Left, middle and right graphs correspond toA = l,0,—1 . 

Analyticity domain of E{X) in the complex A plane. The inner product is not defined 
for all vectors and it has a nontrivial analytic structure. Therefore these eigenvalues 
evade the usual proof for the reality of the eigenvalues of the Hermitian operators . 
Possible exit points in a multi-dimensional potential. 
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2.8 The dominant configuration for thermal tunneling in a 2D potential is jumping to a 

saddle point with one negative eigenvalue of the Hessian matrix. [T3] 

2.9 True and false vacua in a scalar field potential. For convenience, we set V{4>f) = 0 . [15] 

2.10 The motion in Euclidean time is equivalent to motion in an upside down potential . [161 

2.11 Bubble nucleation as seen in Euclidean time from top left to the bottom right. In the 

beginning the whole space is filled with the false vacuum (shaded light blue region) 
and then a small bubble of true vacuum (green region inside the circles) with a wall 
(red region which is the boundary of the circle ) appears and grows. The top right 
configuration shows the configuration at the moment of the bubble nucleation. For 
large Euclidean times, it bounces back to a space hlled with the false vacuum. [T9| 

2.12 Another view of the bounce. In this picture the Euclidean time progresses upwards, 
the false vacuum is pictured as a white background and the true vacuum is in red. 

Any horizontal line cuts the bounce at a spatial hypersurface. In the far past (the 
lower thick dashed line), the field is almost at its false vacuum, and the cut through 
the center of the bounce (upper thick dashed line) shows the field configuration right 
after the nucleation of the bubble. . M 

2.13 The Euclidean equations of motion have drastically different behaviors for paths 

starting from the true and false vacua. On the top, the two potentials with minima 
at xi and at the bottom their upside-down versions which are more relevant in 
Euclidean spacetimes. In the left potential, xi is located at a true vacuum and there 
is no solution for the Euclidean equations of motion except for x{t) = xi. But for 
the left, in addition to this constant solution, there is another solution that starts 
from xi, moves to X 3 and then returns to xi after an infinite time. This solution is 
called a bounce. [ 22 ] 

2.14 Left is a single bounce centered at to and right a multi-bounce configuration. Both 

of these are approximate stationary points of the action and will contribute to the 
path integral in the semiclassical level. (The left bounce would be an exact stationary 
point if the bounce was centered at t = 0 ). [ 22 ] 
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2.15 The zero eigenfunction of the second variation of the action. This eigenfunction has 

a node so there must be another lower (negative) eigenvalue in the spectrum of these 
variations. 

2.16 Left, a class of paths parametrized by z and right, their Euclidean actions as a 

function of z. The z = 0 path corresponds to the constant path x{t) = xi and has 
zero action. The z = 1 path is the bounce. We chose the paths in a way that the 
tangent to the z = 1 path is the negative mode and therefore we are sampling the 
divergent region of the path integral. 

2.17 Euclidean action as a function of z for the potential shown on the left graph in 

Fig. 2.13 This function is positive and the path integral is well-dehned. To get the 
correct result for the potential on the right side of Fig. |2.13 we analytically continue 
it to the potential on the left. 

2.18 The contour of integration chosen for the path integral along z which continues from 

—oo on the real axis and extends to the saddle point at z = 1 where the bounce is 
located. Then it distorts to the upper half plane where the Euclidean action is again 
positive. 

2.19 At a nonzero temperature we can distinguish three different paths for tunneling. The 

horizontal arrow shows a purely quantum mechanical tunneling. The dashed arrows 
show a thermally assisted tunneling in which the particle first jumps thermally from 
1 to 2 and then quantum mechanically tunnels to 4. The dotted dashed arrows show 
a purely thermal tunneling in which the particle jumps to a stationary point with 
one negative eigenvalue. In all cases, after tunneling to the right of the barrier, the 
particle evolves classically towards the lower minimum. 
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2.20 The same as Fig. |2.19[ but for a field theory. This time the horizontal axis represents 

the (infinite-dimensional) field conhguration space. Points (pi{x) and 4’q{x) show the 
two homogenous configurations where the field is spatially uniform and rests in its 
vacua. 4 > 3 {x) is a configuration which has only one negative eigenvalue and it is 
in fact a bubble of true vacuum inside the false vacuum. Its negative eigenvalue is 
along the expansion/contraction directions for the bubble. (I) 2 {x) shows the bounce 
conhguration at r = 0 and r = /3 in Fig. 2.23 and (j)i{x) is the conhguration after 
the nucleation for a thermally assisted bubble. 

2.21 A quantum bubble appears with zero energy, as expected from the conservation laws. 

A thermal bubble appears with energy absorbed from the heat bath and in fact its 
energy is a maximum with respect to the radius. A thermally assisted tunneling is 
between these two extremes. 

2.22 Bounces representing thermal tunneling are cylinders in four dimension with a period 

of (3. Each slice through this bounce gives a three-dimensional ball which is a bubble 
of true vacuum separated from the false vacuum by a wall (the solid vertical edges 
of the rectangle). 

2.23 The left picture shows a low-temperature thermally assisted tunneling and the right 

a high-temperature thermally assisted tunneling. The Euclidean time r progresses 
upwards. The field conhguration at r = 0 and r = /3 is given by 02 (^) and a spatial 
slice through the center gives 04(x). (See Eig. 2.20). In the limit r —oo, this 
pictures reduces to the graph in Eig. 2.12[ as expected. 

2.24 One possible mode of decay is the huctuation of the held in every point to 0top- It 

has a hnite action and nonzero probability in de Sitter space. 

2.25 Left is CDL instanton that describes tunneling from a false vacuum to a true one. 

Right is the same instanton, but the poles interchanged by a rotation. It describes 
the CDL instanton for the decay of the true vacuum. 

2.26 A CDL bounce as it is seen in the static patch. The horizontal dashed lines are the 

constant r hypersurfaces, the shaded region the true vacuum and r ranges from 0 at 
the center of the bubbles to A at the horizon. 
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2.27 The bounce of Fig. 2.26 is shown in the Ui-y^ plane on the left. The lines r = —vrA 
and T = ttA are identified and they show the radial directions. The configuration 
after tunneling is given by half of the hyperplane that goes through the center of 
the bounce. The other half shows the configuration before the nucleation. On the 
right the same bounce is shown after a rotation. This picture shows a r-independent 
bounce which resembles a purely thermal tunneling. 


3.1 A schematic view of the thin-wall solution. The wall is located at the radius R and 

inside is the true vacuum. The horizon is located ru and between the wall and 
the horizon is the false vacuum. The bounce solution is constant in the periodic 
Euclidean time r and apparently rotationally invariant. 

3.2 The induced geometry of the static instanton in two different planes. Left is the r — r 

plane where it is described by a periodic angular coordinate r and the coordinate 
along the vertical axis is r. Horizon is at the tip of this object rn and the wall is 
at the location that there is a discontinuity in the first derivative due to the energy 
density of the wall. Right, the same picture, but in r — (/> plane. Again (p is the 
periodic coordinate and r is along the vertical axis. To avoid a conical singularity, 
the period of r must be a very specific number. 

3.3 Two views of an 0(3) x 0(2)-symmetric bounce. These are from the same viewpoints 


as those of the CDL bounce in 2.27 (The views on the right would be the same if 
viewed from along any axis in the X 4 -X 5 plane.) On the left-hand view the radial 
dashed lines correspond to the similar lines in Fig. |2.27 The bounce is independent 
of the Euclidean time. The three- dimensional configuration corresponding to any 
of these lines contains a single true vacuum critical bubble in the center of a false 
vacuum horizon volume. 
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3.4 Increasing G leaves less and less room for the wall. As a result cf) varies less over the 

horizon volume. In this graph the markers on the left of the barrier show (;/!)(0) and the 
markers on the right of the barrier (f){rH) for different G. Black circles, red diamonds, 
green triangles and the orange square show 4){0) and (pirn) in order of increasing 
G. For small G (black circles), the field starts very close to the true vacuum at the 
center and ends very close to the false vacuum at the horizon. Increasing G makes the 
change in (p smaller. Eventually, for a critical value of G which is calculated in this 
section, the two points merge(orange square). This corresponds to a Hawking-Moss 
solution. 

3.5 Oscillating bounces for different values of N. Green solid line, blue dashed, red 

dotted and black dot-dashed lines represent solutions for N = 1,3,5 and 7. Only 
the blue line represents a tunneling instanton. 

3.6 Evolution of the bounce as the gravitational constant is varied. The first two panels 

show (p as a function of (a) r and (b) r/rfj. In both cases the black dot indicates the 
value at the top of the barrier. The third panel shows / = 1/A as a function of r. 
Reading from left to right, the short-dashed purple, long-dashed brown, dotted red, 
dot-dashed blue, and solid black lines correspond to 8ttGv‘^ equal to 0.711, 0.704, 
0.628, 0.251, and 0.126. In all cases Uq = 2, a = 3, and e = 0.3. 

3.7 Evolution of the bounce as the true vacuum energy Uq is varied. The first two panels 

show (p as a function of (a) r and (b) r/rn- In both cases the black dot indicates 
the value at the top of the barrier. The third panel shows / = 1/A as a function of 
r. Reading from left to right, the short-dashed purple, long-dashed brown, dotted 
red, dot-dashed blue, and solid black lines correspond to Uq equal to 2, 5, 10, 2, and 
26.8 In all cases S'kGv'^ = 0.126, a = 3, and e = 0.3. 

3.8 Behavior of the bounce as e is varied. Again, panels (a) and (b) show <p as a function 

of r and of r/rn, with the black dot indicating the value at the top of the barrier. 
Panel (c) shows the metric function / = 1/A. Reading from left to right, the dashed 
brown, dotted red, dot-dashed blue, and solid black lines correspond to e equal to 2, 
1, 0, and —1. Eor all cases Uq = 2, a = 3, and SttGv^ = 0.126. 







3.9 Variation of the bounce with barrier height for small a with the short-dashed purple, 
long-dashed brown, dotted red, dot-dashed blue, and solid black lines corresponding 
to a equal to 1, 3, 5, 12, and 15. For all of these Uq = 2, e = 0.3, and SttGv'^ = 0.126. 

As in the previous figures, the black dot indicates the value of the field at the top of 

the barrier. [65] 

3.10 Variation of the bounce with barrier height for large a, with the short-dashed purple, 
long-dashed brown, dotted red, dot-dashed blue, and solid black lines corresponding 


to a equal to 15, 20, 25, 40, and 100. For all of these C/q = 2, e = 0.3, and 
SttGv'^ = 0.126. As in the previous figures, the black dot indicates the value of the 
field at the top of the barrier. EH 


3.11 Comparison of the tunneling exponents for the CDL bounce (red diamonds) and the 
0(3) xO(2)-symmetric bounce (black circles) as one parameter is varied with the 
other three held fixed. When held fixed, C/q = 2 , a = 3, e = 0.3, and SttGv'^ = 0.126. ESI 


4.1 


4.2 


4.3 


4.4 

4.5 


The blue (longer) and red (shorter) arrows represent the vector field at the two vacua. 
The thick green line is the domain-wall. From left to right, we show a longitudinal 
wall, a transverse wall, and a wall with orientation 9. The orientation is defined such 

that for a longitudinal wall 9 = 0, and for a transverse wall 9 = tt/2 . 

When the flat, longitudinal wall on the left figure is too massive, it will spontaneously 
breaks into zigzag segments in the right figure. Although the total wall area increases, 

the reduced tension still reduces total energy. 

The right figure is the energy stored in the wall connecting points xi and X 2 (left 
figure) vs the angle 9. The solid blue, dashed purple, dot-dashed red and dotted 
green lines correspond to clIct = 1.5, 2, 2.5 and 3. The dots show the angle with 
minimum energy. Clearly, for high ratios of cl/ct, the zigzag walls (0min 7^ 0) are 

energetically favorable. 

Bubble shape in the first quadrant. 

The bubble profile given by Eqs. (4.33) and (4.18) when cl > 2ct- The correct 
profile of the critical bubble is simply the middle portion. 
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4.6 Picking two pieces of the wall (top left panel) and glning them at angle 9 (top right 


panel). Each piece locally satisfies Eq. (4.33) and therefore gives the minimum energy 


for the fluctuations parallel to the wall. The energy of the bubble vs the gluing angle 
6 is shown in the lower panel for different values of cl /ct for a tension that satisfies 


Eq. (4.18). Blue dashes, green dots, red dot dashes and solid purple lines correspond 


to clIct = 1.5, 2, 3 and 4. The correct spot for gluing the two is the angle which 
gives a local maximum (negative mode). For c^/ct < 2 this gives a smooth bubble. 
For larger ratios it gives a kink. 


4.7 The left figure depicts Eq. (4.35), where we attempt to make a larger bubble by 


inserting true-vacuum regions and extra interpolation walls. The right figure depicts 
Eq. (??), where we try to make a smaller bubble by removing part of the walls and 
part of the true-vacuum region. Both result in smaller total energy, which shows 
that the kinky shape is indeed a saddle point. 


5.1 In a theory with one scalar field (left picture), to go from the false vacuum (right 
minimum) to the true vacuum (left minimum) we have to cross the highest barrier. 
But for the case of a field theory with two fields (right picture), to go from the false 
vacuum to a truer one, we do not need to pass the highest barrier (central peak) and 
can choose to go through many other lower barriers. This decreases the instanton 
action and enhances the tunneling rate. 


5.2 Left panel shows the potential introduced in Eq. (5.13). The progression from top 


to bottom shows the potentials corresponding to 5 = 0,0.2, 0.4,0.6 and 0.8. Case 
6 = 0 is a symmetric double well and the larger values show more and more deviation 
from the Z 2 symmetry. The right panel shows the corresponding field profiles (top to 
bottom) of bounce solution for the values of 6 = 0.4, 0.5,0.6,0.Tand 0.8. The solution 
deviates more and more from a thin wall and at 6 = 0.9 it merges to a Hawking-Moss 

instanton. 

5.3 Left panel, the median value of s for quartic potentials on a log-log plot. Right panel, 
the median value of the lowest barrier heights for quartic potential in units of Xv^. 
The bars in both cases show the 25th and 75th percentile. 
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5.4 Left panel, the median value of the distances to the saddle point which has the 

smallest s for quartic potentials on a log-log plot. Right panel, the median value of 
distances to the saddle points with the lowest barrier height for quartic potentials in 
units of v~^. The bars in both cases show the 25th and 75th percentile. [98] 

5.5 The median value of the number of saddle points around the origin on a log plot 

for quartic potentials. As is apparent from here, this number fits very well to an 
exponential. The bars show the 25th and 75th percentile. |99| 

5.6 Left panel, the median value of s for cubic potentials on a log-log plot. Right panel, 
the median value of the lowest barrier heights for cubic potentials in units of Av^. 

The bars in both cases show the 25th and 75th percentile. 11001 

5.7 Left panel, the median value of the distances to the saddle points which have the 

smallest s for cubic potentials on a log-log plot. Right panel, the median value of 
the distances to the saddle points with lowest barrier heights for cubic potentials in 
units of v~^. The bars in both cases show the 25th and 75th percentile. llOOl 

5.8 The median value of the number of saddle points around the origin on a log plot 

for cubic potentials. As is apparent from here, this number fits very well to an 
exponential. The bars show the 25th and 75th percentile. 11011 

5.9 Left panel, the median value of s for supersymmetric potentials on a log-log plot. 

Right panel, median value of the lowest barrier height for supersymmetric potentials 

in units of Au^. The bars in both cases show the 25th and 75th percentile. 11011 

5.10 Left panel, the median value of the distances to the saddle points which have the 
smallest s for supersymmetric potentials on a log-log plot. Right panel, the median 
value of distances to the saddle points with lowest barrier height for supersymmetric 
potentials in units of v. The bars in both cases show the 25th and 75th percentile. . 11021 

5.11 The median value of the number of saddle points around the origin on a log plot for 
supersymmetric potentials. As is apparent from here, this number fits very well to 

an exponential. The bars show the 25th and 75th percentile. 11021 

5.12 The dependence of the median value of s on the range of parameter 04. The data 

shown are for N = 2 . 11051 
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5.13 The distribution of values of s in the ensembles for quartic (top) and cubic (bottom) 
random potentials for various values of N. The vertical axis represents the natural 
logarithm of the number of values in each bin. For the sake of clarity, the data for 
different values of N have been offset by constants, so only the slopes are meaningful. 
Purple diamonds correspond to = 2, green down-pointing triangles, iV = 4 ; red 
up-pointing triangles, N = 6 blue squares, N = 8 and black circles, iV = 10. ... 

5.14 The distribution of values of lowest barrier heights in the ensembles for quartic (top) 

and cubic (bottom) random potentials for various values of N. The vertical axis 
represents the natural logarithm of the number of values in each bin. For the sake of 
clarity, the data for different values of N have been offset by constants, so only the 
slopes are meaningful. Purple diamonds correspond to = 2, green down-pointing 
triangles, A'" = 4 ; red up-pointing triangles, N = 6 ■, blue squares, N = 8 and black 
circles, A^ = 10. 

5.15 The distribution of values of the distance of the lowest saddle point to the false 

vacuum in the ensembles for quartic (top) and cubic (bottom) random potentials 
for various values of N. The vertical axis represents the natural logarithm of the 
number of values in each bin. For the sake of clarity, the data for different values 
of N have been offset by constants, so only the slopes are meaningful. Purple dia¬ 
monds correspond to A^ = 2, green down-pointing triangles, A^ = 4 ; red up-pointing 
triangles, N = Q ] blue squares, N = 8 and black circles, Ai = 10. 

5.16 The distribution of values of s (top) and lowest saddle point heights (bottom) in 

the ensembles for supersymmetric random potentials for various values of N. The 
vertical axis represents the natural logarithm of the number of values in each bin. For 
the sake of clarity, the data for different values of N have been offset by constants, 
so only the slopes are meaningful. Purple diamond correspond to Ai = 1, green 
down-pointing triangles, N = 2 ] red up-pointing triangles, N = 3 ] blue squares, 
N = A and black circles, N = 5 . 














5.17 Parameter ranges allowing multiverse explanations of the cosmological constant. 

With a = 2.66 and g = 10, a sufficient number of metastable vacua is only pos¬ 
sible for parameters in the region to the left of the solid line. This region is extended 

to the dashed line if g = 100. 11151 

5.18 Like Fig. |5.17[ but for a = 3.16 and <7 = 10. 11151 


6.1 Different behaviors of the potential as we change the number of fluxes. For very large 
n, there is no stationary point (dotted purple line). By lowering n, the potential gets 
a single stationary point (dashed red line) and for smaller n (black line), it gets one 
maximum at a larger value of and one minimum at a smaller ijj. In the last case 
there is a stable compactified phase separated by a barrier from the decompactified 
phase. We are only interested in the last case. 11211 


6.2 Different possibilities for the minimum of the potential introduced in Eq. (6.10) 


As we change h from very small values, the minimum changes from being de Sitter 
(black dotted-dashed line), to Minkowski (red dashed), to an AdS with a positive 
maximum (blue line) to an AdS space which is never positive(dotted purple). These 

graphs correspond to h = 0.467, l/\/2,0.776 and 0.826. 11221 

6.3 The possible situations for the critical points with all the spheres of the same size 


for the effective potential of Eq. (6.25). There is no minimum (vacuum) above the 
black dot-dashed line. The blue solid line specifies the Minkowski vacuum, below 
it are the AdS and above the dS vacua. The transition between the CDL and HM 
bounces happens at the red dotted line. Above the purple dashed line, there is a 

saddle point and a minimum and under it a maximum and a minimum. 11301 

6.4 The phase diagram of the number of minima (-1—1-), maxima(-) and saddle points(-|—) 

for the N=2 case. The two curves satisfy Eqs. and There is no sta¬ 

tionary point above the red line. On the left of the vertical dashed blue line (and 
similarly under the horizontal dashed blue line), a maximum and saddle point anni¬ 
hilate. The derivation of the curves and also the annihilations on the left of these 
lines are derived in the main text. 11321 
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6.5 The phase diagram for the N = 2 case. The green line given by Eq. (6.60) shows the 


location of the Minkowski vacua, above it is the dS and below it the AdS minima. 
There is no minimum above the top dotted-dashed blue line (region 1) which is 


given by Eq. (6.58). The red dashed line, which is given by Eq. (6.59), shows the 


boundary where the CDL bounces vanish and the HM bounces dominate . The most 
interesting part of this landscape is the region between the top red line and the green 
line where all the dS minima live in (regions 2 and 3). As proven in the text, its 
area, and therefore the number of dS vacua, is finite. 


A.l A schematic view of the de Sitter space. Each of the points here is a (d-2)-sphere. 
The dashed line shows a (d-l)-sphere at the end of the contracting and beginning of 

the expanding phase. 

A.2 Penrose diagram for the de Sitter space. The vertical solid lines are time-like points 
corresponding to the points on the north and south poles. As seen, it takes an infinite 
amount of time for a signal to travel between the poles and the north pole is only 
causally connected to the shaded region which is called the static patch. 


B.l The Cl and ct axes are the directions in which the field has purely longitudinal and 
transverse sound speeds. The two dots represent the two discrete vacua. In the limit 
/3 —0, the dashed lines through them are almost in the vacuum, too. The important 
portion of the domain-wall is the red (thick) path from one line to the other, which 

is free to pick the best orientation cj) . 

B.2 The numerically calculated values of the tension for a double-well potential are shown 


in dots. The three analytical fits: Eq .(4.18) is the dashed line, Eq. (4.19) is the 


dot-dashed (blue) line, and Eq. (B.14) is the solid (red) line. We can see that for 


large /3, Eq. (4.19) is a good fit, and for small /3, Eq. (B.14) is a good fit 


B.3 The numerically calculated values of the tension for a double-well potential are shown 


in dots. Again the three analytical hts: Eq. (4.18) is the dashed line, Eq. (4.19) is 


the dot-dashed (blue) line, and Eq. (B.14) is the solid (red) line. The two figures use 


the same potential but different sound speed ratios. 
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B.4 The more complicated potential introduced in Eq .(B.16), with q = 0.5 , ri = 

—2.5 , r 2 = 2.5 . We will use two values of S, 1.1 and 0.9, but that makes no visual 
difference. 11501 


B.5 From the potential given by Eq .(B.16), we again compare the numerical (t{9) with 
the three equations. In the left figure we have S' = 1.1. In the right figure we have 
S = 1. For some orientations the domain-wall does not exist because the path runs 
away toward the —4>x direction. I15II 
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“The desire for knowledge, I could not forego 
Few secrets remained that I did not know 
For seventy-two years, I thought night and day 
Until I came to Know, I had nothing to show.” 
Omar Khayyam 
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Chapter 1 

Introduction 


The vacuum has to be the most boring place on Earth (in the Universe), completely empty with no 
interesting events. This was the notion of people living in the pre-historic era of the 19th century, 
before they had learned to understand quantum mechanics. 

Over the last quarter of century we started discovering how fascinating and interesting this 
place is. In this previously dull environment, all of a sudden a bubble is born. It grows fast, in fact 
very fast, almost at the speed of light, and before we get alarmed it engulfs us. This bubble gives 
birth to new bubbles and we end up with a universe full of universes, a multiverse (See Figj^. 
In this thesis I will explain some topics on this “road to the multiverse”]^ 



Figure 1.1: An artistic depiction of the multiverse. There is a whole universe inside each bubble. 
Let’s start from the simple case of a particle moving in one dimension under the influence of 


^Family Guy, season 8, episode 1. 
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a potential V{x) shown in Fig jl.2[ If the particle is at rest at any of the minima of the potential, 
classical mechanics predicts that it will remain there forever. 

Similarly the vacua of a field 4>{x) are the local minima of its potential V{(f>). If V{(p), the field 
potential, takes its local minima at (/>i, 4>2 and ^3 (see Fig^, the minimum energy configurations 
are the states where the field is spatially homogenous and takes values 4>i, (1)2 or </> 3 . If the field starts 
at rest in any of these vacua, according to classical mechanics it remains there forever. However the 
situation is different when we include quantum mechanics. Quantum tunneling allows the particle 
or field to tunnel through the potential barriers and move out of a metastable minimum. The 
same is true if we put these closed systems in thermal contact with a heat bath. They can absorb 
enough energy to overcome the barrier and move towards a more stable minimum. In fact they 
will also have the opportunity to go from a more stable to a less stable minimum, but at a lower 
rate. For a single particle, there is a nonzero probability for these events to happen. For a field in 



Figure 1.2: Left, the potential for a particle in one dimensions. By vacuum here we simply mean 
xi,X 2 and X 3 . Right, the field potential for a single scalar field. Here by vacuum we mean a spatially 
homogeneous configuration where the field is resting at (/)i, (p 2 or 


an infinite universe, the situation is trickier. The chance that the field tunnels through or jumps 
over the barrier everywhere at once is zero. From the WKB perspective the tunneling exponent is 
infinite and thermal jumps would need an infinite energy, so according to Boltzmann the transition 
rate ought to be zero. What happens in reality for quantum tunneling is that a bubble of the new 
phase is formed and expands. This is the familiar prescription for the first order phase transitions 
we see in the everyday life. When we boil water, it does not evaporate all at once. Instead small 
bubbles form near nucleation centers and then expand until the whole liquid phase disappears. 
The same thing will happen in flat space vacuum decay and after a long enough time no trace 
of the old vacuum is left. When we add gravity to the picture, the situation gets more obscure 
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and meanwhile more interesting. Because of existence of a horizon there is a nonzero Hawking 
temperature which makes thermal tunneling feasible. The old vacuum loses some of its volume to 
the decay and meanwhile it gains more volume because of the cosmic expansion. Depending on the 
rate of the bubble nucleation and the cosmic expansion, we may have different scenarios[I]. If the 
tnnneling rate is very fast, the transition will complete in a hnite time and the whole space converts 
into the new truer vacuum. On the other hand, if the expansion rate is much faster than the decay 
rate, more and more of the true vacuum bubbles nucleate and chip away parts of the false vacuum. 
These bubbles grow and there may be new bubbles created inside each of them. Meanwhile the false 
vacuum itself grows much faster and the overall picture is a universe filled with the false vacnum 
and with a large number of isolated true vacuum bubbles in the false vacuum. The transition never 
completes. Here there are a variety of interesting problems about “what fraction” of the space is in 
the new or old vacnum. Since the space is infinite, defining a measure on this set is an immediate 
question necessary before giving any probabilistic or anthropic description of nature. We do not 
know the answer to this question yet[2HS]- 

In this thesis I will cover several topics regarding vacunm decay. 

In Chapter I will review the basics of tnnneling and describe the formalism needed for the 
following chapters. I start by describing the meaning of the decay in a simple quantnm mechanical 
system and then generalize this concept to a field theory. There are two approaches available on 
the market to deal with the problem of vacuum decay. One is based on the WKB approximation. 
The other uses a Euclidean path integral approach. I will briefly review both and then describe 
thermal tunneling. 

The tunneling rate is determined by the Euclidean action of the solution that carries the decay. 
The higher the action of this solution, the less probable the decay described by it. We know that 
generically the highly symmetric solutions have lower action. It was shown that in flat spacetime 
the solution with lowest action and therefore highest rate has an 0(4) symmetry [6] . However there 
is no snch proof for field theories in curved spacetime. The notion that the higher the symmetry, the 
lower the action seems to be wrong in some cases. For example the solution with the highest possible 
spacetime symmetry, the Hawking-Moss solution which has an 0(5) symmetry, has a higher action 
than the 0(4)-symmetric solutions in most cases. In Chapter]^ I will describe solutions with the 
next highest symmetry, solutions periodic in Euclidean time and spherically symmetric in spatial 
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directions. These solutions have an 0(3) x 0(2) symmetry. I will obtain such solutions numerically 
and analytically. I show that for a wide range of potentials these solutions have higher action than 
the 0(4) symmetric solutions and therefore are subdominant. 

The usual scenario for vacuum decay is bubble nucleation in a scalar field cj) whose potential 
has more than one local minimum. In Chapter]^ I study bubble nucleation in a spatial vector field 
theory. This vector field has different speeds of sound in the longitudinal and transverse directions. 
I show that in this case the bubble does not have a spherical symmetry. This work is motivated 
by the phenomenon of A-B phase transition in the liquid ^He. I show that if the ratio of the 
longitudinal to transverse speeds of sound goes beyond a threshold, the critical bubble develops a 
kink and the flat domain wall breaks into zigzag segments. 

In the last few years anthropic principles have become very popular as a paradigm for explaining 
the cosmological constant problem. The general belief is that it is possible to have an enormous 
number of vacua in compactifications of string theory. It is assumed that there are huge number of 
vacua, perhaps (500 is a number quoted frequently, it could be several thousands as well). 

If there are so many vacua, we should not be very surprised to see a cosmological constant which is 
120 orders of magnitudes smaller than the Planck mass. Having 10^^® uniformly distributed vacua, 
we will get something like 10^®^ vacua with cosmological constants as small as the observed value. 
But it is not proven that all of these vacua are stable. Recently there has been work that shows that 
many of these “vacua” are indeed saddle points rather than local minima. Of course even if “most” 
of these vacua were unstable, let’s say 99.99% of them, still we would end up with 10^®® stable ones. 
However we need to be worried about the lifetime of these vacua. Our Universe is an old one with a 
long lifetime. On the other hand, using string theories with many moduli means that the effective 
potential depends on many variables. In a potential with many variables, there are many different 
directions for tunneling. Some of these directions are crossing low barrier heights. This makes it 
possible for the metastable vacua to tunnel easily. I present the analysis of the lifetime of the vacua 
with large number of fields in Chapter Because performing the exact calculations in the string 
theory context is a nightmare, we used a simple model of a theory of N scalar fields. We looked 
at the behavior of their potential near the vacua. Our conclusion is that the vacua become highly 
unstable with a large number of fields. In fact we show that the lifetime drops like an exponential 
of a power of N. Therefore if our analysis applies to the landscape of string theory, we cannot find 
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enough diversity of vacua to describe the smallness of the cosmological constant in a natural way. 

Our analysis in Chapter is based on a model of N scalar fields with a random potential. 
However we need to be very conservative in applying it to the landscape of the string theory. There 
are many points that need to be clarified and taken into account before such an application. Among 
them the non-canonical kinetic terms, the possibility of brane nucleation and keeping higher order 
terms can be listed. In Chapter I present a simple model of compactification of the Einstein- 
Maxwell model which presents many of the features we expect from the landscape of string theory. 
I describe the landscape of the vacua of this theory which includes N moduli which are the radii of 
the extra-dimensional spaces. I present some of the nice results we obtained about this landscape. 

Finally I conclude in Chapter 
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Chapter 2 

Vacuum decay, background 

2.1 Vacuum decay in one-dimensional quantum mechanics 

Although the minima of a potential are considered to be stable states in classical mechanics, the 
story changes drastically when we deal with quantum or thermal physics. In quantum mechanics, 
the phenomenon of barrier penetration allows escaping from a classically forbidden region. In 
thermal systems, excitation through absorption of energy from a heat bath makes this jump possible. 
Figure |2.1| shows a simple potential for a particle moving in one dimension. If the particle is 
originally at the right minimum, it can penetrate through the barrier and go to the left side and 
emerge at D. The exact treatment of such a system is utterly difficult and we have to use the WKB 
approximations to solve it. The tunneling probability is given by 

T ^ ^ ^ 2 . 1 ) 

The outcome of the tunneling is a particle at rest emerging at D. Similarly, thermal fluctuations 
may render the right minimum unstable by boosting the particle to C with a probability 

V(C)-V(A) 

P^e . ( 2 . 2 ) 

The outcome will be a particle emerging at rest at C that may then classically roll down to the 
more stable minimum B. 
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Figure 2.1: An example of a potential with a metastable vacuum at A and a stable vacuum at B. 


2.2 Imaginary energy states and decay 

The time evolution of a particle with an imaginary energy is not unitary. The probability of its 
survival decays exponentially over time. It is unstable. This is similar to the situation we have in 


hand in the decay of a false vacuum. Let’s start with a simpler system shown in Fig. 2.2 If V 2 
was inhnite, the left side would admit its own ground state with a wave function 


2 TTX 

'®“ = VT“T' 


(2.3) 



0 b b+w 


Figure 2.2: An inhnite well sectioned into two pieces, left a false vacuum of width b and potential 
Vi separated by a wall of width w and potential V 2 from the rest of the well. 

However the hniteness of the barrier makes this state an approximately stationary one. Let’s 


choose the values of the potential in units of 


Vi = Vi 




2mb‘^ 


(2.4) 
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We can use the Schrodinger equation to evolve the approximate ground state [Eq. (2.3)] in time and 
calculate the chance of finding the particle on the left side of the wall, i.e x < b. These probabilities 


are shown in Fig. 2.3 for the cases where Vi > V 3 (left) and Vi < V 3 (right). As is seen, the 
former is unstable and the latter is quite stable. In addition, for the left case we see that the 
chance of survival of the particle drops exponentially with time (at least for a short period of time), 
which is reminiscent of the decay of unstable particles with imaginary energies. Of course, this 
exponential decay cannot continue forever in this simple system because after a while the particle 
starts oscillating between the left and right sides of the well. We can see this in Fig. |2.3| where 
the linear behavior on the log plot gets obscured after some time has passed. To understand this 
better, let’s repeat the same simulation for fixed potential parameters and different values of a, the 
location of the right part of the well. 


Probability 
1.00 k 

0.70 

0.50 


0.20 

0.15 


Probability 

1.00 r 


0.50 

0.20 

0.10 

0.05 

0.02 


Time 


Time 


Figure 2.3: Chance of finding the particle at the left of the barrier as a function of time. Left is for 
the case vi = 4 ,V 2 = 40, V 3 = 0, Right for the case vi = — 8 , V 2 = 40 and V 3 = 0. Time is measured 
in units of 2mb‘^h~^ and the vertical axis has logarithmic scale. In both cases, a = 20, w = 0.2 and 
6 = 1 . 


The survival chance for different values of a is shown in Fig. 2.4 Interestingly enough, the 
exponentially decaying part does not depend on a and its slope is the same for all of them. But 
the bigger a is, the longer the time that the exponential decay approximation is correct. Therefore 
the imaginary part of the energy only depends on the barrier and the state on the left side and is 
independent of whatever happens on the right side of the wall. But making more space available on 
the right of the barrier results in a higher density of states which in turn increases the time it takes 


^To be more accurate, we need to make a comparison between the ground state energies on the left and right sides 
of the barrier. This may or may not be the same as the comparison of the potentials. 
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Figure 2.4: Probability of finding the particle on the left side of the barrier as a function of time. 
The potential is described by vi = A, V 2 = 40, us = 0, 5 = 1 and w = 0.2. Time is measured in 
units of 2mb^h~^ and the vertical axis has a logarithmic scale. The red filled circles, blue squares 
and black empty circles correspond to a = 40, 20 and 5. 


for the particle to return to its initial state. For the case that we will be considering the most, 
quantum field theory, the number of states available is inhnite and therefore the decay process 
will be a real decay with no return to the original state. This means that we are dealing with 
states with an imaginary part of energy. But wait! The Hamiltonian in the form /2m + F is a 
Hermitian operator and its spectrum can only admit real energy states. How can it accommodate 
an imaginary energy? To understand this, let’s look at the quantum anharmonic oscillator. The 
Hamiltonian for this system is 

with the associated boundary conditions 


lim 4'(x) = 0 . (2-6) 

|ir|^oo 


For simplicity, we set m = i,a; = l,^=l. This potential for different values of A is shown in Fig. 


2.5 . If there are unstable states for local but not global minima, we should expect to get states 


with imaginary energies for negative A. 
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V(x) 


V{x) 




V(x) 



Figure 2.5: The anharmonic oscillator potentials defined in Eq. (2.5) for different values of A. Left, 
middle and right graphs correspond to A = 1, 0, —1 . 


Let’s dehne E{X) to be the energy eigenvalues of the Hamiltonian in Eq. (2.5) constrained to 


the boundary conditions of Eq. (2.6). Using Symanzik transformations, we can show that |E(A)| 


grows like |A|^/^ for very large A. Loeffel and Martin [7] showed that E{X) is analytic in the complex 
A plane except for the cut and the dashed region in Fig. |2.6[ Loeffel, Martin, Simon and Wightman 
[8] showed that the LC’th energy level can be written as an asymptotic series 


E^(A) = iL+- + ^HfA^ 


(2.7) 


n=l 


It is apparent from Fig. 2.6 that moving from positive to negative A, which corresponds to moving 
from a stable to an unstable potential requires landing on a branch cut. We should expect that 
something nontrivial happens here. This makes the physical Hamiltonian for negative A the limit 
of a non-Hermitian operator and the branch cuts and discontinuities make the transition to the 
real potential non-smooth. This causes the energy levels to be complex numbers. Bender and Wu 
in a series of papers [9Hn] showed that E{X) acquires an imaginary energy part and were able to 
calculate this imaginary energy for small negative A = — e to be 


ImE^(A) = 


4^+5 


-K-- 
ze 3ee 2 


( 2 . 8 ) 


K\V^ 

This number is very tiny for small e, which corresponds to a tall wall. For a general one-dimensional 


potential the exponential part of the energy in Eq. (2.8) can also be calculated using the standard 


WKB methods for hnding the transmission coefficients. 
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Figure 2.6: Analyticity domain of E{X) in the complex A plane. The inner product is not defined 
for all vectors and it has a nontrivial analytic structure. Therefore these eigenvalues evade the 
usual proof for the reality of the eigenvalues of the Hermitian operators 

2.3 Vacuum decay in multi-dimensional quantum mechanics 

Vacuum decay in many dimensions has slightly different qualitative features. Let’s look at the 
two-dimensional potential in Fig. |2.7[ When the particle tunnels out of the central minimum, it 
can go to any of the points on the (blue) line which have the same energy as the central minimum 
and respect the conservation of energy. After the tunneling, the particle can emerge with zero 
kinetic energy at any of these points. In addition, any path that goes from the minimum to the 
exit point is an acceptable path. These paths do not contribute equally and they are exponentially 
suppressed. When we add them, we only need to take into account the ones which have the least 
damping (resistance). 

In a series of papers. Banks, Bender and Wu (BBW) |12l [T3] generalized the one-dimensional 
WKB methods to N-dimensional systems of coupled harmonic oscillators. Their idea is to break 
down the tunneling problem to a one-dimensional WKB approximation along the classical path to 
get the exponential part of the tunneling rate. The prefactor is calculated using a small tube that 
encompasses this one-dimensional path. The amplitude for a particle to take a path is proportional 
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Figure 2.7: Possible exit points in a multi-dimensional potential. 


iS 


to e'ft , where S is the action along the path. This exponent has an imaginary part for the case 
of tunneling and therefore the probability is exponentially small. Therefore to get the largest 
contribution we can minimize the action along the path. Using the Maupertuis principle, the 
problem of minimizing the action becomes equivalent to the variational problem 


5 j ds{V -EY^^ = 0 . 

Here s is the path length along the trajectory. The corresponding variational equations are 

dxi dxn dV dV 

2{v -E)-^ + -7^^ = 7 ^ . 

^ ds^ ds ^ ds dxj dxi 


(2.9) 


( 2 . 10 ) 


If we had chosen time as a parameter, the sign of the potential gradient on the right side would 
be the opposite of what we would expect. In ordinary cases, we expect Xi = —diV. This means 
that the potential should be treated as an upside down potential or, equivalently, we should solve 
the problem in a Euclidean spacetime. We should have seen this if we were more careful in the 
context of barrier penetration. There the momentum became purely imaginary and to get a normal 
momentum, we should have made a transformation t —)> ir. This is again changing to Euclidean 
spacetime. Solving these equations, we get a set of MPEPs, the most probable escape paths, which 
carry most of the decay. Then we look at the variations of the wave function in narrow ’’tubes” 
around these trajectories. Here we should distinguish between the two levels of approximation 
involved. The first is treating the amplitude of the wave function as a constant and only looking 
at the damping phase of it along the trajectory. This is similar to geometrical optics and zeroth 
order WKB and it gives the exponential damping. Then we look at the variation of the amplitude 
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around this trajectory which is the first order WKB that gives the non-exponential part of the 
decay rate. Using this method, BBW could calculate the imaginary part of the energy of coupled 
anharmonic oscillators with straight MPEPs could give a recipe for a general MPEP. However the 
calculation is very difficult and generalization to many dimensions is a challenge. 

As mentioned earlier there is another type of decay, thermal tunneling, which is possible in the 
existence of a heat bath. In thermal tunneling, the path of least resistance will be a path which has 


the lowest barrier height. As is clear from Fig. 2.8, the point which has the lowest barrier height 
is a stationary point where the potential is increasing in all directions except for the one that leads 
to the neighboring minimum. The particle jumps to this saddle point and emerges there. It may 
then classically roll down to the other minimum. 



Figure 2.8: The dominant conhguration for thermal tunneling in a 2D potential is jumping to a 
saddle point with one negative eigenvalue of the Hessian matrix. 


2.4 Vacuum decay in field theory 

Vacuum decay in field theory was studied by Langer [13] in the context of statistical physics of 
saturation of droplets and also by Kobzarev, Okun and Voloshin [15] . However the modern accepted 
view appeared in two seminal works of Coleman m and Callan and Coleman m- In this section. 
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we follow the approach of the latter two works. The decay rate is given by T = Ae~^^^, where B is 
shown to be the Euclidean action of the configuration that leads to tunneling. In first part of this 
section, a WKB generalization to field theory is presented and the exponential part of the tunneling 
is calculated. In the second part, we use a path integral approach to calculate the prefactor of the 
tunneling rate. 


2.4.1 WKB approach 


Let’s start from a scalar field theory with a Lagrangian 

C = ^d^cpd^4>-Vicp), ( 2 . 11 ) 

where V is shown in Fig. |2.9| and we have shifted the potential so that the false vacuum has a zero 
potential. We want to find a configuration in which the field rests in its false vacuum at a long 
(Euclidean) time ago and a bubble emerges at rest at time r = 0. These correspond to 



( 2 . 12 ) 

(2.13) 


Figure 2.9: True and false vacua in a scalar field potential. For convenience, we set V{(j)F) = 0 . 


To get the dominant tunneling path, we need to find the stationary points of the Euclidean 
action 

SE = Jdri>x K^) ^'"*'** ■ 

Therefore we have to solve the Euclidean equations of motion 
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Because this equation is invariant under time reversal, the first boundary condition in Eq. (2.12) 
can be consistently written as 


lim 4>{t,x) = 4)f 

T ^±00 


(2.16) 


This means that we look for field configurations which start from the false vacuum, evolve to the 
field configuration for tunneling at r = 0 and then bounce back to the false vacuum. These are 
called bounce solutions. To get a finite Euclidean action and therefore a nonzero tunneling rate, we 
need to make sure that at any moment of time. The field rests at its false vacuum for large spatial 
distances 

lim = cj)F ■ (2-17) 

|:r|^oo 

The center of the bubble can appear anywhere in space and, by shifting the time origin, we can 
arrange the tunneling to occur at any desired Euclidean time. This means that the nucleation 
nucleation probability obtained for a specific bubble is in fact the nucleation probability per unit 
time per unit volume and the total nucleation rate takes a factor of spacetime volume. Coleman 
proved that the configuration with the least Euclidean action has an 0(4) spherical symmetry. 
Let’s assume that the field is only a function of the four-dimensional distance p = (r^ + |xp) ' . 


Equation (2.15) simplifies to 


(Pcj) 3 d(j) dV 
dp'^ p dp d(j) 


(2.18) 


This equation resembles the motion of a classical particle moving in an upside down potential 


shown in Eig. 2.10 Boundary conditions in Eq. (2.16) correspond to a motion starting at time 

(p 



Eigure 2.10: The motion in Euclidean time is equivalent to motion in an upside down potential . 


p = 0 at somewhere close to the true vacuum cpF with zero speed (^ = 0) and landing on the false 
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vacuum (pp at p = oo. This motion is affected by the presence of a time-dependent friction force 
proportional to velocity (^^)- If </>(0) is very far from pp, it will not have enough energy to reach 
pp. If it starts too close to pp, it will stay there for a long time. This suppresses the friction term 
so due to the difference between the potentials and lack of enough friction, it will pass pp m a finite 
time. Therefore there must be some point in between for which the particle does not overshoot or 
undershoot and reaches pp va. aa infinite time. This proves that the Euclidean equations of motion 


always admit a solution. Except for special cases, Eq. (2.18) does not have a closed form solution. 


However if the potential difference between the true and false vacua e = V{pp) — V{Pt) is small, 
we can solve this equation using a thin-wall approximation. If e is small, in order for the field to 
make it to pp, the effect of friction should be very tiny. Therefore the field should remain close 
to pT until p becomes equal to some large radius R and therefore friction gets suppressed. This is 
only possible if it starts very close to pp- Then for p > R we can neglect the friction term and Eq. 


(2.18) simplifies to 


d?p 

dp^ 


dV 

dp 


This equation has a solution of the form 


P = 


f 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


/o V2^V{P) - V{Pp) ' 

The Euclidean action for the spherically symmetric solution is 

1*00 1 / j \ 2 

Sp = 2 p^ 

Jo 

The decay exponent B is the difference between the Euclidean action of the bounce and the pure 
false vacuum solution. We can break this integral into three parts. The first part is where the 
field stays very close to the true vacuum {p < R) and we can neglect the derivative terms. After 




subtraction of the pure false vacuum action we get 

rR 


Bi = SiE — Sip = 27r 


dp p-^ [V{Pp)-V{Pp)] = -^-p^R^e. 


to 


( 2 . 22 ) 


In a short transition region, which is the wall region, we can treat p to be constant and we get 


B2 = S2E — S2P = 2 p 


2 l^'dpp^ 

>0 




2 p‘^R^ 


rR 


dp 


l{%] 


= p^R^a . 


(2.23) 
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Here d is a number that depends on the explicit form of the potential and represents a surface 
tension. For R < p, where the field get very close to its false vacuum value (in fact exponentially 
close for any potential), there is no difference between the false vacuum action and the Euclidean 
action. Therefore the net contribution to the decay exponent is only due to the wall and the interior 
region of the bubble and is given by 

B = + 27r‘^R^a . (2.24) 


We are left with one free parameter R which determines for how long the field stays close to its 
true vacuum value. We can vary this number to get the largest possible decay rate. Doing so we 
get 


0 = ^ = -27r^Rh + 67r^R^a 
uK 


(2.25) 


This leads to 


R 

B 


3a 

277rV^ 


(2.26) 

(2.27) 


It is important to notice that the transition happens in the configuration space of the field and 
not through a straight path through the barrier. The term potential barrier may be misleading 
since the spatial gradient term is effectively a potential if we think about the field as points on a 


lattice. We will get back to this point in more detail in Sec. 2.5 The family of field configuration 


that leads to tunneling is shown in Fig. 2.11 where a bubble of true vacuum pops out of the false 
vacuum and reaches to a maximum size at r = 0 and again disappears when r gets very large . 


Stacking all of these spatial slices together, gives the picture in Fig. 2.12 


2.4.2 Growth of the bubble. 

Until now we have described the formation of a bubble through quantum tunneling. Inside a space 
filled with a false vacuum, a bubble of true vacuum pops out and the field configuration is given by 


(f>{x, t = 0) = (/>(x, r = 0) 

§^(p{x,t = 0 ) = 0 . 


(2.28) 
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Figure 2.11: Bubble nucleation as seen in Euclidean time from top left to the bottom right. In 
the beginning the whole space is filled with the false vacuum (shaded light blue region) and then a 
small bubble of true vacuum (green region inside the circles) with a wall (red region which is the 
boundary of the circle ) appears and grows. The top right configuration shows the configuration at 
the moment of the bubble nucleation. For large Euclidean times, it bounces back to a space filled 
with the false vacuum. 



Eigure 2.12: Another view of the bounce. In this picture the Euclidean time progresses upwards, 
the false vacuum is pictured as a white background and the true vacuum is in red. Any horizontal 
line cuts the bounce at a spatial hypersurface. In the far past (the lower thick dashed line), the field 
is almost at its false vacuum, and the cut through the center of the bounce (upper thick dashed 
line) shows the field configuration right after the nucleation of the bubble. 
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If we assume that the size of the bubble is not too small, we can describe the subsequent evolution 
of the field by its classical (Lorentzian) equations of motion 






(2.29) 


Comparison between Eq. (2.15) and Eq. (2.29) shows that the solutions of the Lorentzian 


equations of motion are the analytic continuations of their Euclidean counterparts and can be 
obtained by replacing t ^ it in the Euclidean solutions. Therefore 


t) = 4 >{x, ir) = 4>{p) = 4 >{\x\‘^ - t^) 


(2.30) 


As seen from Eq. (2.30), the 0(4) symmetry of the Euclidean solution translates into the 
0(3,1) symmetry for the consequent expansion of the bubble. Again we can get intuition about 
the solution by going to the thin-wall limit. In this approximation, the location of the wall is at 


| f |2 -^2 ^ 


(2.31) 


The bubble radius R is determined by the potential as described in the previous subsection. It 
should be of the same order as the energy scales of the scalar field and therefore a relatively short 
length compared to macroscopic lengths. This means that immediately after the nucleation of the 
bubble, the wall moves almost at the speed of light and starts eating away more and more of the 
false vacuum. The wall’s Lorentz factor 7 = (1 — from Eq. (2.31) is 


7 = 


R 


(2.32) 


During the conversion of the false vacuum into the true vacuum, some energy is released. This 
energy is spent on accelerating the wall. It is easy to show that in the thin-wall limit, all the energy 
gained from the transition is exactly converted into the kinetic energy of the wall. This means that 
when the bubble passes a point, there is no ripple or radiation left behind it and only the true 
vacuum region at rest is created. 


2.4.3 Path integral approach 

Coleman and Callan m used a path integral approach in Euclidean spacetime to calculate the 
tunneling exponent and prefactor for the decay of metastable vacua. In this section we follow 
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their calculation. For simplicity let’s start again from a one-dimensional quantum mechanics in 
imaginary time and then generalize it to field theories. The amplitude for a particle to move from 
Xi at time —Tf2 to Xf at time T/2 is 

{xf\e-^^/^\xi) = N j[dx\e-^^/^ , (2.33) 

where in the path integral side the sum is over all the paths that satisfy the boundary conditions 
x{—T/2) = Xi and x{T /2) = Xf and is a normalization factor. The left hand side of this equation 
has a simple expansion in terms of the energy eigenstates of the Hamiltonian: 

{xf\e~^'^/^\xi) = . (2.34) 

n 

In the limit where T —)• oo, the ground state contribution dominates and 

{xf\e~^'^/^\xi) = e~^°'^/^{xf\t)){t)\xi) . (2.35) 


This gives a simple expression for the ground state energy Eq. In this section we evaluate the 


ground state energy Eq by calculating the right side of Eq. (2.33) using standard path integral 
methods in the semiclassical approximation. We will see that indeed it has an imaginary part if 
the energy eigenstate is localized around a metastable vacuum. This imaginary part of the energy 
gives the decay (nucleation) rate. 

In the limit h —)• 0, the main contribution to the path integral comes from regions close to paths 
x{t), the stationary points of the Euclidean action where 


5Se d‘^x , 


(2.36) 


This is the equation of motion of a classical particle moving in an upside down potential —V. We 
use the method of steepest descent to evaluate the path integral in a neighborhood of the classical 
path. But before proceeding further, let’s look at two different cases where Xi = xj = Xmin) one 


for a global minimum and one for a local minimum. These two cases are shown in Fig. 2.13 If xi 


is a true (global) minimum, the only solution to Eq. (2.36) is a constant solution in time. But if 
it is a local minimum and not a global one, there is another solution which starts from xi at time 
—oo and then moves to the left and after an infinite time returns to xi. For obvious reasons this 


solution is called a bounce and is shown in Fig. 2.14 At the semiclassical level, we can evaluate 
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Figure 2.13: The Euclidean equations of motion have drastically different behaviors for paths 
starting from the true and false vacua. On the top, the two potentials with minima at xi and at 
the bottom their upside-down versions which are more relevant in Euclidean spacetimes. In the 
left potential, xi is located at a true vacuum and there is no solution for the Enclidean equations 
of motion except for x{t) = xi. But for the left, in addition to this constant solution, there is 
another solntion that starts from xi, moves to X 3 and then returns to xi after an infinite time. 
This solution is called a bounce. 



X 



Eigure 2.14: Left is a single bounce centered at to and right a multi-bounce configuration. Both of 
these are approximate stationary points of the action and will contribute to the path integral in the 
semiclassical level. (The left bounce would be an exact stationary point if the bounce was centered 
at t = 0 ). 
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the path integral by looking at paths close to xi. This can be done easily by expanding these paths 
in terms of a complete set of orthonormal functions Xn{t) which satisfy the boundary conditions 
Xn{±T/2) = 0: 


x{t) = x{t) + ^ CnXnit) , 
n 

fT/2 


l-T/2 


dt Xji(t^ XjYiiUj — Sj, 


(2.37) 


The integration over all paths now simplifies to an integration over the coefficients Cn- We define 
the measure to be 

[dx] = J^(27rh)“^/^(ic„ . (2.38) 

n 

We choose x„(t) to be the eigenstates of the second variation of the Euclidean action 

d?Xr 


"■ + V"{x)xn = XnXn ■ 


dt'^ 


(2.39) 


Now the path integral simplifies to Gaussian integrals over c^’s 




[det(-a2 + E"(x))] ^'\l + 0{h)] 


(2.40) 


We made an important assumption here that the A^’s are all positive. Otherwise the integrals would 


diverge. For the case of a true vacuum, shown in the left part of Fig. 2.13 the only contribution 
comes from the path x = 0 and therefore V"{x) = V''{xi) = is a constant number. A simple 
calculation (for example [18]) leads to 


N 


[det (-52 + F"(x))]-'/' = e-'^/2. 


(2.41) 


The energy of the ground state has shifted by the familiar zero energy of a harmonic oscillator. 


The situation is different for the false vacuum depicted in the right panel of Fig. 2.13 Here we have 
to take into account the contribution from the other stationary point of the action, the bouncej^ 
For the bounce solution, we have x(±T/2) = xi and at some time to we have x{to) = X 3 and 


^ The bounce solution is only an exact stationary point of the action in the limit T —>■ 00 and we should justify 
using this approximate stationary point in the sum over all paths. However there is a more accurate treatment of this 
issue in page 275 of [18] by comparing the approximate stationary points with having stationary points at infinity. 
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x{to) = 0. This is the turning point. First we calculate the action for a single bounce. In Euclidean 
time, the energy is a constant of motion 


E = 


1 f dx'' 


2 \dt 


- V{x) . 


(2.42) 


We shifted the potential so that V{xi) = 0. This makes E vanish for the bounce. Let’s call the 
Euclidean action of this bounce B. 


B= f dt 

J —OO 


1 / dx^ 






dxy/2V(x) . 


(2.43) 


’X3 


We also need to sum over multi-bounce configurations. If a multi-bounce configuration is composed 
of n separate bounces, its action is nB and the det[—Of -|- V''{x)\ is the product of the determinant 
of n single bounces which are separated by large time intervals and the determinant is 


\TThJ 




(2.44) 


Here K is a, factor chosen so that this expression is correct for a single bounce. For multi-bounce 
solutions, the turning points can take different values and we have to sum over the location of all 
of them. This sum contributes a factor 


rT/2 

l-T/2 


Cin — 


dti 


l-T/2 


dt 2 


n-1 1 

dtn = -.T^ 


l-T/2 


n\ 


Now we can sum over all the contributions from bounces: 


/cux 



1/2 / 
1 exp ( 

\TTh) 

n! ' 

Wh/ 


n=0 

Therefore the energy of the ground state becomes 

Eo={^- [1 + 0{h)] . 


(2.45) 


(2.46) 


(2.47) 


Please notice that the second term in the parentheses is much smaller than h?, so in general it does 
not make sense to keep it here. However, as we will see shortly, K is imaginary and this term is in 
fact the first nonzero contribution to the imaginary part of the energy. Before calculating K, we 


should clarify some points about the eigenvalues of Eq. (2.39). In order to perform the Gaussian 
integrals in the path integral, we needed to assume that all A^’s were positive numbers. But because 
of the time invariance of the equations of motion, there is one zero eigenvalue, corresponding to 


xi 


= 


dx 

dt 


(2.48) 







CHAPTER 2. VACUUM DECAY, BACKGROUND 


25 


The factor B is introduced here to satisfy the convention in Eq. (2.37). Integrating over this zero 


eigenvalue yields a factor proportional to T. Fortunately, we already have taken care of it when we 


integrated over the location of the turning point to- The eigenfunction xi is shown in Fig. 2.15 


Because the zero eigenfunction has a node, there must be a lower eigenvalue which is negative. This 
is worrisome, because now the path integral diverges. This should not be very surprising, because 
from the beginning we were trying to calculate the energy of a state localized near the false vacuum. 
We know that such a state should not be stable and therefore it should not be part of the spectrum 
of the Hamiltonian. The correct way to treat this energy is by analytic continuation. To make the 



Figure 2.15: The zero eigenfunction of the second variation of the action. This eigenfunction has a 
node so there must be another lower (negative) eigenvalue in the spectrum of these variations. 


computation as simple as possible, we restrict ourselves to a subspace of paths parametrized by a 


real parameter z. These paths are shown in Fig. 2.16 The z = 0 path is the constant path that 
stays at the false vacuum and therefore has zero action. The z = 1 path is the bounce and we chose 
the paths in such a way that the tangent to it is the negative mode. For large z, where the path 
remains in the negative potential region near X 2 for a long time, the action can get arbitrarily large 
and negative. If we evaluate the path integral over z, which is 


J = {2Trh)-^ / 


(2.49) 


we will end up with a badly divergent integral for z S> 1. To remedy this divergence we use an 


analytic continuation of the potential. If instead of the potential on the right side of Fig. 2.13 


we had used the potential on the left side, the action as a function of z would look like Fig. 2.17 


where the integral is convergent. Now let’s devise an analytic continuation of the potential on the 


left side of Fig. 2.13 which maps it to the potential on the right side. We choose this continuation 
in such a way that along the real z axis, we recover the right potential and on the upper complex 
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X s 




Figure 2.16: Left, a class of paths parametrized by z and right, their Euclidean actions as a function 
of z. The z = 0 path corresponds to the constant path x{t) = xi and has zero action. The z = 1 
path is the bounce. We chose the paths in a way that the tangent to the z = 1 path is the negative 
mode and therefore we are sampling the divergent region of the path integral. 


z plane we continue to the potential on the left. To take the z integral, we can follow a contour 
which is shown in Fig. |2.18[ This contour extends from —oo along the positive z axis to the saddle 
point at z = 1 and distorts to the upper half plane. The main contribution comes from the region 
near z = 1 and we can use the steepest descent method to evaluate the path integral: 


1+200 1 

-l/2h[2Sil)+S"{lKz-lY] ^ ^ ^-S(l)/h 

2vW^ 

Because the integration took a path only half-way around the saddle point, a factor of 1/2 appeared 
here. By generalizing this analysis to the path integral over the whole function space, we get 


(2.50) 


ImJ = {2'Kh) ^Im 


S 




Figure 2.17: Euclidean action as a function of z for the potential shown on the left graph in Fig. 


2.13 This function is positive and the path integral is well-defined. To get the correct result for 

we analytically continue it to the potential on the left. 


the potential on the right side of Fig. 


2.13 


Im 




one bounce 


2 



1/2 


T |det' [-+ + F"(a 


1 - 1/2 


(2.51) 
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Figure 2.18: The contour of integration chosen for the path integral along z which continues from 
—oo on the real axis and extends to the saddle point at z = 1 where the bounce is located. Then 
it distorts to the upper half plane where the Euclidean action is again positive. 


where det' means evaluating the determinant ignoring the zero mode, which we have already taken 


care of by integrating over the turning points. By comparing this with Eq. (2.44), we find that the 
imaginary part of K is 

- 1/2 


ImK = - 
2 


The decay rate per unit of time is 


T = -2Im^ = 
h 


/ B y/2 

det' [—+ V"{x)\ 

\2Trh ) 

det [-5| + a;2] 


(2.52) 



det' [—+ E"(x)] 

\2TTh J 

det [—+ a;2] 


- 1/2 


(2.53) 


The generalization to field theory is very straightforward. We need to notice that there are four 
modes, corresponding to four translations in spacetime. This causes four factors of 


zero 


in Eq. (2.53). The decay (nucleation) rate in this case is per unit volume and we get 


T 

V 






.-B/h 


det' [—□ + x))] 


det [-□ + E"(</>fv)] 


- 1/2 


(2.54) 


where 4>{t, x) is the bounce solution. In the case of a scalar field theory in four-dimensional spacetime 
with 0(4) symmetric solution, Coleman proved that there is a single negative mode. It was crucial 
for our argument that there is a single negative mode. If there are an even number of negative 
modes, there will be no imaginary part contribution to the ground state energy and it cannot 
predict a decay. Also a higher number of negative modes for a bounce is usually an indication 
that there is a nearby bounce with lower Euclidean action which is going to be the dominant path. 
Therefore the bounce with more than one negative mode will be a subdominant decay mode. 
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2.5 Tunneling at finite temperature 


In this section we consider the effect of a nonzero temperature on tunneling. We assume that the 
system is in equilibrium with a thermal bath at a temperature T. One of the interesting cases that 
we will encounter in the future is the nonzero temperature due to the existence of a horizon. 

To treat the problem properly we have to use the finite temperature effective potential I4ff (</>) T). 
To simplify the notation we drop T in the text, but have it implicitly in mind. Also we need to 
take into account the new ways for tunneling at nonzero temperature |19[ [20] . Again our starting 
point is a one-dimensional quantum mechanical system and then we will look at the generalization 
to field theory. There are two new possibilities for tunneling. We can call them thermal tunneling 


and thermally assisted tunneling. These are shown in Fig. 2.19, The horizontal arrow shows the 
purely quantum tunneling that we covered in the previous sections. The dashed arrows show a 
thermally assisted tunneling in which the particle jumps thermally to point 2 which has energy Et 
and then quantum mechanically tunnels to point 4. The dot-dashed arrow shows a purely thermal 
jump from point 1 to 3, which is a stationary point with one negative eigenvalue. Purely thermal 



Figure 2.19: At a nonzero temperature we can distinguish three different paths for tunneling. 
The horizontal arrow shows a purely quantum mechanical tunneling. The dashed arrows show 
a thermally assisted tunneling in which the particle first jumps thermally from 1 to 2 and then 
quantum mechanically tunnels to 4. The dotted dashed arrows show a purely thermal tunneling 
in which the particle jumps to a stationary point with one negative eigenvalue. In all cases, after 
tunneling to the right of the barrier, the particle evolves classically towards the lower minimum. 


tunneling is the easiest one to explain. The particle jumps from xi to X 3 and then rolls classically 
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down to xq. In more than one dimension, the point X 3 is replaced by a stationary point of the 
potential which has a single negative mode. This point is the lowest point on a ridge that surrounds 
the false vacuum. The rate for this tunneling is given by 

Ttunn ~ ^ (2.55) 


To calculate the tunneling rate, we do not need a path in configuration space. We only need 
an escape point. The easiest way to generalize this concept to a field theory is using the energy 
functional 


E[m] 



1 

2 




(2.56) 


The potential energy part is 


U[^{x)] 



^(Vcpf + VicP) , 


(2.57) 


which is different from V{4>). This is depicted in Fig. 2.20 


U[0(T)] 



Figure 2.20: The same as Fig. 


2.19 


but for a field theory. This time the horizontal axis rep¬ 


resents the (infinite-dimensional) field configuration space. Points 4>i{x) and 4>6{x) show the two 
homogenous conhgurations where the field is spatially uniform and rests in its vacua. <^ 3 (x) is a 
configuration which has only one negative eigenvalue and it is in fact a bubble of true vacuum 
inside the false vacuum. Its negative eigenvalue is along the expansion/contraction directions for 


the bubble. (j) 2 {x) shows the bounce configuration at r = 0 and r = /3 in Fig. 2.23 and (l) 4 {x) is the 
configuration after the nucleation for a thermally assisted bubble. 













CHAPTER 2. VACUUM DECAY, BACKGROUND 


30 


As mentioned earlier, the tunneling happens through an infinite-dimensional configuration space 
with a potential U[4>(x)] not V{(p). The correct tunneling configuration is a stationary point 
of U[(f){x)] in the configuration space with a single negative mode. This stationary point is a 
static solution of the full equations of motion. A solution 4>{x) to these equations is not spatially 
homogenous and in fact is a three-dimensional bubble of true vacuum which is separated from the 
false vacuum by a wall. This solution is time independent and it has an 0(3) x 0(2) symmetry 
where the 0(2) symmetry comes from the periodicity in the time direction. The negative mode 
corresponds to expansion/contraction of the bubble. The energy of the critical bubble in the thin- 
wall limit is given as the sum of the contributions from the surface tension of the wall and the 
energy density inside the bubble. 

47T 

E = 47ri?^cj-’ (2.58) 

o 

where e = 17((^fv) — U{(f)u). We can find the energy of the critical bubble by maximizing this 

E 

R 


Figure 2.21: A quantum bubble appears with zero energy, as expected from the conservation 
laws. A thermal bubble appears with energy absorbed from the heat bath and in fact its energy 
is a maximum with respect to the radius. A thermally assisted tunneling is between these two 
extremes. 



expression with respect to R. This situation is shown in Fig. 2.21 The bubble emerges with 
nonzero energy that it absorbs from the heat bath. It is clear from the picture that the radius of 
the three-dimensional thermal bubble is smaller than the radius of the four-dimensional bubble for 
quantum tunneling. We can think about this static solution as a cylinder in four dimensions where 
each cross section corresponds to the critical bubble. The shape of this bounce solution is given in 


Fig. 2.22 


Now we get back to the thermally assisted tunneling for a single particle. The probability for 
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T=/J' 

T = 0 



Figure 2.22: Bounces representing thermal tunneling are cylinders in four dimension with a period 
of /?. Each slice through this bounce gives a three-dimensional ball which is a bubble of true vacuum 
separated from the false vacuum by a wall (the solid vertical edges of the rectangle). 


thermal transition from 1 to 2 is given by the Boltzmann factor 


Pi^2 = , 


(2.59) 


B{Et'), the Euclidean action for the transition from 2 to 3, is given by 

B{E't) = 2 / dxy/V(x) — Et . (2.60) 

Jx2 

Therefore the total rate of the transition is an integral over Et 

Ftunn ~ f ^ cIEt _ (2.61) 

JE{^ 

The main contribution comes from the region close to the Et* that maximizes the integrand 
(minimizes the exponent). Minimizing the exponent with respect to Et leads to 


/3 = -2 


= 2 


d 

dEr Jx2 

dx 


dx-\/V{x) — Et = 2 


dx 


1X2 


^V{x) - Et 


rX4 

J X2 


= 2 / dr = 2 (r 4 — T 2 ) 


(2.62) 


This is the Euclidean time for the motion from point 2 to point 4. But the full bounce takes twice 
this Euclidean time. Therefore, along a bounce the change in Euclidean time should be the same 
as (3. This means that the bounce solution must be periodic in Euclidean time with a period of /3. 
The same argument goes through for the case of quantum mechanics with n degrees of freedom. 
Similarly the solutions for field theory should have a period of /?. This is what we expect from a 
finite temperature field theory. Keeping in mind that the integral for the bounce action is over a 
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Euclidean time /3, we can interpret P{Et^* — Efv) in Eq. (2.61) as part of the Euclidean action and 
we get the familiar result for the tunneling rate 


"P. ^ [‘^bounce 5'fv] 

r tunn e 


(2.63) 


The shapes of the thermally assisted bounces are different from their quantum counterparts shown 


in Eig. 2.12 They do not have the same symmetry as the 0(4)-symmetric quantum bubbles. 


In addition, the time direction is periodic. However, at very low temperatures the two bounces 


converge to each other. The analog of Eig. 2.12 is shown in Fig. 2.23 



Figure 2.23: The left picture shows a low-temperature thermally assisted tunneling and the right 
a high-temperature thermally assisted tunneling. The Euclidean time r progresses upwards. The 
field configuration at r = 0 and r = /3 is given by <p 2 {x) and a spatial slice through the center gives 
(/> 4 (x). (See Fig. 2.20). In the limit r —)> oo, this pictures reduces to the graph in Fig. 2.12, as 
expected. 


At low temperature the thermal bounces acquire multiple negative modes and they do not 
correspond to tunneling processes. 


2.6 Vacuum decay in curved space time 

Vacuum decay in curved spacetime is a much more formidable problem. Despite all the progress 
made over the last few decades, this phenomenon is not completely understood. On top of the 
usual difficulties associated with the large number of new metric variables and the complexity of 
the Einstein equations, the lack of a well-defined energy in curved spacetime is the main source of 
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the difficulty. The problems of the number of the negative modes, the symmetry of the solution 
with lowest action and also the interpretation of the bounce solutions are not very clear. On the 
other hand, this feature also makes a rich new variety of phenomena possible which are not allowed 
in flat spacetime. For example, not only can the false vacuum decay to a truer (lower) one, it is also 
possible to tunnel from a true vacuum to a false one|21j. The existence of a finite volume horizon 
makes it possible for the whole field to jump to a saddle point of the potential V{4>) (not U[4>{x)]) 
and roll down [22]. This solution is called a Hawking-Moss bounce. 

Initially the field is at its false vacuum at every point. The vacuum energy in this state is 
V{4>{v) and it is constant in spacetime. Therefore the energy-momentum tensor and the curvature 
of the spacetime are constant. Depending on the sign of the corresponding space is called 

a de Sitter, Minkowski or Anti-de Sitter space. These spaces and their basic properties are briefly 
described in Appendix [A| 

The first step towards a calculation of tunneling rates in curved spacetime was taken by Coleman 
and De Luccia (CDL) They used the analogy with the previous flat space results and borrowed 
the formalism developed there. The only change they made was to add the Einstein-Hilbert term 


to the action in Eq. (|2.21|). The action of a scalar field in curved spacetime is then given by 

1 


S = I (fx ^/\g\ 




(2.64) 


In flat space, the tunneling is dominated by solutions with 0(4) symmetry. Although unproven, 
CDL conjectured that the dominating curved spacetime solutions have the same symmetry. This 
means that the metric takes the very simple form 


ds^ = df + , 


(2.65) 


where is a radial parameter and is the metric of a unit three-sphere. The scalar curvature of 
this metric is 

7^ = 4(l-p/-p'') , (2.66) 

where the primes denote differentiation with respect to The scalar field also only depends on ^ 
so we can write it as The variational equations for cf) and p are 
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= l + - v'j , (2.68) 

We still have the freedom to shift ^ by a constant and we choose it so that a zero of p occurs at 
^ = 0. If /9 only has one zero, then ^ ranges from 0 to infinity. In this case the boundary conditions 
are 


</.'( 0 ) = 0 , 

p(0) = 0, 

(/)(cx)) = (?ifv . (2.69) 


The first condition is to avoid a singularity at ^ = 0 in Eq. (2.67), the second one is the way we 
chose the origin of ^ and the third condition is to ensure the finiteness of the Euclidean action. 
These boundary conditions completely determine both fields. The solution obtained this way has 
the topology of If p has a second zero, let’s say at ^max, the third boundary condition is replaced 


by </>^(?max) = 0 to avoid a singularity in Eq. 2.67 This space has the topology of a four-sphere. 

Again, in general it is not possible to solve these equations in a closed form, but fortunately in 
the weak gravity limit, they reproduce their flat space counterparts. 

Let’s look at the main qualitative differences between the flat space and curved space solutions. 

For the case of tunneling from a false to a true vacuum in flat space, the argument using 
undershoot and overshoot solutions ensures the existence of a solution. But in curved spacetime 
there is no guarantee for the existence of these solutions. One famous example is the case of a 
potential E((/>) which gets too flat on top. In this case the CDL bounces disappears [2H - I26] . 

When the CDL bounce has a four-sphere topology, the instanton that describes the transition 
from a true vacuum to a false vacuum has a region close to the true vacuum in its north pole 
(assuming the center of the bubble is at the north pole). It approaches the false vacuum at the 
south pole. We can rotate this instanton and get another instanton that describes the tunneling 


from false to true vacuum|21j. These are shown in Fig. 2.25, Because the instanton action is the 
same, the ratio of the tunneling rates is given by 

r(True —)• False) Q-{Shounce-SiA 


r(False —)■ True) 


('^'bounce Stv) 


= = e“§(''^tv“^fv) 


(2.70) 


Again in the limit where Atv —Af^ <C Atv and Af^, it simplifies to the familiar thermal expression. It 
should not be surprising that in the curved spacetime, the true vacuum can also decay. In contrast 
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to the flat space where a bubble of false vacuum has to collapse because of the surface tension and 
the pressure from the outside, the bubble of false vacuum in curved spacetime can be saved from 
collapsing by the cosmological expansion. Another way to think about it is the existence of the 
thermal bath which provides the energy needed for creation of the bubble of higher energy (false 
vacunm bnbble). But this upward decay is only allowed for a de Sitter to de Sitter transition. 
Bubbles of Anti-de Sitter or Minkowski have an infinite S'fv and therefore cannot up-tnnnel. 

Another difference between flat spacetime and cnrved spacetime tunneling is the existence of a 
new type of solutions which are called Hawking-Moss bounces. Because of the infinite volume of 
the flat space, the field cannot fluctuate at once in every point to reach the top of the barrier. The 
Euclidean action of such a solution is infinite and the associated decay rate is exactly zero. However 
in curved spacetimes with positive energy density, namely the de Sitter spaces, a horizon volume 
is finite and this transition is possible. The field fluctuates to the top of the barrier shown in Fig. 


2.24 and then classically rolls down to the true vacuum. The Euclidean action of this solution is 

^Horizon 


Se = —- 


4G 


= --A^ 


(2.71) 


where Atop is the horizon radius of a de Sitter space with energy Ttop- It is the inverse of the Hubble 
parameter of the space and is given by 


Atop — -f^top ~ 


SvrGHtop 


(2.72) 


V(0) 



4 > 

Figure 2.24: One possible mode of decay is the fluctuation of the field in every point to 4>top- It 
has a finite action and nonzero probability in de Sitter space. 

After subtracting the Euclidean action of the false vacnum, we get the tunneling exponent B 

B=l (A?op - Ai) . 


(2.73) 
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In the limit that Atop — Afv <C Atop and Afv, this expression takes a very simple form 


B 


-[Vtop-Vfv]x 


Horizon Volume 


AU 

= e , 


(2.74) 


where AU is the difference of the energy in a horizon volume of (ptop and Tfv is the temperature 
of the de Sitter space filled with the false vacuum and this suggests that the Hawking-Moss has a 
thermal nature. 

Since the calculations are much easier in the thin-wall limit, we study them here to get a better 
understanding of the CDL solutions. So again we assume that the bounce has three regions. Region 
one extends from the center at ^ = 0 to a value ^ where the wall is located. In this region the field 
is exactly at its true vacuum. There is a narrow transition region at ^ ^ and outside of this region 


the bounce is filled with the false vacuum. Let’s first use Eq. (2.68) in Eq. (2.64) and integrate by 
part. This leads to 


Se = ^7r^ I . (2.75) 

In the thin-wall limit, the action of the bounce inside the wall is the action of a de Sitter or Anti-de 
Sitter space and the action of the wall comes from the surface tension. The general expressions are 
not very simple. But in the simple case where the initial space is a low-energy de Sitter space and 
the transition is to a Minkowski vacuum, the tunneling exponent is 


B = 


Hjiat 



(2.76) 


The tunneling exponent for decay from a Minkowski space to a low-curvature Anti-de Sitter space 
is given by 


B = 


.Rflat 



(2.77) 


These show that the tunneling exponents for the decay from the de Sitter space are smaller than 
the flat space case. Therefore gravity enhances the decays from de Sitter space, as we might expect 
from the new possibility of thermally assisted tunneling. On the other hand it decreases the decay 
rates from Minkowski or Anti-de Sitter spaces. 

Until now all the formalism developed for curved spacetime was by analogy with the flat space 
case. But a more accurate treatment of the solutions is necessary for the interpretation of the 
bounces. Eor large bounces the interpretation is not very clear. In the flat space case shown in 
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Figure 2.25: Left is CDL instanton that describes tunneling from a false vacuum to a true one. Right 
is the same instanton, but the poles interchanged by a rotation. It describes the CDL instanton 
for the decay of the true vacuum. 


Fig. 2.12, at very large negative Euclidean time we get the configuration before the tunneling and 


at r = 0 the configuration afterwards. But in curved spacetime with large bubbles, there is no 
such time which shows the configuration before the tunneling. Brown and Weinberg m clarified 
the interpretation for the case of a fixed background approximation!^ In this approximation the 
variation of the potential in the horizon volume is much smaller than the potential itself. Therefore 
we can assume that the metric is the one of a de Sitter space. This metric in the (Euclidean) static 
patch is 


= P(r)^dT^ + hijdx^dx^ 


= ( ^ ^ j i ^ ~ A2 ) dr^ + r^{d6^ + sin^ 9d4>^) 


-1 


(2.78) 


In the static patch, r starts from 0 and ends at the horizon radius A and the Euclidean time r 
has a period of 27rA. In the fixed background approximation, the constant r surfaces are three- 
dimensional balls bounded by two-spheres at the horizon. In this static background, we can define 
an energy functional with respect to the time-like Killing vector dr 

2 


E= d^xVdeth 


lP(rW^d,^d,4> + P{r)V{4>) 


(2.79) 


where P{r) and hij are defined in Eq. (2.78). Now we can treat the tunneling problem as a flat 


^Here we follow the approach in [2S] 
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space case with a nonzero temperature (the time coordinate has a finite period) and there is no 
ambiguity. The bounce is shown in Fig. 2.26[ 



Figure 2.26: A CDL bounce as it is seen in the static patch. The horizontal dashed lines are the 
constant r hypersurfaces, the shaded region the true vacuum and r ranges from 0 at the center of 
the bubbles to A at the horizon. 

These bounces have the topology of a four-sphere, as expected from the Euclidean version of 
the de Sitter space. Let’s use a coordinate system which shows the four-sphere structure of the 
bounce more clearly 


yi = r sin 6 cos (j) , 

2/2 = rsin^sint/) , 

2/3 = rcosO , 

2/4 = \/h? — r'^ cos (t/A) , 

2/5 = \/— r'^ sin (r/A) . ( 2 . 80 ) 


Here r = —ttA and r = vrA are identified. If we look at the bounce of Fig. 2.26 in the 2 / 4-?/5 plane 


we see the left picture of Fig. 2.27 Rotating it we get the right picture of Fig. 2.27 This one shows 


a r-independent bounce which resembles the purely thermal bounces for the flat space. This means 
that the CDL bounce can be viewed as a thermal bounce. This should not be very surprising, since 
the thermal nature of the de Sitter space is a consequence of quantum mechanics. 
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Figure 2.27: The bounce of Fig. 


2.26 


is shown in the Ui-y^ plane on the left. The lines r = —vrA 


and r = vrA are identified and they show the radial directions. The configuration after tunneling is 


given by half of the hyperplane that goes through the center of the bounce. The other half shows 
the conhguration before the nucleation. On the right the same bounce is shown after a rotation. 
This picture shows a r-independent bounce which resembles a purely thermal tunneling. 
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Chapter 3 

Bounces with 0(3) x 0(2) symmetry 


We studied tunneling in curved spacetime in Sec,2.6 by starting from the argument for the flat 
spacetime case m and generalizing it by analogy to curved spacetime [23]. In this picture, tunnel¬ 
ing is carried by bounce solutions which are the solutions of the Euclidean equations of motion of 
gravity coupled to matter fields. In flat space and zero temperature decay of a single scalar field, 
the solutions with the minimum Euclidean action, and therefore the highest rate, have an 0(4) 
symmetry |6|. There is no such statement for decays in curved spacetime or decays in theories with 
more than one scalar field. However it is generally believed that 0(4) symmetric solutions are also 
dominant in curved spacetime or multi-scalar field theories. In this chapter, we study solutions 
which have the next highest symmetry, 0(3) x 0(2), where 0(3) is the rotational invariance in 
spatial coordinates and 0(2) shows that the solutions are independent of the periodic Euclidean 
time r. 

This chapter is based on |29|. We provide more evidence for the dominance of 0(4) symmetric 
solutions in curved spacetime and also shed more light on the reinterpretation of the Coleman- 
De Luccia bounces in de Sitter space introduced in m where these bounces were understood as 
tunneling in a finite horizon volume and finite temperature. This interpretation was described in 


length in Sec. 2.6 


Static self-gravitating rotationally symmetric bounces have been previously studied in various 
contexts, such as the false vacuum decay in the presence of a black hole (for example [30H3l] and 
the references in these papers ). Our interest in these bounces was first raised by the work of 
Garriga and Megevand |35|. Motivated by the analogy with the thermal production of bubbles 
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in flat spacetime, they studied bounces with 0(3) x 0(2) symmetry and interpreted them as the 
nucleation of a pair of bubbles. Their discussion is in the context of brane nucleation |36l I37j . in 
which the bubble walls are two-branes that separate the true vacuum at the center of a horizon 
volume from the false vacuum that extends to the horizon. This calculation can be applied to the 


case of a field theory in the thin-wall approximation. In Sec, 3.1 we discuss the thin-wall solution 


and review the results of Garriga-Megevand. In Sec 3.2 we set up the formalism that we use for 
studying the decay of a single scalar field in the context of the new symmetry. Our analysis is not 


restricted to the thin-wall limit. The limiting cases are studied in Section 3.3 and we show that when 
the bubble gets large, the thin-wall approximation eventually breaks down. We then summarize 


in Sec. 3.5 and conclude that the 0(4) symmetric solutions are the dominant configuration of 
tunneling in all the cases we studied. 


3.1 Review of the thin-wall results for the case of brane nucleation 

The thin-wall solutions we are considering here are static solutions with rotational symmetrj0 We 
can think about them as a three-ball in spatial coordinates. There is a wall at a radius R and the 
horizon is located at a radius rn > R- This geometry is shown in Fig |3.1[ Inside the wall, the 
metric is that of a de Sitter (dS) space. Because the mass inside the wall is not zero, the metric on 
the outer part is a Schwarzschild-de Sitter (SdS) metric. Therefore the metrics inside and outside 
of the wall (in static patch coordinate^ are 

dsi = C^fi{r)dT^ + f-^{r)dr^, 

dsl = fo{r)dT‘^ + f~^{r)dr‘^+ r^dn‘^ , (3.1) 

where 

fi{r) = 1 - Hfr^ , 

Ur) = l-Hy-‘^, (3.2) 

Here dQ‘^ = dO'^ sin^ 0 dcjP' is the element of area on a unit two-sphere, Hi and Hq the Hubble 
parameters inside and outside of the wall, and M the mass parameter in the SdS metric which is 

^This section and the figures in it are following |35| . 

^These coordinates are described in Appendix A . 








CHAPTER 3. BOUNCES WITH 0{^) x 0(2) SYMMETRY 


42 



Figure 3.1: A schematic view of the thin-wall solution. The wall is located at the radius R and 
inside is the true vacuum. The horizon is located rn and between the wall and the horizon is the 
false vacuum. The bounce solution is constant in the periodic Euclidean time r and apparently 
rotationally invariant. 


related to the total mass of the inside bubble plus the mass of the wall. C is a parameter to make 
the metrics smooth at the tangential directions at the wall. The projections of this geometry in 
the r — T and r — 4> planes are shown in Fig )3.2[ 

The two metrics are glued smoothly at the location of the wall using the Israel junction 
conditions [38]. Knowing a, the surface tension of the wall, and i/tv and H{^, the Hubble pa¬ 
rameters of the true and false vacua is sufficient to determine the mass and radius of the region 
inside the wall. In terms of the extrinsic curvature Kah, the junction conditions can be written as 


[Kab] = -AnCa-iab , 


(3.3) 


where [Kab] is the difference between the extrinsic curvatures on both sides of the wall and 'jab is 
the world-sheet metric. The extrinsic curvature is 


Kab = dr gab 


(3.4) 


We can rewrite Eq. (3.3) in terms of metric components in Eq. (3.1) and obtain 


Qo- gi = -47rGcJ 


9o 9i d ? 


(3.5) 





CHAPTER 3. BOUNCES WITH 0{^) x 0(2) SYMMETRY 


43 



Figure 3.2: The induced geometry of the static instanton in two different planes. Left is the r — r 
plane where it is described by a periodic angular coordinate r and the coordinate along the vertical 
axis is r. Horizon is at the tip of this object rn and the wall is at the location that there is a 
discontinuity in the first derivative due to the energy density of the wall. Right, the same picture, 
but in r — ^ plane. Again (p is the periodic coordinate and r is along the vertical axis. To avoid a 
conical singularity, the period of r must be a very specific number. 


where 

9 o[r) = - , gi(r) = ^- , (3.6) 

r r 

and all the derivatives are with respect to r. These equations completely determine M and R in 
terms of u. Ho and Hi: 


R 


-2 




M = 


AiraR 

3x 


where x is defined by 


e 3a 
la^ IGMl 


+ 


e 3cj \ Hf 
16M |,) 


1/2 


(3.7) 


(3.8) 


In this equation, e is the difference between the true and false vacuum energies. These equations 
imply that the tension cannot get too large, which is what we would naively expect. To compensate 
for the increase of the tension, we have to increase the radius of the bubble. But because of the 
existence of a horizon, there is not enough room for the radius to get too large and therefore there 
should be a maximum allowed tension ajsr which is shown to be 


aN = 2MIJ3H2 - Hf . 


(3.9) 
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When a approaches the outside metric starts developing a double root at the horizon and 


beyond un, there is no solution. We will see in Sec 3.2.5 that this picture gets modified in the 


context of a scalar field theory and that before reaching this limit the thin-wall approximation 
breaks down. 

The periodicity of r is a measure of the temperature felt by the system. It can be determined 
by demanding that there be no conical singularity at the horizon. Near the horizon, we can expand 
the metric in powers of distance from the horizon 


/o(r) = B'^irn-r) , 


(3.10) 


where B^ is positive to make the metric positive inside the horizon. It is straightforward to calculate 
it: 


B^ = - — . (3.11) 

th 


Using the standard change of variables 


2 _ b"^ 

P= -r , 7 = —t , 


(3.12) 


we can rewrite the metric near horizon as 


+ dp^ + rjjdCL^ . 


(3.13) 


To avoid a conical singularity, 7 must have a period of 27r and therefore the period of the time 
coordinate is given by 


/3 = 


AirrH 


2-Kr‘jj 


2>Hlr‘jj - 1 th- SGM ■ 
In this approximation, the action of the bounce solution is 


c _ 


(3.14) 


(3.15) 


We will show in Sec. |3.2| that this result is in fact exact. Numerical calculations in [35] show that 
the action for the 0(3) x 0(2) in the thin-wall approximation is higher than their 0(4) symmetric 
counterparts in thin-wall limit. 
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3.2 Inclusion of a scalar field 


In this section, we set up the formalism for the rest of this chapter. Our goal is to study the decay 
of a scalar field cj) with a potential V{4>) in de Sitter space. The initial Hubble parameter is set by 
Tfv (or Vtv for decays from a true to a false vacuum) and the bounce solution has an 0(3) x 0(2) 


symmetry. Later in Sec. 3.4 we will scan over possible H’s. The scalar field and the metric are 
independent of the Euclidean time r and only depend on the spherical coordinate radius r. The 
most general r-independent metric compatible with spherical symmetry can be written as 


= B{r)‘^dT^ + A{r)dr‘^ + r^dO^ 


(3.16) 


We still have the freedom to rescale r. This will modify B{r) by a rescaling. The normalization 
for B is explained in the next paragraph. Each spatial slice is a three-ball and the horizon radius 
rH is located at the zero of 5 ^'’ = 1/A. 

In order to get a non-singular solution at the origin, we need ^(0) = 1 to avoid a conical singu¬ 
larity. We choose the normalization of B by demanding that A{rH)B[rH) = 1 • This configuration 
corresponds to thermal nucleation of a bubble at the center of the horizon volume. The bounce 


solution is drawn in Eig. 3.3 on a four-sphere. The analysis of m predicts the production of a 
single bubble, while the conventional interpretation is the production of a pair of bubbles. The 
former seems to be a more natural interpretation. 


3.2.1 Periodicity of Euclidean time 


In thermal tunneling in flat space, the temperature determines the period of the Euclidean time 
and can be chosen arbitrarily. In curved spacetime, however, we expect the temperature and hence 
/?, the period of the Euclidean time to be determined by the metric]^ Therefore we assume that j3 
is determined by the surface gravity at the horizon of the Lorentzian counterpart of our Euclidean 
metric: 


/3 


— = dvr 
T 




/I - 1/2 


J rn 


(3.17) 


only enters in our calculations as an overall multiplicative factor in the action and does not otherwise affect the 
solution. It would be a straightforward matter to extend our results to an arbitrary temperature, but the physical 
origin of the temperature would be less compelling. 
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Figure 3.3: Two views of an 0(3) x 0(2)-symmetric bounce. These are from the same viewpoints 


as those of the CDL bounce in 2.27, (The views on the right would be the same if viewed from 
along any axis in the X 4 -X 5 plane.) On the left-hand view the radial dashed lines correspond to the 
similar lines in Fig. 2.27[ The bounce is independent of the Euclidean time. The three- dimensional 
configuration corresponding to any of these lines contains a single true vacuum critical bubble in 
the center of a false vacuum horizon volume. 


where primes denote differentiation with respect to r and the subscript indicates that all the 
quantities are to be evaluated at the horizon. Our normalization condition that A{rH)B[rH) = 1, 


the fact that at the horizon ^ ^ = 0, and Eq. (3.17) lead to 

d 


B' = 


dr 


(AB) 


= {AB)' 


+ {AB) 


(3.18) 


Therefore (3 simplifies to 

P = -ATi[B'{rH)]~\ (3.19) 

with the minus sign coming from the fact that B'[th) is negative. If the space were pure de Sitter, 
(3 would be 2'KrH- It should remain the same order of magnitude for the bounce solution. 


There is another way to deduce Eq. (3.17). At the horizon 1/A = 0 and AB = 1. This means 
that B{rH) = 0. Therefore both (^) and B vanish at the horizon. Assuming that they do not 
develop double roots, to the first order of Taylor expansion 


B{r) = -B'{rH-r), 

= C‘^{rH - r) , 


A{r) 


(3.20) 


where = — (^)^ In terms of the new variable y = metric near the horizon is 

B'C^ 


ds^ = dy^ -|- 




4 


(3.21) 
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Avoiding a conical singularity in the y-r plane leads to 


/? = dvr {B'C^) = dvr 




h - 1/2 


which is consistent with Eq. (3.17) . 


(3.22) 


3.2.2 Equations of motion and boundary conditions 

The Euclidean action for a scalar field coupled to gravity i^ 


S = j (^- 


1 


IhvrG 


n + c 


matter 


(3.23) 


where IZ is the scalar curvature and ^matter is 


>Cmatter = -d^cjyd^^cj) + V{4>) . 


(3.24) 


The nontrivial components of the Einstein tensor and the energy-momentum tensor for the metric 


in Eq. (3.16) and the action in Eq. (3.23) are 

B 


{-A + A^ + rA') , 


Grr = 


Trr = -Bf^^ + V{cP)] , 

Trr = A{^-V{ct>)] , 


j,2^2 

B-AB + rB' 


= 


AA^B'^ + rBA'{AB + rB') + A (-4B^ + r^B'^ - 2rB{2B' + rB")) 


(3.25) 


2r2^2^2 

The gravitational and matter contributions to the action inside a horizon volume are (k = SttG) 

rrn rBA'{AB + rB') - AAB^ + r^AB'^ - 2rAB{2B' + rB")] 

Tidr- 


S, 


gravity 


5, 


JO 

r^H 


k{AB)1 


matter — 


_ / M2 

drA'irr'^y/AB i — + V{4>) 


(3.26) 


^Please notice that the ^ is a positive quantity in Euclidean spacetimes. 
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The equations of motion arising from varying the action in Eq. (3.26) with respect to (p, A and B 
ar^ 


» - i 




+ kC 




/^2 


2 A 


+ V(P) 


0 = 


^ d / r'^^JABp'\ 


dr 


A I 


-r‘^VAB'^ 


We can integrate the second of these equations and get 


nu 


A[r)B{r) = A{rH)B{rH) { —K I dr r' 


We chose A{r}{)B{r}{) = 1 in 3.2.1 Therefore 


r^H 


A{r)B{r) = exp ( —k I dr' r' {p')^ 


(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 


We can eliminate B from these equations so the only remaining fields to determine are A{r) and 
4>{r) . Before proceeding further, it is convenient to define new variables /(r) and AA{r) by 


■'' ' yl{r) r 


We can rewrite Eq. (3.27) in terms of / as 

0 = ^ [^ (/ - 1 )] + 

This can also be written in the familiar form 




M' = dvrr^ 


{p'Yf + v{p) 


Substituting Eq. (3.28) in Eq. (3.29) gives 


0 = / 


p" 4 — p' + ATvGr^p'Y 

r 




(3.32) 


(3.33) 


(3.34) 


(3.35) 


This equation together with Eq. (3.33) gives two equations involving only / and p. They are second 
order in p and first order in / and need three boundary conditions. To avoid a conical singularity 
at the origin and guarantee the smoothness of p, we would need 

/(O) = 1, 

p'{0) = 0. 


“These equations can also be obtained by plugging the ansatz in Eq. (3.161 into the full field equations. 


(3.36) 
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The first boundary condition prevents a conical singularity at the origin (since dO,'^ is the standard 


two-sphere metric) and the second one is a consequence of Eq. (3.35). To completely determine the 
fields, we need one more boundary condition. At the horizon, / = 0 and if we assume that 4>" and 
4>' are finite, it imposes 


/V- 


dV 


= 0 


(3.37) 


r=rH 


In fact this condition holds even for singular cp' if cp' grows slower than (rn — as we approach 

the horizon. Such a divergence does not imply a divergent action density and is merely a coordinate 
singularity. In fact, choosing the proper distance from the horizon as a new coordinate makes <p' 
finite. 

There is another subtlety that we have to address here. With the periodicity introduced in 


Eq. (3.17), the four-sphere obtained by identifying all of the boundary two-spheres is a smooth 


manifold, just as in the CDL case. However, as in CDL, we do not require that (p to be smooth 


on this two-sphere. Smoothness of (p on the four-sphere requires (p'{r}{) = 0. Equation (3.35) 
implies that dV/dcp = 0 at the horizon which in turn means that the field must be precisely at a 
vacuum value (or else exactly at the top of the barrier, as in the Hawking-Moss solution). This 
is only possible in the (unattainable) limit in which the thin-wall approximation is exact. Having 
(p'{fH) 7^ 0 causes a discontinuity in slope along a path passing through the two-sphere. However, 
this discontinuity does not cause any cost in action and therefore it is fine for cp' to be non-smooth. 


Knowing (p and / and using Eq. (3.31) determines B. 

3.2.3 Bounce action 

One of the interesting features of these 0(3) x 0(2)-symmetric bounces is the remarkable coincidence 
that the on-shell Lagrangian density is a total derivative and therefore the action only depends on 


the boundary values of the fields. Starting from Eq. (3.26) and integrating by parts, we can obtain 

Aicp'f 


*^“40 


TD! 1 r=rH pth 

'^= +i7r/3 : 

yAB\r=Q Jo 


r d 



9 

\ Kdr 

r{ 




2 A 


+ V(cP) 


. (3.38) 


Upon using Eq. ( |3.27 ), the terms in the curly brackets vanish. Now using Eq. (3.17), we can 
simplify the first term 


Q — _ZL 2 

'^bounce — ^ 


(3.39) 
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The action only depends on the horizon area. The tunneling exponent B is the difference between 
the bounce action and the false vacuum action. 


i3 = 5b, 


(Ounce 


5fv = ^(A?.-4) , 


(3.40) 


where 



(3.41) 


It agrees with the thin-wall results of [35] and is expected on more general grounds. As shown in 
|39l HQ], the action of a static Euclidean solution with periodicity (I is Se = 13E — S, where E is 


the total energy and S is the entropy. For a solution without a boundary the total energy vanishes 
and the action depends on the entropy, which in turn depends on the horizon area. 

3.2.4 Comments on the thin-wall approximation 

The thin-wall approximation is a useful tool for gaining an intuitive understanding of tunneling 


processes. In this section we describe the conditions for validity of this approximation in 0(3) x 


0(2)-symmetric solutions. In this approximation, the bounce is approximated by a region of pure 


true vacuum inside the wall, a transition region which is the wall and the outside region which is 


in the pure false vacuum. To use this approximation, we require two essential conditions. First, 


the thickness of the wall (transition region) must be much smaller than the radius of the wall so 


we can approximate the wall as being locally planar. The second condition concerns the position 
dependence of the surface tension a of the wall. Although in most of the previous tunneling 
calculations, a was considered to be position-independent, it can depend on the position. To make 
the thin-wall approximation valid, the fractional variation of a must be small through the wall. 

In flat spacetime, the thickness of the wall is determined by the shape of the potential barrier 
separating true and false vacua. The wall radius is determined by the ratio of cr to e, the difference 
between the false and true vacuum energy densities. Reducing e and meanwhile keeping the shape of 
the potential almost unchanged increases the wall radius and therefore it is always possible to obtain 
a thin-wall by making e small enough. This is the case for both 0(4)-symmetric zero-temperature 
and for high-temperature solutions that describe a critical bubble in three dimensions. 


The situation is different in the presence of gravity because it introduces a new length scale, 
the horizon radius. Therefore the radius of the wall gets an upper limit equal to the horizon length 
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and the true vacuum region cannot get arbitrarily large. To even have a possibility of a thin-wall 
approximation, the shape of the potential must be chosen in such a way that the natural width 
of the wall is small compared to the horizon radins. If this criteria is satisfied, then by making e 


small, we can achieve a thin-wall solution. In Sec,3.2.5 we argue that this is not always the case 
for our 0(3) x 0(2) solutions. 


3.2.5 No Nariai limit in the thin-wall approximation 


In the context of brane nucleation it was claimed in [35] that as the surface tension a approaches 
the critical snrface tension cjjv, the wall approaches the horizon and can get arbitrarily close to 
it. Meanwhile, the metric develops a donble-root at the horizon. This is the Nariai limit [411142j . 
Here we show that it is not possible in a field theoretical set np to achieve a thin-wall solution 
which approaches the Nariai solution. If there is such a solution, near the horizon the metric can 
be approximated as f{r) = B{Rh — r)^, where H is a positive constant. Because of the thin-wall 
approximation, the field near the horizon should be very close to the false vacuum. Expanding the 
potential near the false vacunm {P and Q are determined by V), 

= Vi, + ip(0 - + ... . (3.42) 


It is more convenient to write Eq. (3.33) and Eq. (3.35) in terms of y = Rh — r and ^ = cj) — (/){, 


In the remainder of this subsection the primes denote differentiation with respect to y we obtain 

s /2 


1 + Rnf -yf'-f = SnCiRn - yf 


$"/ + $ 7 ' - 


-^c^7 - ^7tG{Rh - yW^f = ^ . 
Rn-y 


(3.43) 

(3.44) 


The boundary conditions are 




dV 


y=0 


<I>'{Rh) = 0 , 


/{Rh) = I . 


(3.45) 

(3.46) 

(3.47) 
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Notice that y = Rh is the same as r = 0. At the horizon, the equations of motion simplify to 


/2 


1 — SttGR^ (/ — h Vh ) , 


a'2 


BRjj = -1 + AirGRij [f^ + V 


(3.48) 

(3.49) 

(3.50) 


Here we get a second order and a first order differential equation and two unknowns, B and Rh, and 
five boundary conditions which make it possible to solve the equations near the horizon. Expanding 
‘h in terms of y to the first nonzero order, we get <h(y) = A^/y where 


A = 


ttGR 


P 

B 


Plugging back this solution for into the metric equation gives 


BRir = -I + AttGR 


H 


/^2 


+ U] = -1 + AttGR^A'^ 


B P 


This leads to 


B 


AP\ 


1 


5r 1-^ +5^=0, 


B J 


R 2 

Hr 


(3.51) 


(3.52) 


(3.53) 


where H is a positive number. This shows that the metric cannot develop a double root at the 
horizon. Because the left side is always a nonzero positive number, adding a very small linear term 
to the metric does not make the situation better. Therefore not only can the metric not develop a 
double root in the thin-wall limit, but also it cannot get arbitrarily close to having a double root. 


3.3 Limiting cases 


We could not find closed form solutions of the field equations Eq. (3.27)-Eq. (3.29) and had to 
use numerical methods to understand the behavior of bounce. However there are limiting cases in 
which we could get more analytic insight. We present them in this section. 


3.3.1 Weak gravity limit 

If the gravitational constant G is small, we should expect to recover the flat space results with small 
corrections due to the nonzero temperature of the horizon. If G is sufficiently small, cj) varies from a 
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value near the true vacuum to one which is exponentially close to in an interval 0 < r < f, where 
f is much smaller than the horizon radius Af^. In region r < r < Af^, the held almost assumes the 
false vacuum value and the space is well approximated by a SdS metric. In region 0 < r < f, the 


held does not feel the curvature of the spacetime and Eq. (3.34) should be well approximated by 


the hat spacetime equations. Outside this region, we can integrate Eq. (3.34) to get 


M Eflat + 2^ ^2^, 


(3.54) 


where E\^ is the energy of the critical bubble in hat spacetime. The horizon appears at the place 
where / vanishes or, equivalently, 2GAA.{r) = r. To hrst nonzero order in G, we get 


ru = Afv - GEflat . 


(3.55) 


Since the bounce action only depends on the horizon radius, the tunneling exponent Eq. (3.40) is 
given by 


B = 27rAfvEflat = 


Eflat 


(3.56) 


where from, Appendijij^ Tfv = l/27rAfv. This is the expected Boltzmann exponent that we would 
obtain for nucleation of a critical bubble in hat spacetime at a temperature T = T^s. However, 
because T <C l/r, the high temperature bounce has multiple negative modes and must be discarded. 


3.3.2 Strong gravity and oscillating bounces 

Increasing the gravitational constant shrinks the horizon radius Afv and leaves less and less space 
for the wall separating the two vacua. As Afv becomes comparable to or smaller than the natural 
width of the wall, cp is restricted to an increasingly narrow range of values around the top of the 
barrier. Eventually i?!)(0) and <p{rH) merge and the held becomes spatially homogeneous. This is 
the usual Hawking-Moss bounce whose 0(5) symmetry contains both 0(4) and 0(3) x 0(2) as 
subgroups. This is illustrated in Eig |3.4[ This is very similar to the behavior of 0(4)-symmetric 
bounces described in [231126] . 

When the solution approaches the Hawking-Moss bounce, 6(p = (l){r) — (ptop becomes small and 


this enables us to use an expansion in powers of 6(p in Eq. (3.27)-Eq. (3.29). Because the variation 
of the held is small, the background metric remains very close to a de Sitter metric with A = A^^^p. 
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V(0) 



Figure 3.4: Increasing G leaves less and less room for the wall. As a result cj) varies less over the 
horizon volume. In this graph the markers on the left of the barrier show (/)(0) and the markers 
on the right of the barrier 4>{rH) for different G. Black circles, red diamonds, green triangles and 
the orange square show (/>(0) and (pirn) in order of increasing G. For small G (black circles), the 
held starts very close to the true vacuum at the center and ends very close to the false vacunm 
at the horizon. Increasing G makes the change in cj) smaller. Eventnally, for a critical value of G 
which is calculated in this section, the two points merge(orange square). This corresponds to a 
Hawking-Moss solntion. 

The analysis presented here is along the lines of |25] . To zeroth order in 54>, the solution to held 
equations are 


</’(^) — (ptop ) 

/(r) = 1 - , where H'^ = ^^Utop ■ 

O 

Expanding cp and / as 

(p{r) = (ptop + S(p{r) , 

fir) = Mr) + 6fir) = l-H\^ + 5f{T) , 
VicP) = Vtop + H^I^-^cP^ + ^cP^ + ^cP^^ . 


The boundary conditions for 6f and 5(p are 


dV 

dcp 


A 


top 


(3.57) 


(3.58) 


6cP'{0) = 0 , ,5/(0) = h/(Atop) = 0 , h,/'(Atop) [-2Htop + <5/'(Atop)] 


(3.59) 
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Using Eq. (3.58) in Eq. (3.33) and Eq. (3.35) leads to 


d 


r^dr 


[r^dr^'ifo + df)] + A7TGr{6(j)'f{fo + df) = H^-^dc^ + b{dcj))^ + X{dcpy 


^ [r (1 - /o - df)] = Kr^ 


^{df)')'^f + Vtop + h'^ (-^{d4>)'^ + -^{d4>)^ + -^{d4>y 


(3.60) 

(3.61) 


Equation Eq. (3.61) implies that the first correction in df is of order dcj)'^. This simplifies the first 


equation immensely: 

d(j)”(l - H^r^) + d(j)'(-AH^r + -) + 47TGr(d(j)')^(l - H^r'^) = H^(-^d(p + b(d(f>)^ + \(d(j))^) . (3.62) 

r 

Defining y = rHtop, which brings the horizon to y = 1, and keeping only the first nonzero order 
in df) (the primes for the rest of this section denote differentiation with respect to y) leads to 


dcp'fl — y^) — Aydcj)' H— df)' + fldf) = 0 . 

y 

To bring this equation into a more familiar form we change variables to 


(3.63) 


g{y) = ydf>{y) . 


(3.64) 


This leads to 


(1 - y^)g” - 2yg' + (/3 + 2)5 = 0 , 


(3.65) 


which is the Legendre equation of order — i + 4^4/3 + 9. The boundary conditions Eq. (3.59) i 


m 


terms of y are 


df)'{0) = 0 


'f'l = - 
■' 'y=^ df) 


(3.66) 


y=l 


After some simple manipulation using Legendre’s equation, the boundary condition at y = 1 
becomes 

(1 - y^)g" = 0 , (3.67) 

As long as g" is not more divergent than 1/(1 — y), this is satisfied automatically. The most 


general solution to Eq. (3.65) is a linear combination of Pn{y) and Qn{y)- Since all the Qn’s are 
logarithmically divergent at y = 1, they can’t satisfy this equation. Therefore the solution must be 
proportional to Pn{y)- Let’s look at the expansion of g around zero: 


g{y) = ao + aiy + a 2 y^ + C>(y^) . 


(3.68) 
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The boundary condition at y = 0 in terms of g and y is 


= 0 . (3.69) 

y=0 

To satisfy this boundary condition we need oq = 02 = 0 . The Taylor expansion of Pn{Q) is 



Pn{y) = 


TT 


+ 


n(n + l)y' 7 r n{n + 2){n? — 1)^ 


r(i^)r(2^) ' 2r(^)r(3f^)^+ 8r(3^)r(^) 


+ 0{y'^) 


(3.70) 


2 y 2 

The only possible way for uq and 02 to vanish is to choose n to be an odd integer. This makes 


the Gamma function infinite. This restricts j3 defined in Eq. (3.58) to (3 = n{n + 1) — 2 where n 
is an odd integeij^ Therefore, if we drop all the non-linear terms in Eq. (3.62), the solutions exist 
only for very special values of (3. If we add the nonlinearities, it relaxes this condition. Since the 
Legendre polynomials form an orthogonal basis, we can expand any solution of the full nonlinear 


equation in this basis. We can rewrite Eq. (3.62) in terms of g 


(1 - y‘^)9" - W + (/? + 2)g - bg^ - \g^ = 0 


(3.71) 


Expanding the solution in terms of Pm{Q) (M is an odd integer): 


Vl'^l M 


(3.72) 


and plugging back into Eq. (3.71) leads to 


PM{y) (/? + 2 - M(M + 1)) Am 


M 


Wfrr ~ ^2 qijkmAjAjAk ) = 0 > 

V l-^l I,J I,J,K j 

(3.73) 


where the p’s and g’s are defined by 


PiJM = f2 dyPi{y)Pj{y)PM{y) , 

qijKM = ^^^2^^^ flidyPi{y)Pj{y)PKiy)PMiy) ■ 


(3.74) 

(3.75) 


Since the nonlinear terms are small, we expect the solution to the field equations to be close to 
the solutions for the linearized equation. Let’s assume that the solution is very close to one of 

®This is very similar to the results obtained for 0(4)-symmetric solutions in |25| . There /f had to be equal to 
-I- 3) for integer N's. 
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the PN{y) modes where A = /? + 2 — N{N + 1) is small (N is an odd integer). Therefore all the 


expansion coefficients Am{M ^ N) in Eq. (3.72) are very small compared to A^. Solving for Am 
to the lowest nonzero order leads to 


^m(/3 + 2 - M{M + 1)) - -^J^bAlfPNNM = 0 , M ^ N , 

An^ - {Aj^pNNN + “^AnAmPnmn) - sgn{X)A^qNP^NP^ = 0 . 

From the first equation, we can solve for Am- Plugging back into the second equation we have: 


(3.76) 

(3.77) 


cAff + ti^pnnnA j\f — A + 0{A^) — 0 , 


where 


2^2 


c = sgn{X)qNNNN + -jyr ^ 


Pnnm Pnmn 


|A| ^/3 + 2-M(M + l) 


(3.78) 


(3.79) 


Since N is an odd integer, P/v(?/) is an odd function and therefore pnnn vanishes. Solving for An 
we have: 

[a 

An = \—. (3.80) 

V c 

Let’s assume that 6 = 0. The N = 1 solution corresponds to A = /3, which is positive.This makes 
Ai a real and acceptable coefficient. In this case g{y) = Pi{y) = y which corresponds to a field 
configuration which is constant because (f>{y) = This solution is not acceptable because it is 


inconsistent with the change of variable in Eq.(3.64). So the first useful solution corresponds to 


the case N = 3, which which makes A = /3 — 10. Because A is under the square root, it cannot be 
negative and this sets a minimum value for (3, 


V”(0) 


(3.81) 


These are the solutions that start from one side of the barrier and end on the other side. They 
are acceptable tunneling conhgurations. ^3 gets smaller when (3 approaches 10 from above and 
the bounce solution gets closer and closer to the top of the barrier and merges with Hawking-Moss 
solution at /3 = 10. For larger values of N we will have similar situations. For example the N = 5 
solution exists only for /3 > 28. This solution crosses the barrier twice. It has more than one 
negative mode and therefore not an acceptable tunneling configuration. In general, solutions for 
N = 3 + 2k cross the barrier k times and merge with Hawking-Moss when f3 = 2(A:-|-l)(2A:-|-5).None 
of them are acceptable as tunneling configurations. Some of these solutions are shown in Fig. |3.5[ 
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Figure 3.5: Oscillating bounces for different values of N. Green solid line, blue dashed, red dotted 
and black dot-dashed lines represent solutions for N = 1, 3, 5 and 7. Only the blue line represents 
a tunneling instanton. 


3.3.3 Large vacuum energy 


Another limiting case that we can study is to keep G fixed and raise the vacuum energies by adding 
a large constant to V{^). On a pure vacuum solution, increasing the vacuum energy has the same 
effect on the metric as increasing G and shrinks the horizon radius. This leaves less room for the 
wall and therefore the field spends more time near the top of the barrier and the analysis of the 
previous section applies. If the vacuum energy gets too large, the 0(3) x 0(2)-symmetric solution 


merges with Hawking-Moss at the critical value of /3 found in Eq. (3.81). 


3.4 Numerical results 

Because it was not possible to solve the field equations in a closed form, we tried to understand the 
qualitative behavior of these bounces by choosing a quartic potential. If we look at the expansion of 
F((/)) near iji = 0, we need to at least keep terms to quartic order to capture the tunneling features. 
We considered quartic potentials with two minima. By shifting (p, we set up the two vacua to be 
equally separated around the field space origin, with cptv = v and = —v. Any such potential 
can be written in the form 

(^Gov‘^ - kv^(p - + ^v(p^ + 


V{(P) = A 


(3.82) 
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with 0 < A: < 1. The top of the barrier separating the two vacua is located at (j) = —kv. The theory 
is therefore characterized by four dimensionless quantities: Cq, k and 


h = . 


(3.83) 


The dependence on A is very simple and we can scale it out. If 4>{x) = g{x) is a bounce solution 
for a given value A, then ())(x) = gi'yx) will be a solution with A replaced by y^A. Therefore the 
strong coupling limit (large A) is related to the weak coupling limit by a rescaling of distances. 
The action S of the bounces replaced by 'y~^XS . Therefore the classical action has a nontrivial 
dependence on only the three remaining dimensionless parameters.We found it more useful to work 
with a different set of alternative parameters, 


e = ^ (Vfv - Ttv) = ^Xk , 

3 

« = ^(k^top-I"fv) = ^(3-A:)(l + /^)^ 

1 / 1 2A' A 

Uo = ^^tv = A(^---- + CoJ . (3.84) 


These parameters characterize the splitting between the energies of the true and false vacua, and 
between the top of the barrier and the false vacuum and the energy of the true vacuum, all in units 
of v^. In the subsequent sections we study the dependence of the bounce solution on each of these 


parameters. We change one while keeping the others fixed and solve the bounce equations (3.27)- 


(3.29) numerically. Such a variation will correspond to a more complicated path in the X — Cq — k 
plane and may extend to regions where A is large and therefore the coupling is strong. However, 
we are not concerned with this strong coupling since any large coupling solution can be mapped 
onto a weak coupling solution by the rescaling introduced in the beginning of this section. 


3.4.1 Varying G 

As mentioned in Sec. |3.5[ increasing the gravitational constant causes the solution to approach the 
Hawking-Moss. In Fig j3.6t the numerical solutions for different values of G are shown. In these 
graphs I/q, e and a are held fixed. Increasing G causes the field to get close to the Hawking-Moss 


^This is explained in details in chapter 12 of |28j . 
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(purple dashed line) and eventually the 0(3) x 0(2) merges with the Hawking-Moss at a value of 


/3 = 9.998 which is in a very good agreement with what we found in Eq. (3.81). In this limit cj) 


becomes spatially homogenous and / corresponds to a pure de Sitter space with A = Atop. 


3.4.2 Varying Uq 


We also argued in Sec 3.3.3 that the effect of increasing V{(l)) by a constant density is qualitatively 


very similar to that of increasing G. Our numerical solutions show that this statement is correct and 


the solution merges with Hawking-Moss at a value given by Eq. (3.81). These numerical solutions 
are shown in Eig|3.7[ 


3.4.3 Varying e 

In flat space, changing e, the dimensionless difference between the false and true vacuum energy 
densities, changes the radius of bubble. Eor smaller e the bubble must get bigger to compensate for 


the surface tension of the bubble wall. In Eig. 3.8 we showed the field profile and metric for different 
values of e. We also included graphs for e = 0, which corresponds to two degenerate minima, and 
also graphs for negative e’s which show tunneling from a true vacuum to a false vacuum. Neither of 
these cases are possible in flat spacetime at zero temperature, where the bubble radius approaches 
infinity as e goes to zero. However, the situation is different in the presence of gravity. If the true 
and false vacua are both de Sitter, it is possible to tunnel upward [21]. What saves the bubble 
from collapsing due to both surface tension and vacuum energy differences in up tunneling is the 
Hubble flow of the de Sitter space provided that the initial bubble size is large enough. 

There are qualitative differences between the CDL (0(4) symmetric) and 0(3) x 0(2)-symmetric 
cases. In the CDL case, tunneling from true to false vacuum proceeds by nucleation of bubble of 
true vacuum (which we can assume to be centered about the “north pole” of the four-sphere by 
using the de Sitter symmetry groupie and the false vacuum region around the south pole. To 
describe the tunneling from false to true vacuum, we only need to to interchange the south and 
north pole labels. The situation is different for 0(3) x 0(2)-symmetric bounces. The spacelike 
slices are three-balls with the bubble of the new phase at the center of a horizon volume and the 


®The geometry of de Sitter space and the terms north and south poles are explained in Appendix. 
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Figure 3.6: Evolution of the bounce as the gravitational constant is varied. The first two panels 
show (f> as a function of (a) r and (b) r/rn- In both cases the black dot indicates the value at 
the top of the barrier. The third panel shows / = 1/^ as a function of r. Reading from left to 
right, the short-dashed purple, long-dashed brown, dotted red, dot-dashed blue, and solid black 
lines correspond to SttGv'^ equal to 0.711, 0.704, 0.628, 0.251, and 0.126. In all cases Uq = 2, a = 3, 
and e = 0.3. 
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(a) (b) 



V r 


(c) 


Figure 3.7: Evolution of the bounce as the true vacuum energy Uq is varied. The first two panels 
show (j) as a function of (a) r and (b) r/rn- In both cases the black dot indicates the value at 
the top of the barrier. The third panel shows / = 1/^ as a function of r. Reading from left to 
right, the short-dashed purple, long-dashed brown, dotted red, dot-dashed blue, and solid black 
lines correspond to Uq equal to 2, 5, 10, 2, and 26.8 In all cases SttGv'^ = 0.126, a = 3, and e = 0.3. 
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old vacuum at the outer edge, extending to the horizon. There is no longer a symmetry between 


the two regions and therefore the bounce solutions are different, as can be seen in Fig. 3.8 


3.4.4 Varying a 

The most interesting behavior was that which corresponds to changing the barrier height a. In¬ 
creasing the barrier height has two effects. First, it increases the surface tension of the bubble 
which in turn (as in flat spacetime) increases the bubble radius. Second,it makes the wall thinner. 


The situation for bubble nucleation in de Sitter spacetime is shown in Figs. 3.9 and 3.10 For small 
values of a, illustrated in Fig j3.9t the field starts very close to the true vacuum at the center and 
approaches the false vacuum close to the horizon. However, in contrast with the flat spacetime 
case, the wall gets closer to the center of the bubble as a and therefore a increases. This is due 
to the shrinkage of the horizon radius which makes all of the evolution happen on a shorter scale. 
However, when the field profiles are viewed as functions of r/r/f, the results resemble the flat space 
case. In both cases, whether as a function of r or r/rn, the walls get thinner when a gets larger 
and they get closer to the horizon. 

We may expect from this pattern that by increasing a further we might make the wall arbitrarily 
narrow and close to the horizon. However, this trend does not continue when the wall starts reaching 


the horizon for larger values of a. For very large values of a, shown in Fig 3.10 the wall is cut 
by horizon and there is no false vacuum region near the horizon. This means that the thin-wall 


approximation has already been broken, as proved in Sec 3.2.5 ^{ru) falls visibly short of the false 
vacuum value as a gets bigger. [For the largest a shown here, = —0.85u.] The wall does not 

get thinner and, as shown in panel b of Fig j3.10[ the shape of the wall does not change much. The 
field prohles as functions of r jrn are almost indistinguishable. In this region the behavior of the 
field as a function of r near the horizon is 


= (pirn) + , 


(3.85) 


for a constant B that depends on the potential and therefore (j)' diverges. But this divergence is a 
matter of a bad coordinate system. Choosing the physical distance defined by 


J’phys = / dr'sjA{r') , 

Jo 


(3.86) 
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V r rpH 

(a) (b) 



V r 


(c) 

Figure 3.8: Behavior of the bounce as e is varied. Again, panels (a) and (b) show cp as a function 
of r and of rjrn, with the black dot indicating the value at the top of the barrier. Panel (c) shows 
the metric function / = 1/A. Reading from left to right, the dashed brown, dotted red, dot-dashed 
blue, and solid black lines correspond to e equal to 2, 1, 0, and —1. For all cases Uq = 2, a = 3, 
and SttGv'^ = 0.126. 
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(a) (b) 



V r 


(c) 


Figure 3.9: Variation of the bounce with barrier height for small a with the short-dashed purple, 
long-dashed brown, dotted red, dot-dashed blue, and solid black lines corresponding to a equal to 
1, 3, 5, 12, and 15. For all of these Uq = 2, e = 0.3, and = 0.126. As in the previous figures, 

the black dot indicates the value of the field at the top of the barrier. 
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Figure 3.10; Variation of the bounce with barrier height for large a, with the short-dashed purple, 
long-dashed brown, dotted red, dot-dashed blue, and solid black lines corresponding to a equal to 
15, 20, 25, 40, and 100. For all of these Uq = 2, e = 0.3, and SvrGu^ = 0.126. As in the previous 
hgures, the black dot indicates the value of the held at the top of the barrier. 
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as a coordinate removes this divergence. Even when (p'(r) develops a singularity, the boundary 


conditions introduced in Eq. (3.36) and Eq. (3.37) still hold. The actions of these bounces, which 
only depend on the horizon radius, remain finite and therefore we can neglect this coordinate 
singularity. 


3.5 Summary and conclusions 


In this chapter we studied bounces with 0(3) x 0(2) symmetry. These bounce solutions show 
the thermal nucleation of a single bubble at the center of a horizon volume. The conventional 
interpretation of these bounces is the creation of two bubbles. However following the interpretation 
of |27j and taking a constant r slice through the center of the bubble (which is topologically a 
three-sphere) gives the configuration of metric and the scalar field inside a horizon volume without 
any reference to quantities beyond the horizon. 

We used general quartic potentials and scanned the parameter space of these potentials. We 
showed numerically and analytically that increasing the gravitational constant G or adding a large 
positive constant to the field potential V{4>) drives the bounce towards the Hawking-Moss solution 
and calculated the threshold for G at which the independent 0(3) x 0(2) solutions ceases to exist 


in Eq. (3.81). 


The effects of decreasing e, the energy difference between the true and false vacua, is very 
different from their flat space counterparts. In flat space decreasing e makes the bubble radius 
large. As e tends to zero, this radius goes to infinity and there is no tunneling from false to true 
vacuum. However, in these 0(3) x 0(2)-symmetric solutions, because of the existence of the horizon 
and a nonzero temperature, very small e is not a very special point in the parameter space and 
the solutions for both transitions from true to false vacua and vice versa are possible. If e gets too 
large (either negative or positive), it loses its significance because after a while the field does not 
start close to the true vacuum and therefore does not see a very large e. 

The effect of increasing the barrier heights is quite notable. In flat space, increasing the barrier 
height increases the surface tension and, as a result, pushes the bubble radius outwards and makes 
the wall thinner. But this is not completely true for our solutions. Increasing the barrier height 
increases the surface tension and as a result decreases the horizon radius. Therefore the wall is 
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pushed inwards (in r coordinates). But if we look at the field profiles as functions of r/ru, the 
situation is more similar to the flat space case. The wall is pushed outwards and gets narrower. 
This behavior continues until the wall reaches the horizon and in this regime increasing the barrier 
heights further does not make the wall thinner. The scalar field rapidly changes near the horizon 
and (j)' develops a singularity near the horizon. However it is a coordinate singularity and none of 
the physical quantities blow up. 

For our solutions to be relevant to tunneling processes, it is not sufficient for them to satisfy the 
field equations. The fluctuations around these solutions must have one and only one negative mode. 
If there are more than one negative mode, it is an indication that there is another path in field 
space that has a lower Euclidean action and that dominates the decay. There are potentially three 
different types of negative modes for our 0(3) x 0(2)-symmetric solutions. The hrst one is the usual 
r-independent negative mode corresponding to increasing or decreasing the radius of the bubble. 
This is what we expect for a critical bubble and we expect it to always be present. The second 


type of negative modes is similar to the negative modes explained in Sec. 2.5 which correspond 
to modulating the amplitude of the expansion-contraction mode with a sinusoidal variation in 


imaginary time. As explained in 2.5 these additional modes are present when the radius of the 
bubble is much smaller than the horizon radius and therefore we have to discard these very small 
bubbles. The third type of negative modes is closely related to the four-sphere geometry of the 
bounce. The easiest way to visualize this mode is by looking at the shaded region in Fig j3.3[ If 
this region gets too narrow, we can reduce its length and therefore its action by moving along the 
r direction . This mode may be avoided if the radius of bubble is not too small and the the wall is 
not very thin. This is the same parameter region which we may need to be in to avoid the second 
type of negative modes. 

Even if these 0(3) x 0(2) symmetric bubbles have only one negative mode, they will not be very 
relevant if their action is higher than the CDL action. One of our main motivations for studying 
these solutions was checking whether the 0(4)-symmetric solutions continue to be the dominant 
solutions in curved spacetime or not. In Fig j3.1l| we plotted the action of our solutions compared 
to the CDL bounce.lt is clear from the graphs that in all of these cases the CDL bounces have lower 
action and therefore are the dominant path of tunneling. This is in agreement with the thin-wall 
approximation and brane nucleations presented in [33]. However we did not restrict ourselves to 
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Figure 3.11: Comparison of the tunneling exponents for the CDL bounce (red diamonds) and the 
0(3)x0(2)-symmetric bounce (black circles) as one parameter is varied with the other three held 
fixed. When held fixed, Uq = 2, a = 3, e = 0.3, and SttGv'^ = 0.126. 


the thin-wall limit and by scanning over a much larger portion of the parameter space, although 
not conclusively, we provided more evidence that the 0(4)-symmetric bounces are the dominant 
tunneling paths even in curved spacetime. 
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Chapter 4 

Bubble nucleation in spatial vector 
fields 


In Chapter we gave a plausible argument that even in the presence of gravity the bubbles with 
0(3) X 0(2) symmetry are subdominant and 0(4)-symmetric bubbles dominate. But all of the the¬ 
ories we studied were Lorentz invariant field theories. In this chapter we study a three-component 
spatial vector field theory which prefers non-spherical bubbles as the dominant mode of the decay. 
The model we study has a field that transforms under spatial rotations and that has different lon¬ 
gitudinal and transverse speeds of sound. This chapter is based on |l3]. First we show that the 
domain wall tension depends on the orientation of the wall and describe analytical and numerical 
methods to calculate this orientation-dependent domain-wall tension. We then show how to calcu¬ 
late the shape of the critical bubble for thermal tunneling. The tunneling exponents get modified 
by a factor of the ratio of the speeds of sound in the longitudinal and transverse directions. If this 
ratio goes beyond (a model dependent) threshold, the flat wall becomes unstable to the formation 
of zigzag segments of wall. This causes a kink in the shape of the critical bubble. The scaling of of 
the tunneling exponents with the speeds of sound is very different from what we expect for smooth 
bubbles. 
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4.1 Motivations and outline 


As mentioned earlier, most first-order phase transitions happen through thermal or quantum nu- 
cleation of a bubble of true vacuum in the surrounding false vacuum. In this chapter we focus on 
thermal nucleation of bubbles. If the bubble is too small, the surface tension makes it collapse to 
nothing. If it is too large, it expands to complete the phase transition. In this sense, we need to 
create a critical bubble which corresponds to the lowest saddle point of the energy functional with 
only one negative eigenvalue, along the direction of expansion-contraction. For thermal tunneling, 
the dominant solutions are constant in time. Therefore time variable is treated separately and in 
general easier to study. This makes the non-relativistic theories which have vectors transforming 
under the spatial rotation group interesting. Studying them may shed light on many subtleties 
which are difficult to see in other contexts. For the case of scalar fields in N spatial dimensions, it 
has been shown that the 0(N) symmetric solutions have the lowest action. The probability 

F of creating a bubble of radius R is 

where E{R) is given by 

E{R) = anN-iR^~^ - AVQnR^ , (4.2) 


where AV = Vfv — Vtv and Htv is the area of a unit N-sphere. The surface tension a depends on 
the path 7 in the field space that takes the field from false vacuum to near the true vacuum 


a = vp / d4> . 

J 'y 


(4.3) 


The scaling with the Fermi velocity vp follows from the equations of motion. For the critical bubble 


we choose the path that minimizes Eq. (4.3) and then we minimize Eq. (4.2) with respect to R. 


This leads to the commonly used estimate for the tunneling rate. 

.N 1 


logF 


a 


(4.4) 


(AV)iV-i k^T ’ 

In this chapter we will generalize the theory to include vector fields. Our motivation arises from 
condensed matter systems such as liquid crystals. Helium 3 and Langmuir monolayers [l5H52]0 We 


^In cases like the famous A-B transition in liquid ®He, the system shows different longitudinal and transverse 


speeds of sound. However, people used to study their nucleation properties using spherical bubbles. 
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study models which show non-spherical bubbles. To do so, we use spatial N-dimensional vector 
fields in an (N+l)-dimensional spacetime. These vectors transform under the spatial rotation group 
and have different speeds of sound in the longitudinal and transverse directions. The Lagrangian 
for this model is 


~ ordi(l>jdi(l>j - (cl - CT)di(j)idj(f)j ) - V{(j)i) . 


(4.5) 


The transitions occur between the two minima of the potential V, which we denote by 4>±. To 
avoid ghost-like instabilities we need cl > ct ■ When cl ^ ct, the potential can minimally break 
the spatial SO{N) symmetry. We will consider this minimal breaking. In this set up, {4>+ — (j)-) 
specifies a direction which is the direction of the longitudinal wall and causes non-spherical bubbles. 


We first study the planar domain-walls in 4.2 and set up numerical and analytical results to 
determine the orientation dependence of the tension. We show a variety of different behaviors of 
< 7 ( 6 ), the tension of the wall oriented at the polar angle 6 , in appendix [Bj An instability in the 


flat wall is studied in 4.2.2 and the corresponding effects on bubble shape are shown in 4.3 When 
Ct is not much smaller than cl, the bubble shape is close to a sphere. By increasing this ratio, 
the bubble shape gets more and more skewed and after reaching a critical value, it becomes a 


multivalued function of angle. In 4.3 we study these effects and give the correct interpretation of 
the singular shape of the bubble and find the corrections made to the tunneling exponents and 
transition rate. These techniques are identical to the one used for equilibrium bubbles, known as 
the Wulff construction [53]. It has been applied to “soft matter” systems like liquid crystals and 
Langmuir monolayers |45l H9ll52| . Our result agrees with the major conclusions in these earlier 
works. 


4.2 Orientation Dependence 


The equation of motion derived from the Lagrangian in Eq. (4.5) is 


ipi - - {cl - (^)didj4>j 


dV 

d(t)i 


(4.6) 


Here it is more apparent that ct and cl correspond to the transverse and the longitudinal sound 
speeds. To make the transitions possible, the potential V must have at least two isolated local 
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minima. There are many potentials which show this feature. For example, if a vector field cj) is put 
in an external field H, the potential to the second order in the external field H is 


m 


A 


y{4>) = + ^W + a{H • <t>) + h{H • cj^Y + c{H X <j>Y + f{HY{<i>Y ■ 


(4.7) 


2 ' 4' 

Here we assume that c = / = 0. Neglecting a and assuming 6 < 0 makes the direction parallel to 
H preferable. If + m?/2 < 0, we get two degenerate vacua at 


(t)± = 


'm2 + 26#2 H 


(4.8) 


A \H\ 

The degeneracy in the minima can be removed by allowing a small a to allow a first order phase 
transition. This example shows how generic this model is. Our further analysis in this chapter 


will be independent of the form of the potential in Eq. (4.7). Later we show other examples. We 


start by considering a potential with two almost degenerate vacua. This allows for a thin-wall 
approximation which simplihes the problem to that of finding the tension for different orientations. 
Since the interpolation between the two vacua connects two vectors in the field space, it breaks the 

4>-, which is the field at the true vacuum (interior 


spatial rotational symmetry, as shown in Fig. 


4.1 


of the bubble), is shown as red arrows and the vector field at the false vacuum (exterior of the 
bubble), is shown as blue vectors. The transition between them in different orientations corresponds 
to different contributions from gradient terms in the longitudinal and transverse directions which 
in turn makes the action along the walls (which per unit area is the tension) dependent on the 
orientation. The longitudinal ones are the most massive domain walls. The field wants to spend 
less time in those orientations. This skews the shape of bubble and creates non-spherical bubbles. 


Here we will provide the general formalism to find cr{6), and then in Sec 4.3 we will use this to find 
the bubble shape. 


4.2.1 Domain walls in two spatial dimensions 

The easiest cases to study are flat domain walls in a theory with vector fields in two spatial 
dimensions. The domain wall solutions for the case of degenerate vacua are static solutions of the 


equations of motion obtained in Eq. (4.6) 


BV 

- + dl)(l)x - (4 - cl)dxidx(j)x + dy4>y) = , 


-C^{dl + dl)(t)y - {cl - cY)dy{dx(t>X + dy(l)y) = - 


BV 

Bcpy 


(4.9) 
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Figure 4.1: The blue (longer) and red (shorter) arrows represent the vector field at the two vacua. 
The thick green line is the domain-wall. From left to right, we show a longitudinal wall, a transverse 
wall, and a wall with orientation 9. The orientation is defined such that for a longitudinal wall 
6 = 0, and for a transverse wall 9 = 7r/2. 


If the wall is along the x direction, the boundary conditions can be written as 


dx(f = 0 , 

lim (j){x,y) = $± , (4.10) 

y±^oc 

where are the true and false vacua. Similarly if the wall is parallel to v and orthogonal to u the 
boundary conditions can be written as 


< 

II 

0, 

(4.11) 

lim (f>{Xu) = 

A— 


(4.12) 


For a given wall, we can always rotate the coordinate systems (and at the same time </>) to set u to 
be along the y-axis. If the rotation is through an angle 9 and the rotated vector field is (j), we can 
define the rotated potential V0{4>x,4>y) as 


V{^x, ^y) = V((px cos 9 + (t)y sin 9, ipy cos 9 - ipx sin 9) = V0{<4)x, 4)y) ■ 


(4.13) 


In terms of this new field (f), the wall is along the x axis and therefore x-independent. After this we 
will drop the tildes for simplicity. The field profile only depends on y. This simplifies the equation 
of motion to 


- c^dlcjix 
-cld^(py 


m 

d(l)x 

m 

d4>y 


(4.14) 
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The tension of the domain wall is given by the total energy per unit x (assuming the potential 
of the vacua is set to be zero). 


a{e) = dy 




= 2 dyVe . 


(4.15) 


The last step came from integrating the equations of motion in Eq. (4.14). The practical way to 
find the domain-wall solution is to numerically minimize this tension [531 - 156] . which is what we do 
in Appendix [B| 

For a general potential, the orientation dependence of a can be complicated. Here we start 
from a simple and, in some sense, the typical case. Imagine the situation where at 9 = 0 the 
interpolation is purely longitudinal(/ij, = const. Since a rotation of 7r/2 just exchanges (l)x and 
(j)y, the interpolation will become purely transverse with = const. It is then easy to work out 


from Eq. (4.15) that 


cr(0) = CL 


f path 


^/w\d$\ , 




= CT 


path 


VW\d(P\ , 


where the two integration paths are the same, so 

O-(O) CL 




CT 


(4.16) 


(4.17) 


Potentials given by Eq. (4.7) when > 0 satisfy Eq. (4.17). So do the simpler potentials we use 
in Appendix They not only show a good agreement with Eq. (4.17), but also demonstrate an 


excellent fit to a naive interpolation, 


a{9) = cj(0) cos^ 9 + ^ > 


(4.18) 


in the regular range of parameters. From the symmetry of the problem, it seems natural to expand 
cr{9) as a polynomial of sin^ 9 and keep the lowest order terms. 

We also analyze two extreme choices of parameters in Appendix |B| One of them corresponds 


to Eq.(4.18) and the other satisfies 

(t{9) = ^ (t(0)2 cos^ 9 + ^ • 


(4.19) 


^Note that we talk about a particular solution, instead of imposing some symmetry on V. This is necessary. 
One might try a rotational symmetry on V along the vector (0+ — (f)-). That tnrns out to be neither necessary nor 
sufficient to guarantee that is constant. 
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It turns out that Eq. (4.18) and (4.19) are quite representative for our further analysis. Despite 


their simple forms which will simplify the calculation, they can actually have dramatically different 
behaviors. 


4.2.2 Flat Wall Instability 


Because cl > ct and (/>+ — cj)- is along the y-axis, it is apparent from Eq. (4.17) that the wall 


in the longitudinal direction has the highest tension o'(O) . Since the energy of a domain wall is 
proportional to its length, we would naively expect straight flat walls to have the lowest energy. 
However, if in some direction the wall gets too heavy, it may be favorable to create a wall which is 
longer (not straight-line) but with lower energy. Therefore, even if the boundary condition is set 
up to preserve the x translational symmetry, the minimum energy interpolation can spontaneously 


break that symmetry. The full treatment of this problem is to remove the condition (4.11) and see if 


a symmetry breaking configuration can further minimize the total energy. However, this approach 
requires a complicated numerical work and we will not pursue it in this chapter. We will simply 
demonstrate this possibility in the thin-wall approximation. 

The total energy is a functional of the domain-wall shape, y{x). 

^^2 ai 9 ) 


CX2 


E[y{x)] = / (T(-tan ^ y')v^l -h y'‘^dx = 


cos 9 


dx . 


(4.20) 


JXi JXi 

Here y' denotes the derivative with respect to x. Now, given a symmetric boundary condition 
y{xi) = y{x 2 ), we can study whether the flat wall is stable to perturbations that deform the wall 
and, if it is unstable, we can look for other configurations which have the lowest possible energy. 
Expanding E near 9 = 0 gives (assuming (t( 0) is a stationary point ) 


^ ^ rwid'^a 

E ~ Efiat + / ( 7 ; “^777 

J Xi 


2 




cj(0) 


B=0 


59‘^dx . 


(4.21) 


Thus, the perturbative stability condition is 

1 d^a 


cj(0) d9^ 


> -1 . 


(4.22) 


9=0 


Next, the non-perturbative instability can occur if there is a 0 7 ^ 0 such that 


a{9) < cr(0) cos 9 . 


(4.23) 
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Figure 4.2: When the flat, longitudinal wall on the left figure is too massive, it will spontaneously 
breaks into zigzag segments in the right figure. Although the total wall area increases, the reduced 
tension still reduces total energy. 


When either or both instabilities exists, there is a critical angle Oc such that (t{9c)/ cos 9c is the 
global minimum, and the wall prefers to settle into the zig-zag configuration in which every segment 
is oriented at 9c, as shown in Fig j4.2[ It is favorable to not have many corners which will cost more 
energy. It is straightforward to see that domain walls more massive than a{9c) always break into 
zigzags, while less massive walls are unaffected. 

The edges of the zigzag will contribute to the energy of the wall and the thin-wall approximation 
will break down at the sharp corners. A more careful treatment using thick-wall effects is necessary. 
However, if the distance between points xi and X 2 is large enough and there are not many kinks, 
the contribution from the corners will be negligible and the zigzag walls will have lower energy and 
still be favorable. For the case of bubbles, we will not necessarily have the possibility of making 
the energy of the kinks negligible. However if the difference between the energy of the true and 
false vacuum is small, the bubble size will become very large and we can restore the results for the 
breakdown of flat walls El 


The two examples whose cr{9) is given by Eq. (4.18) and Eq. (4.19) show very different behaviors 
for different ratios of c^/cr. While the latter always gives stable flat walls, the former develops 


^There are two different scales for nucleation of kinky bubbles. One of them is the thin-wall approximation which 
demands that the physical length in which the transition happens be small compared to the radius of the bubble. 
We can achieve this approximation by making the energy difference of the two vacua small. The other hierarchy we 
need is the smallness of the energy of the kinks to the energy stored in the surface tension of the wall. Fortunately 
this also can be achieved by creating a large bubble. Therefore, in order to make both of these approximations valid, 
we just need to make the two vacua almost degenerate. 
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kinks once 


cj(0) 


= ^> 2 . 


cr(7r/2) CT 

A flat wall with 9 = 0 wonld break and settle into zigzag segments with 


Or = sin 


-1 


la{ 0 ) - 2a(f) . 

= sin 


o-(O) - o-(f) 


Cl — 2ct 

CL - CT 


(4.24) 


(4.25) 


The wall shape and energy between two points xi and X 2 for the latter are shown in Fig. 4.3 


Energy 



Figure 4.3: The right figure is the energy stored in the wall connecting points xi and X 2 (left figure) 
vs the angle 6. The solid blue, dashed purple, dot-dashed red and dotted green lines correspond to 
clIct = 1.5, 2, 2.5 and 3. The dots show the angle with minimum energy. Clearly, for high ratios 
of cl/ct-, the zigzag walls (0min / 0) are energetically favorable. 


4.3 Bubble Shape 


In order to allow a phase transition, we have to break the degeneracy between the two vacua and 
in order to keep the thin-wall approximation valid, the energy difference between the two vacua 
should be kept small. The phase transition via thermal bubble nucleation has a rate 


exp 


hT 


(4.26) 


where Eg is the saddle point energy of the bubble. We can find this saddle point by treating E as 
a functional of the bubble shape y{x), 


E[y{x)\ = (surface integral of u) — AI/(volume) 

a (- tan“^ y) + - y AV 


A dx 


(4.27) 
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We put the origin at the center of the bubble. The bubble has a reflection symmetry along both 
the X- and y-axes. Therefore the contributions from different quadrants to the bubble energy are 
equal and we can keep only the first quadrant. In this quadrant y' = — tan 9. This configuration is 
shown in Fig. |4.4[ To minimize the energy functional, we need to specify the boundary conditions. 
The ordinary boundary conditions are y'{B) = 0 and y(xmax) = 0 at some unknown value of Xmax, 
which will be determined later by the energy difference between the two vacua. But we learned 


from Sec,4.2.2 to replace y'{t)) = 0 by y'{0) = — tan0c. Despite that y{x) is not smooth, it does 
eliminate the boundary variation and is indeed what we get from the Euler-Lagrange equation. A 
more formal argument is to write down E[x{y)] instead, for which the standard choice, x = 0 and 
x' = 0 does not exclude the kink. The resulting Euler-Lagrange equation is essentially the same 
and we can just use it. 


y 



Eigure 4.4: Bubble shape in the first quadrant. 


We will keep it simple and solve the Euler-Lagrange equation for Eq.(4.27) 

y' 


d 


^^ = Tx 


^^/l+y'^ + cr 
dy' 


\/l +y'^ 


d (.. da „ 
= — cr sm 9 + —— cos 9 
dx V d9 


(4.28) 


The general solution for x is 


x(0) = const. + ( a sin 9 + ^ cos 9 

A V V d9 


(4.29) 


The quantity in the parentheses is zero at both 9 = 9c and 0 = 0. Therefore, solutions starting at 
either value will eliminate boundary variations, as promised, and also set that integration constant 
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to zero. We can then integrate to find 

^ aV ^ ^ 

For the case of higher dimensions, the structure of the potential can further break the SO{N — 1) 
symmetry, since the exact interpolation between (p± may still involve nontrivial profiles of the 
transverse fields. But if the solutions retain SO (A — 1) symmetry the generalization to higher 
dimensions is straightforward. We can choose the XAr-axis along — </>_) and denote it by 
and xi to xn-i as the transverse directions xt- In these cases a is only a function of the polar 
angle 6 and the energy is 


E[xl{xt)\ = J dxTX^ ^ ^(T(tan ^ x'jP)^! + x'1 — AVxl^ , (4.31) 


where LIn -2 is the area of a unit {N — 2)-sphere and xt = xf. 

This leads to the general solution, 

^ (Rasing + ^cosg 

^ fo-cosg- ^sin6> 


(4.32) 

(4.33) 


Although this is a naive generalization of m, we should take a closer look. Note that by 
symmetry we have ^ = 0 at 0 and 7r/2. Also, for simplicity we can treat cr{9) as a monotonically 
decreasing function. So we can see that x^v is positive definite, but there is a risk of xt being 
negative. Since x{tt/2) is still always positive definite, and x goes to zero exactly at 6c, if we 
naively plot Eq. ( |4.33 ) for the range 6 G [0, 27r], we may get something like a wrapped candy, as 
in Fig j4.5[ We will provide a simple argument using the number of negative modes to prove the 
following statement: 


Eq. (4.33) always gives the correct critical bubble profile. When the fiat longitudinal 


domain-wall is stable, this solntion is valid for 0 < 0 < 7r/2. When the flat longitudinal 
domain-wall is unstable, we should take the largest Oc such that xt{0c) = 0 and use 
the portion Oc < 0 < 7r/2. Therefore in Fig j4.5| we only keep the middle part of the 
candy shape region. 


Any valid tunneling configuration must have one and only one negative mode, in the expansion- 
contraction direction. The energy should increase for any other perturbation. This means that 
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Figure 4.5: The bubble profile given by Eqs. (4.33) and (4.18) when cl > 2ct- The correct prohle 


of the critical bubble is simply the middle portion. 

fluctuations of the wall shape should still correspond to positive modes and every segment of the 
wall should still settle into a minimum energy conhguration. However, in the contraction-expansion 
direction we must have a maximum energy. Let’s go back to the (2-|-l)-dimensional case and pick 


two pieces of wall and glue them together. Since these pieces satisfy Eq. (4.33), they are indeed 


local minima with respect to fluctuations that change the shape of the wall. However, to get one 
negative mode, at the location of gluing, we must have a maximum of energy for the critical bubble. 


This situation is shown in Fig. 4.6 We see that, at least for the case that satisfies Eq. (4.18), for 


large values of c^/ct this angle is not at zero and therefore there are kinky bubbles. To make 
this more rigorous, let’s look at Fig j4.7[ We can glue the two pieces, as in the right panel, using 
a smaller portion of the two shells. Also, can we take the two shells further apart and interpolate 
between them with zigzag walls. The former possibility is making the bubble smaller, while the 
later is making it bigger. Through a pictorial argument, we can show that they both make the 


total energy smaller, establishing that Eq. (4.43) is really the saddle point with this unique negative 
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E(0) 



Figure 4.6: Picking two pieces of the wall (top left panel) and gluing them at angle 6 (top right 
panel). Each piece locally satisfies Eq. ( 4.33| ) and therefore gives the minimum energy for the 
fluctuations parallel to the wall. The energy of the bubble vs the gluing angle 6 is shown in the 


lower panel for different values of clIct for a tension that satisfies Eq. (4.18). Blue dashes, green 


dots, red dot dashes and solid purple lines correspond to c^/ct = 1.5, 2,3 and 4. The correct spot 
for gluing the two is the angle which gives a local maximum (negative mode). Eor c^/ct < 2 this 
gives a smooth bubble. For larger ratios it gives a kink. 
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Figure 4.7: The left figure depicts Eq. (4.35), where we attempt to make a larger bubble by inserting 
true-vacuum regions and extra interpolation walls. The right figure depicts Eq. (??), where we try 
to make a smaller bubble by removing part of the walls and part of the true-vacuum region. Both 
result in smaller total energy, which shows that the kinky shape is indeed a saddle point. 


mode. Erom Eq. (4.33) we get 


AVxLiOc) = {N - 1) 


crjOc) 
cos 9r 


(4.34) 


which means that the energy difference due to volume for a cylinder, an (A^ — l)-sphere times height 
xl{9c), is equal to the energy in the domain wall that covers the {N — l)-sphere by a zigzag profile 
with orientatiorl^^c- 


Then, as shown in the left portion of Fig. 4.7, the energy lost due to the green (shaded) region 
is equal to the contribution from the dotted domain-wall. After these cancel each other, the two 


"^There are two different scales and approximations involved in this problem. The first is using a flat-wall approxi¬ 
mation to derive the orientation dependence of the tension. The second is the ratio of the contribution of the kink to 
the wall in the bubble energy which we neglected here. Both of these approximations are valid if we choose a small 
energy difference between the two vacua [for example by choosing a very small a in Eq. ( |4.7[ )] which makes the bubble 
bigger and flatter and, by increasing the area of the bubble and keeping the kink contribution unchanged, makes the 
kink contribution negligible. 
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extra triangular regions still contribute —AV, so the total energy is indeed less. 

(zigzag wall) = Al/(rectangle) < Al/(extra false vacuum region) . (4.35) 

In the right portion of Fig |4.7| we try to make a smaller bubble by removing the true vacuum 
region and domain walls covered by the green (shaded) rectangle and then patching the remaining 
two shells together as a smaller bubble, with a kink angle larger than Oc- Since cr{9c)l cos 9c is a 
minimum of a{9)/ cos 9, losing those wall segments overcompensates for the energy gain even if we 
remove —AV of the entire green (shaded) rectangle. 


4.3.1 Smooth Bubbles 


We can get more intuition by solving the exact bubble shape for a specific cr{9). In the first 


example we will use Eq. (4.19), where no spontaneous symmetry breaking should occur. Thus, we 


are expecting a smooth bubble. Plugging into Eq. (4.33), we get 

(T(7r/2)^ sin0 


xt{9) = 
xl{9) = 


N - 1 


^1^ \/ ct(0)2 cos^ 9 + sin^ 9 

N — \ (t(0)^cos6< 


\J ct(0)2 cos^ 9 + (j{'K 1‘lY sin^ 9 
Obviously, the bubble takes the shape of an ellipsoid. 




+ ^ = (iV-l)^r^, 


where 


cj(0) (T(7r/2) 


1 




comes from Eq. (4.16). 


It is then straightforward to calculate the saddle point energy. 


Es = S'at-; 


(A-l) 


A^-1 


path 


N 


vwi#!)" 


(4.36) 

(4.37) 

(4.38) 

(4.39) 


(4.40) 


Comparing this to the usual form people use assuming a spherical bubble, Eq. (4.1), the difference 
can be characterized by an effective Fermi velocity. 


VF (CLC^ ^)^/^ , 


(4.41) 


which is a weighted geometric average of sound speeds. 
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4.3.2 Kinky Bubbles 


Now we turn our attention to Eq. (4.18). Plugging it into Eq. (4.33), we get 

(iV-1) 


xt{0) = 

xl{0) = 


AV 

(N-l) 

AV 


2o-(7r/2) - cj(0) 
2o-(0) - cj(7r/2) 


sin0 + 
cos 6 — 


(t( 0) — cr(7r/2) 
it(0) — a{Ti/2) 


sin^ 6 
cos^ 9 


(4.42) 

(4.43) 


As expected, when cr(0) < 2o{jxjT)^ we still have a smooth bubble profile. When (t( 0) > 2cr(7r/2) 


we just use the portion 9c < 9 < 7r/2. This is proved in Sec,4.3 


The expression for Eg is quite complicated in arbitrary dimensions, so we only present the 
results for “realistic” dimensions. For N = 2, we have 


E. 


N=2 


-^ (lOCLCT - cl - cl) - , for CL < 2ct , 


(4.44) 


AAV 


(IOc^ct — cl — cl) cos 


-1 


Cl - ‘^ct 
CL - CT 


+ {cl + 13ct)-\/(cl^^^ct)ct^ , for cl > 2ct 


(4.45) 


This is quite complicated. We should again compare it to the spherical bubble and think in 
terms of the effective Fermi velocity, especially in the limit cl 3> ct- 


VF 

Vp 


^ lOcLCr - cl- c l'^ 


for Cl < 2ct , 


— {cLcl)^^"^, for cl:>ct 


(4.46) 


For N = 3, we have 


eI^=^ = Af 


= 47r 

Similarly we have 


path V^\d(j)\^ 4 |-g 2 _ iOc|cr + 52clcI — 8cl) 

AF2 105 

path 32 cI(7cl - 6ct)^ 


for Cl < 2ct , 


AF2 


105Vcl - Ct 


(8 1/2 5/2\ r _ 

( Cl:^ CT . 


, for Cl > 2ct ■ 


VF 


(4.47) 


(4.48) 


Comparing these to Eq. (4.41), we find that the effective Fermi velocity is still a weighted 


geometric mean, but the weight on is always reduced by half. This is quite understandable since 
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the critical bubble approaches a thin slit. A major portion of its domain-wall is aligned in the 
transverse direction. It is straightforward to show that this limit generalizes to N dimensions as 

vp ~ . (4.49) 


4.4 Summary 


In this chapter we studied the vacuum decay problem in a vector field theory which has different 
longitudinal and transverse speeds of sound and constructed the critical bubbles for thermal tunnel¬ 
ing. Because the spherical symmetry is broken due to a combination of different sound speeds and 
having a singled out direction ((/>_|_ — </>_), the surface tension has a dependence on the orientation 
of the wall. This makes the critical bubbles non-spherical. We derived the shape of the bubble 


using a simple formula, Eq. (4.33), and showed a modification in the bubble energy by a factor of 
speeds of sound ratios. Since the tunneling rate depends exponentially on this energy, this implies 
a large modification to the rates. We showed that if the ratio of speeds of sound goes beyond a 
critical value, the flat walls become unstable and break spontaneously into zigzag segments. The 
same effect is seen for the nucleation of bubbles, and so the critical bubble develops kinks. We gave 
a recipe for making these non-smooth bubbles and explained the correct way to interpret them 
according to the number of negative modes. We also developed numerical ways to calculate the 
orientation dependence of the tension and showed examples for different behaviors of the surface 

tensiorl3 

Our analytic and numerical study shows that the freedom to take different paths in a multi¬ 
dimensional field space is essential for the instability. If we choose parameters that reduce the 
number of dynamical fields down to one, the longitudinal wall is always stable. That is however 
an extreme choice. For typical choices of parameters, at cpjcT > 2 the longitudinal wall becomes 


^Analysis of the bubble shape for nucleation is identical to that for equilibrium bubbles, known as the Wulff 
construction. Earlier works |45H52 | have qualitatively similar results. We further show that the flat wall settles into 
zigzag segments with a corner angle of 9c- The recognition and interpretation of the kinky bubble shape is also more 
transparent in our analysis. The tension previously studied is often expanded as o{9) = no + a cos 9 + & cos 20, and 
most analysis focused on the effect of a 7^ 0. In our model there is a reflection symmetry—the domain-wall tension 
does not change when you look at it from the other side. Thus we always have a = 0. This makes our situation closer 
to a two dimensional lattice model[^, where similar bubble shape was observed. 
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unstable, and the critical bubble develops two kinks. Such behavior can appear with an even 
smaller sound speed ratio, clIct > \/2, if we tune the potential to the other extreme limit. This 
range of sound speed ratios is not hard to find in real materials. 


We picked two representative forms of a{0), given by Eq. (4.18) and (4.19), to calculate the 
exact shapes of critical bubbles. This allowed us to observe the scaling property of tunneling rates. 
When we further increase the sound speed ratio the critical bubble is deformed but still smooth, 


we can modify the standard tunneling rate formula, Eq. (|4.1|), through an effective Fermi velocity, 

(4.50) 


VF (CLC^ l)jl/iV _ 


This is quite intuitive, since one particular orientation is longitudinal and all others are transverse. 
They care about the sound speeds in their own orientations. When we increase the sound speed 
ratio the bubble starts to develop kinks and the scaling behavior changes to 


/ i/Z I 

Vf ^ (c^ c. 


1/2 {N-1/2) 


)l/iV _ 


(4.51) 


This is because the kink development removes a large portion of the longitudinally oriented wall, 
so Cl becomes less important. 

On top of modifying the tunneling rate estimation, our result has a practical impact. Typically, 
experimental measurement of domain-wall tension involves measuring the bubble radius |58]. That 
is done by observing a domain wall popping through a partition with holes. When it does, the 
radius of the hole is identified with the bubble radius. Our result shows that for vector fields, 
the orientation of that partition is important. The popping radius can be identified with x given 


by Eq. (4.33) only for a longitudinally oriented partition. For other orientations the hole and the 
bubble do not have common symmetries. Therefore the exact relation between the popping radius 
and the critical radius requires further analysis. 

We have only taken a small step toward a rich phenomenology. Given the new insight here, 
many nontrivial questions arise. How does the domain wall move (bubble expand) given this 
orientation dependence? When there is a kink, can we expect the tip to travel at cl, leaving 
behind a Cherenkov-like tail of domain walls bounded by c^? How does the spontaneously broken 
planar symmetry interact with impurities or other external effects? All these await future study, 
and may lead towards a more practical understanding of some exotic theories of phase transitions 
that rely on the properties of domain walls [5^ 159116(1]. 
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Chapter 5 

Instabilities in large dimensional field 
spaces 


There are important qualitative and quantitative differences between tunneling in theories with one 
scalar field and ones which involve many fields. We studied one specific case in Chapter that has 
more than one field and found dramatically different behavior. In general, finding bounces solution 
in multi-field theories is not easy and there is not much focus on it in the standard literature of 
tunneling. But there is one aspect of these theories which on very basic principles may be worrisome 


and needs further exploration. From Fig. 5.1, it is clear that, in the case with one field, to get from 
the false vacuum to the true vacuum we must pass through the highest barrier that separates them, 
but that for a theory with two fields there are many other passages which do not require crossing 
the highest barrier. This may enhance the tunneling process. The number of saddle points and 
possible tunneling paths increases very rapidly with the number of fields. This may render most 
of the minima of these theories unstable and short-lived due to quantum or thermal fluctuations. 
In this chapter, which is based on m, we show that, indeed, a generic minimum is not protected 
by high enough walls to ensure its stability. Specifically, we will study theories with multiple 
scalar fields and provide numerical evidence that for a generic local minimum of the potential the 
semiclassical tunneling rate, F = Ae~^ increases rapidly as the number of fields increases. As 
a consequence, the fraction of vacua with rates low enough to maintain metastability decreases 
exponentially with the number of fields. We also discuss possible implications for the landscape of 
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string theory. If our results prove applicable to string theory, the landscape of metastable vacua 
may not contain sufficient diversity to offer a natural explanation of dark energy. 



Figure 5.1: In a theory with one scalar field (left picture), to go from the false vacuum (right 
minimum) to the true vacuum (left minimum) we have to cross the highest barrier. But for the 
case of a field theory with two fields (right picture), to go from the false vacuum to a truer one, we 
do not need to pass the highest barrier (central peak) and can choose to go through many other 
lower barriers. This decreases the instanton action and enhances the tunneling rate. 


The discovery that string theory admits a huge number of flux vacua |62H64j played a very 
important role in the development of the theory over the last decade. The landscape of vacua of 
theories with several hundred moduli fields was considered to be a good anthropic framework to 
address the problem of naturalness of the cosmological constant. Many aspects of these vacua, from 
their phenomenological and cosmological properties to their distribution and statistical features, 
have been extensively studied [SSI [SS] . There have been many attempts to study some of the 
statistical features of these vacua in previous works [67H71]. But due to the complexity of this 
landscape of vacua, many important questions have not been answered. One important one is 
whether the vacua obtained in this way are long-lived or no. A direct analysis of this question brings 
formidable challenges in calculation. Therefore, we take a different approach and study generic field 
theoretic models of landscape and focus on how the stability of vacua varies as the dimension of the 
moduli space (the number of fields) increases. Our results suggest that tunneling rates, and hence 
vacuum instability grow so rapidly with the number of moduli fields that the probability of a given 
local minimum being metastable is exponentially small. At the semiclassical level, the transition 
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rate is given by F = Ae~^ where B is the instanton action and A is given by the determinant of 
fluctuations around the classical path. We show that B drops very quickly with the number of 
moduli fields. For tunneling from de Sitter, we also have the possibility of Hawking-Moss solutions. 
The tunneling exponents in the latter case in the thin-wall approximation are proportional to the 
difference between the saddle point energy density and the energy density of the false vacuum. We 
show that this difference drops very quickly with the number of moduli fields, so this tunneling 
channel also gets a big enhancement. In string theory, the complexity of the geometry and the large 
number of moduli fields makes the situation much more complicated. Meanwhile, it opens up a new 
line of attack based on statistical analysis. For example, Douglas and Denef m developed a method 
of calculating the density of flux vacua in the string landscape in terms of the Kahler potential on 
the moduli space of a given Calabi-Yau compactification. Their work showed that the vacua tend to 
accumulate near the conifold locus in moduli space. Dine et al.[69] used a scaling argument to show 
that the vacua with small cosmological constant in the landscape of string theory become unstable 
when the fluxes are large compared to their compactification volume. Chen et al. m studied the 
chance that a given stationary point in the landscape is a minimum. They found that the chance 
that a randomly chosen saddle point is a minimum is suppressed exponentially for large numbers of 
moduli fields. They showed that the suppression becomes more severe as the cosmological constant 


of the vacuum gets higher. In Sec. 5.1, we discuss our approach for estimating tunneling rates in 


a field theoretical context. In Sec. 5.2, we show the numerical methods and results. These results 


show a general feature of high-dimensional field theories and are independent of the application to 
string theory. In Sec. |5.3[ we discuss the possible implications of our results for the landscape of 
string theory and explain some of the considerations that would need to be resolved before applying 


this model to the string landscape. Finally, in Sec. 5.5, we estimate the maximum dimension of 
moduli spaces that allows enough diversity to explain the cosmological constant in a natural way, 
and we suggest future directions for studying these models. 


5.1 Background and approximation method 


We consider the dynamics of a system composed of N scalar fields (pj with a Lagrangian 


1 


N 


C = X ^- I/((/>i,(/>2, ...An) . 


i=i 


(5.1) 
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In general the potential V has many local minima that represent metastable false vacua. We are 
interested in the rate at which the transitions between these vacua occur. Let’s consider one of 
these minima and use the freedom to shift the field by a constant value to bring this minimum to 
the origin. Since we are ignoring gravity at this point, we can shift the value of the potential by 
a constant number and set it to zero at the origin. If the potential is smooth enough around this 
minimum, we can expand the potential in a power series, 


(5.2) 


y = A ^ AjkMj^kv + Y, + ■ 

\ i ijk ijkl 

Here we used the freedom of using an 0{N) transformation to diagonalize the quadratic terms. 

( 2 ) 

Assuming that the origin is a minimum makes AI positive definite. Here we have extracted a mass 
scale V from the potential to make and dimensionless. This scale determines the typical 
distance between the stationary points of V in the field space. Finally, we have extracted an overall 
factor A from the potential to make the quantities inside the parenthesis of order unity. Although 
the potential itself is bounded from below, if we keep a finite number of terms in the series, 
the resulting approximation may not be so. We are not concerned about making the potential 


in Eq.(5.2) bounded form below as long as we care about the immediate vicinity of the origin. 
Therefore the coefficients A^^j, and can be positive or negative without imposing an overall 
positivity condition. 

The tunneling exponent B is the Euclidean action of the bounce solution. Since we are ignoring 

gravitational effects, we can safely assume that the dominating bounce has an 0(4) symmetry and 

the fields are all functions of s = The equation for the bounce solution |I6] is 

dV 

d<j)j ■ 


(r4>j 3 d(l)j 

ds'^ ^ s ds 


(5.3) 


This configuration will have a region on the true vacuum side at the center of the bubble and, 
in order to ensure it has a finite action, at large s it should approach the false vacuum, i.e. 
lims^oo </>(s) = </>fv We do not specify (/>(0), the field at the center of bubble, a priori and since this 
is a second order differential equation, we would need another boundary condition. We demand 
</>^(0) = 0 to avoid a singularity at s = 0. Except for the thin-wall approximation, we do not expect 
the field at the origin be very close to the true vacuum. 

For our study, we will generate a large ensemble of potentials with random coefficients. In order 


to calculate the tunneling rates, we would have to solve Eq. (5.3) for each sample. But these are N 
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coupled nonlinear differential equations and solving them numerically in a reasonable time is not 
possible. Instead, we develop an approximation to make this problem tractable. Our starting point 
is the thin-wall approximation for a single scalar field. In the thin-wall case, the field configuration 
can be thought of as a ball, filled inside with true vacuum while the outside is in false vacuum. 
The two phases are separated and patched smoothly by a transition region which is thin compared 
to the length scale in the problem (the radius of the ball). In this transition region, the potential 
barrier and the gradient terms have an energy density per unit of area which is equivalent to a 
tension. Therefore the net tunneling exponent is [2H] 


B = 2M 


dss^ + 


= 271^R^a - 


(5.4) 


where the factors of vr^ come from the volume elements in four-dimensional spherical coordinates, 
R is radius of the bubble, e = Vfv — I4v comes from the integration over the inside of the bubble 
and a represents the surface tension. By virtue of the equations of motions, a can be rewritten as 


a = 


dct>^2 [V{ct>) - ^fv] , 


(5.5) 


J Nv 

and the integration is taken from the false vacuum outside the bubble to some point on the true 
vacuum side of the barrier which has the same potential as the false vacuum I4((/>*) = I^v For the 
case of multi-field theories, the possible constitute a co-dimension one subset in the field space. 
Depending on the endpoint, we will have different values for tension. 

To find the dominant path, we minimize this action with respect to R (and also the endpoints 
) to get 


R 

B 


3(T 


e 



277rV^ 


(5.6) 


These results are only valid if e is much smaller than the energy scales involved and therefore the 


tunneling rate given by Eq. (5.6) is very slow. Also this result is very sensitive to e which cannot be 
correct outside the thin-wall limit, because the field does not get close enough to the true vacuum 
to see the real energy difference. But we can extract some useful information for the thick wall 


case. 
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Here we follow chapter 12 of [28] to get an approximation for the tunneling exponent in the 


thick-wall limit. Let’s rewrite the potential in Eq. (5.2) as (notice that in this case Efv = 0) 


E(0) = \vW{(t)/v) + Efv 


We can define a dimensionless length and field as 


(5.7) 


s = 


uVA ’ 

= vd> . 


(5.8) 


In terms of these new variables, the tunneling exponent in Eq. (5.4|) becomes 

2 


97r^ /-oo 

B = ^ dss^ 
A Jo 


mi 


(5.9) 


Therefore the tunneling exponent is independent of v and inversely proportional^ to A. The bounce 
solution follows a path 71 in the field space starting from some point on the true vacuum side of 
the barrier and evolving into the false vacuum. At some point Stop, it passes halfway through the 
barrier that separates the two vacua. Similarly to the case of a thin wall, we can assign a surface 
tension to this path. 

d = j #^2 [V{(l)) - Efv] = VXv^ #-^2 [l/(0) - Efv] = VXv^s . (5.10) 

Since the contribution from s = 0 to s = Stop on average is the same as contribution from s = stop 
to the point where it reaches the false vacuum, we can evaluate the integral over a new path 7 ' 
which is the portion of 7 that extends from s = 0 to s = Stop • This causes a factor of 2 to appear 
in front of the integral to compensate for the part that we have neglected. 


a = 2 [vm - ^fv] = 2 \/Au3 / #-\/2 

J'y' Jy 


Vm - Efv 


= VXv^s . 


(5.11) 


Outside the thin-wall limit, we cannot talk about the radius of the bubble, since there is no sharply 
defined boundary between the true and false vacuum. However we can introduce a notion of 
radius to be corresponding to stop ■ This quantity which should be proportional to R = stop = 


^Although B is independent of v, the value of v effects the prefactor in the tunneling rate. 
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^top ( \/A' 




fc(VAf) ^ where A: is a number of order unity. Now we can generalize the tunneling 


exponent for the thin-wall solution calculated in Eq. (|5.6|) to the case of a thick-wall solution, 

R 


\/A Vv 


B ~ n'^R^a ~ . 

A 


(5.12) 


We can check the reliability of this approximation using a one field example which has thick-wall 
solutions. We start from a potential of the form|^ 




A 


— — 2v^(j) — 2v^) 

3 


0 < 6 < 1 


(5.13) 


This potential has two local minima at 4> = and a local maximum at 4> = —bv, while 6 is a 
measure of the breaking of Z 2 symmetry between the two vacua. If 6 « 0, the two vacua are almost 
degenerate and the thin-wall approximation works very well. The potential for different values of 
b and the corresponding exact bounce solutions are shown in Fig j5.2[ For large values of b, the 
wall gets thick and we can see how well the approximation in Eq. ( |5.12[ ) works. To do so, first we 
calculate the exact tunneling exponent by solving the bounce equation and then calculate s from 


Eq. (5.11) by evaluating the integral in the field space from (j) = —v to <f> = —bv. We then calculate 


the corresponding value for k. The values of k are shown in Table 5.1 There is not much change 
in k when we scan through the values of b. A typical value of k is some number between 5 and 


6 . Therefore the approximation introduced in Eq. (5.12) is valid and we will use it in all of our 
numerical calculations. 

Now we can use a prescription to estimate the tunneling action S in a simple form. We choose 


a number between 5 and 6 for k and use the s introduced in Eq. (5.11) as a measure of the tension 
for different paths. Minimizing the action is equivalent to minimizing the surface tension among 
all different paths 7 ' in the field space. However, this is still too complicated if the dimensionality 
of field space, and hence the number of possible choices for path 7 ', are large. But we expect that 
the paths that give the minimum tension should pass through low barriers which are the saddle 
points of the potential. As an upper bound for the tension, we can look at all the paths that go 
along straight lines from the false vacuum to a saddle point of the potential and evaluate the value 


^This example and the related graphs are taken from chapter 12 of |28 |. 
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V/Aa^ 



0/fl 



Figure 5.2: Left panel shows the potential introduced in Eq. (5.13). The progression from top to 
bottom shows the potentials corresponding to 6 = 0, 0.2,0.4,0.6 and 0.8. Case 6 = 0 is a symmetric 
double well and the larger values show more and more deviation from the Z 2 symmetry. The right 
panel shows the corresponding field profiles (top to bottom) of bounce solution for the values of 
b = 0.4,0.5,0.6, O.Tand 0.8. The solution deviates more and more from a thin wall and at b = 0.9 
it merges to a Hawking-Moss instanton. 


Table 5.1: Different values of k calculated over a wide range of potentials which have a thick- 
wall solution. Despite all the changes in the parameters of the potential, the values of k do not 
change much when we scan potentials that give nearly thin-walls solutions to the ones that give 
Hawking-Moss instantons. 


b 

0.4 

0.5 

0.6 

0.7 

0.8 

a 

0.2746 

0.1758 

0.1017 

0.0500 

0.0183 

B 

30.199 

13.973 

7.017 

3.576 

1.735 

k 

6.04 

5.42 

5.17 

5.23 

5.75 


of s along these lines. We keep the lowest of these tensions s = Smin- As a measure of tunneling 
exponent, we can use Eq. (5.12[) (setting A: ss 5) 


B ~ 10" 


A 


(5.14) 


Of course this gives an upper bound for B. The actual value should be lower than this. If the 
lifetime of the metastable vacua calculated in this approximation is not long, it should be even 
shorter when we do not use the approximation. If the tunneling exponents get too small, let’s say 
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of order unity or less, then too many bubbles pop out and this make the dilute gas approximation, 
which is a pillar of the semiclassical approximation, invalid. In this case, the metastability of the 
false vacua is essentially gone. 


5.2 Numerical results 


We numerically studied large ensembles of theories with potentials given by Eq. (5.2) with random 


coefficients. We use s as a measure of stability of the vacua. Since the numerical calculations are 
very time consuming, we only keep terms up to quartic order in fields. We later use different orders 
in the expansion and show that our results do not change much. An ensemble is defined by choosing 
random numbers for the with uniform distributions over ranges defined by 


A, 


( 2 ) 


G [ 0 ,02] , 


ll 

e [-03,03] , 


G [-04,04] , 


(5.15) 


where the A^^) are chosen to be positive to ensure that the origin of the field space is a minimum. 
The coefficients of the cubic and quartic terms are not necessarily positive and therefore the poten¬ 
tials are not in general bounded from below. However, because the important region for tunneling 
is a small neighborhood of the origin, the lack of a lower bound at large field values does not concern 
us and the boundedness of the V at large distance is achieved by higher order terms. 

Because we are interested in the dependence of the tunneling exponent on the number of fields. 


it is important to have a prescription for how to change the ranges for 02,03 and 04 in Eq. (5.15) 


when we change N. Our guiding principle is to keep the variation of the potential inside an iV- 
dimensional ball of radius 4>r independent of the number of fields. This has the advantage that if 
we constraint ourselves to an Nq < Wdimensional hypersurface of the field space, we can recover 
the same variation inside an A'o-dimensional ball of the same radius (pR. This ensures that our 
normalization does not have a peculiar dependence on the number of fields. 

With this assumption, it is an easy job to find the dependence of the on the number of fields. 
The typical variation of each field cpi inside a ball of radius cpR is (pn/^fN. There are N positive 
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quadratic terms and their sum is roughly Therefore the variation of the quadratic terms in 
the potential is independent of N if 02 is. Similarly, there are of order N^ cubic terms. Because 
they can be positive or negative, their sum is effectively N^/'^ times a typical value of which is 
(j)^/N^^‘^. Therefore we need to keep 03 independent of N. Similarly there are N^ quartic terms. 
Because they may be positive or negative, effectively we have N'^ terms. Each of them is roughly 
(p'j^/N'^ and therefore 04 should be independent of N. 

It is important to notice that the range of all the coefficients are not independent and they can 
have trivial effects. For example, if we multiply the potential by an overall factor, the ranges of 
the Uj get rescaled by the same overall factor. Also, by redefining the fields by a multiplicative 
factor, we can set the ratio 03/02 to be any arbitrary number. Therefore if we keep only up to 
quartic order, the only independent range will be the ratio of 04/03 (or the ratio of any chosen pair). 
Therefore, without losing generality, we can assume 02 = 03 = 1. For our numerical calculations, 
we also set 04 = 1. We get back to this shortly. 

For a given N, we generated an ensemble of 10000 random potentials and found all the saddle 
points.For each potential, we calculated the lowest s (as explained in the Sec j5.1[ ), the lowest barrier 
height, the distance to the lowest saddle point and the distance to the saddle point which has the 
smallest s. In Fig j5.3[ the median values of the lowest s and the lowest barrier heights are shown 
for quartic potentials. As the graphs indicate, the medians drop very quickly with the number of 
fields and they follow a power law relation. Not only do the barrier heights and the domain-wall 
tensions drop very quickly, these saddle points get closer to the origin (false vacuum) and this is a 
good indication that our truncation at the quartic level was the right idea. In Fig j5.4[ the distance 
to the saddle points with the lowest s and the saddle points with lowest barrier height for a quartic 
potential are plotted on a log-log plot. It shows very good agreement with a power law relation. 
The two graphs are almost identical. 

To understand this behavior better, in Fig |5.8t the median value for the number of saddle points 
around the origin are plotted on a log graph. Apparently this fits very well to an exponential. 
Although this is a numerical work, there are rigorous proofs that the number of stationary points 
of a polynomial of order D in N variables is an exponential in N [771178]. Now let’s look at the 
lowest barrier heights. We kept the variations of the potential independent of N (this was the 
guiding principle for choosing the If there are K saddle points and they are randomly 
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Figure 5.3: Left panel, the median value of s for quartic potentials on a log-log plot. Right panel, 
the median value of the lowest barrier heights for quartic potential in units of Au^. The bars in 
both cases show the 25th and 75th percentile. 
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Figure 5.4: Left panel, the median value of the distances to the saddle point which has the smallest 
s for quartic potentials on a log-log plot. Right panel, the median value of distances to the saddle 
points with the lowest barrier height for quartic potentials in units of v~^. The bars in both cases 
show the 25th and 75th percentile. 
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distributed (in fact being well scattered is enough and we do not really need randomness here), 
we expect the lowest of them to have a height proportional to K~^. Therefore if the number of 
saddle points K, increases rapidly with the number of fields (in this case an exponential growth), 
we expect the lowest of them to fall very rapidly with N. A similar argument may be used to 
justify that the wall tensions should drop rapidly with the number of dimensions. 
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Figure 5.5: The median value of the number of saddle points around the origin on a log plot for 
quartic potentials. As is apparent from here, this number fits very well to an exponential. The 
bars show the 25th and 75th percentile. 


From these figures, we see that the median value of all the quantities we are interested in drops 
like a power law for randomly chosen quartic potentials. Let’s define this power to be a. For 
example 

'^median — C*tensionA^ tension _ 


The values of a and C for s, lowest barrier heights, and distance to the saddle point with lowest s are 
shown in Tablefor the case of randomly chosen quartic, cubic and supersymmetric potentials 
(they are explained in the next few paragraphs). We may question the validity of cutting the 
potential at quartic order. To have a check, we repeated the same calculation for randomly chosen 
cubic potentials. For each dimension, we chose ensemble of 10000 randomly chosen cubic potentials. 


and using the prescription of Sec. 5.1, we calculated s, the lowest barrier heights and their distances 
to the false vacuum. The medians of the different quantities are shown in Figs|5.6[|5.7|and|5.8[ The 


same power law decrease is observed in the median values of the different quantities of interest. 
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Figure 5.6: Left panel, the median value of s for cubic potentials on a log-log plot. Right panel, 
the median value of the lowest barrier heights for cubic potentials in units of Au^. The bars in both 
cases show the 25th and 75th percentile. 


H 

■s. 


1 

0.48 

0.24 

0.12 

0.06 

0.03 


H 

■s. 


1 2 3 456789 10 

N 


1 

0.48 
0.24 
0.12 
0.06 
0.03 

1 2 3 456789 10 

N 



Figure 5.7: Left panel, the median value of the distances to the saddle points which have the 
smallest s for cubic potentials on a log-log plot. Right panel, the median value of the distances to 
the saddle points with lowest barrier heights for cubic potentials in units of v~^. The bars in both 
cases show the 25th and 75th percentile. 
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Figure 5.8; The median value of the number of saddle points around the origin on a log plot for 
cubic potentials. As is apparent from here, this number fits very well to an exponential. The bars 
show the 25th and 75th percentile. 
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Figure 5.9: Left panel, the median value of s for supersymmetric potentials on a log-log plot. Right 
panel, median value of the lowest barrier height for supersymmetric potentials in units of Au^. The 
bars in both cases show the 25th and 75th percentile. 
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Figure 5.10: Left panel, the median value of the distances to the saddle points which have the 
smallest s for supersymmetric potentials on a log-log plot. Right panel, the median value of 
distances to the saddle points with lowest barrier height for supersymmetric potentials in units of 
V. The bars in both cases show the 25th and 75th percentile. 
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Figure 5.11: The median value of the number of saddle points around the origin on a log plot for 
supersymmetric potentials. As is apparent from here, this number fits very well to an exponential. 
The bars show the 25th and 75th percentile. 
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Another issue to address is the correlation between the coefficients. We chose all the coefficients 


in Eq. (5.2) as uncorrelated random numbers. We may worry that the results we found about the 
rapid drop of different quantities versus the number of fields could be merely an artifact of choosing 
so many random uncorrelated numbers. To show that this is not the case, we started from the 
superpotential of a set of N scalar superfields. This potential is described in details in Appendix. 
[Pj We chose the coefficients mij,gijk and hijmk to be randomly chosen numbers. Using an SU{N) 
rotation in the field space, we can diagonalize the positive definite matrix rriij. We chose these 
coefficients to be uniformly distributed over the following intervals 


rriij 

G 

[0, 02] 

9 ijk 

G 

[0, 03] 

Qijk 

G 

[0, as] 

hijf^i 

G 

[0, 04] 

Im 

G 

[0, 04] 


(5.16) 

(5.17) 


Although we chose uncorrelated numbers for the rriij, gijk and hijki, we can see from Eq. (D.12), 
that the coefficients in the expansion are highly correlated. Please notice that the supersymmetry 


is broken after we truncated Eq. (D.IO) at the quartic order and this example only addresses the 


effect of correlation between coefficients on the median of different quantities. The median values 


of the quantities of interest are shown in Eigs. 5.9, 5.10, 5.11 Because the N chiral superfields 


correspond to 2N real scalar fields, instead of writing the fits to the data as in Eq. (5.16) we write, 

e-g-, 

^median = Ctension(2iV)-" . (5.18) 


Again we get a power law decrease for all the median quantities, but the powers are different. The 
drop with the number of fields is even more severe in the case of supersymmetric potentials which 
have correlated coefficients. 

Now we address one more issue regarding the numerical calculations. We chose the range of the 
coefficients a* to be all equal to one. As mentioned earlier, the values 02 and 03 can be arbitrarily set 
by an overall factor for the potential and rescaling the fields. We studied the effects of changing 04 
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^tension 

Q^height 

Q^distance 

Cubic potentials 

2.73 

3.16 

1.15 

Quartic potentials 

2.66 

3.12 

1.10 

SUSY 

3.16 

3.99 

1.19 



extension 

C^height 

C^distance 

Cubic potentials 

0.26 

0.090 

0.67 

Quartic potentials 

0.22 

0.083 

0.60 

SUSY 

0.25 

0.11 

0.60 


Table 5.2: Best fit parameters, defined as in Eq. (5.16), for power law fits to the data in Figs. 5.3 
and 


5.4 


5.6 


5.7 


5.9 


5.10 


on the median value of s in Fig 5.12 Making the value of 04 larger makes the tunneling exponents 
and barrier heights smaller and the limit of 04 —)• 0 simply reduces to the results we had for cubic 
potentials as expected. Therefore we cannot compensate for the drop of median values by changing 
the value of 04 (or 02 and 03 ). 

Until now, we have seen that the median values of tensions, barrier heights and tunneling 
exponents drop with a power law of the number of dimensions. Therefore if we choose a generic 
vacuum in these theories, it will not have a long lifetime. This by itself does not tell us about the 
lifetime of the “non-generic” points and the number of stable points. To find the population of 
these outliers, we studied the distribution of tunneling exponents, barrier heights and the distances 
to the false vacuum around their median values for different dimensions. In Figs. |5.13 5.14 and 


5.16 we show the distribution of values of s, Uop and saddle point distances for several choices of 
N with quartic, cubic and supersymmetric potentials. They roughly coincide when plotted as a 
functions of s/smedian- Our data suggest that the frequency of finding a vacuum with large s has 
an approximately exponential falloff which we can describe by 

re(s) Ri noexp(- 7 s/smedian) • (5.19) 

Similar results are found for the barrier heights and distances of the lowest saddle points to the 


false vacuum. The values of 7 for various N are shown in Tables, 5.3, 5.4 and 5.5 Now using the 
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Figure 5.12; The dependence of the median value of s on the range of parameter 04 . The data 
shown are for N = 2. 


equation for medians from Eq. (5.16) in Eq. (5.19) 


n{s) « no exp ( — 


7 


-Ar«te 


(5.20) 


Ctension 

We can use this result to find an approximation for the distribution of tunneling exponents using 


Eq. (|5.14|). This suggests that the fraction of vacua with a tunneling exponent greater than some 

f{B) ~ exp , 


value B is roughly 


where 


/? = 10 


-3 


7 A 


(5.21) 


(5.22) 


dtension 

Erom the numerical analysis of this section, we see that 7 /Ctension is of order unity. We extracted 
an overall factor A from the shape of the landscape potential. There is no reason for it to be a 
small coupling constant. It could be well of order unity, which makes /3 ~ 10“^. Since atension ~ 


3, Eq. (5.21) shows a huge suppression in the number of long-lived vacua in high dimensional 
landscapes. 

We also show that the barrier heights that protect the false vacuum got very low, with a very 
similar relation 


/(14arrier) ~ exp - 


7 


/\/’f^barrier TZ 

^ barrier 


(5.23) 


C*barrier 

Therefore the fraction of barrier heights which are high enough to support stable vacua drop with 
an exponential of some power of N. This huge suppression means that only an extremely small 
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portion of barrier heights are high. Until now we have neglected gravity. But in the presence of 
gravity and a de Sitter space with a relatively flat barrier, the Hawking-Moss bounce provides an 
alternative mode of decay. The corresponding decay exponent is 

3 Ub arrier -^fv 


Bum = 


8G2 Ub arrier v-f. 


(5.24) 


arrier ^ iv 


where Vbarrier is the potential on top of the barrier and Gn is the Newton’s constant. Our numerical 
solutions suggest that the barrier heights relative to the false vacuum, (i.e., Ubarrier — 14 v) drop like 
an exponential of a power and therefore this mode of decay also gets enhanced for large number of 
fields and adds to the instability of the false vacuum. 


N 

Cubic 

Quartic 

SUSY 

1 

0.58 

0.46 

0.45 

2 

0.38 

0.39 

0.33 

3 

0.40 

0.35 

0.33 

4 

0.34 

0.33 

0.32 

5 

0.37 

0.35 

0.32 

6 

0.37 

0.34 


7 

0.38 

0.34 


8 

0.38 

0.35 


9 

0.38 

0.37 


10 

0.35 

0.34 



Table 5.3: Best fit exponents for the s distributions of results for cubic and quartic non- 
supersymmetric and quartic supersymmetric potentials with various N. 


5.3 Possible implications for string landscape 

The results we showed so far are general features of quantum field theories of scalar fields. However 
our main motivation in this work is understanding the implications of large dimensions for the 
tunneling rates of metastable vacua of the string landscape. It this section, we contemplate the 
applicability of our assumptions and results for the string landscape. 
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N 

Cubic 

Quartic 

SUSY 

1 

0.49 

0.75 

0.78 

2 

0.32 

0.69 

0.64 

3 

0.34 

0.65 

0.50 

4 

0.31 

0.60 

0.55 

5 

0.29 

0.60 

0.54 

6 

0.32 

0.59 


7 

0.32 

0.58 


8 

0.35 

0.59 


9 

0.34 

0.61 


10 

0.32 

0.62 



Table 5.4: Best fit exponents for the distributions of the lowest barrier heights of results for cubic 
and quartic non-supersymmetric and quartic supersymmetric potentials with various N. 


N 

Cubic 

Quartic 

SUSY 

1 

0.85 

1.42 

1.93 

2 

1.28 

1.52 

1.89 

3 

1.52 

1.98 

3.17 

4 

1.72 

2.04 

2.2 

5 

1.78 

2.34 

2.5 

6 

1.79 

1.95 


7 

1.53 

1.75 


8 

1.54 

1.68 


9 

1.52 

1.67 


10 

1.49 

1.54 



Table 5.5: Best fit exponents for the distributions of the distances of the lowest saddle points to 
the false vacuum for results for cubic and quartic non-supersymmetric and quartic supersymmetric 
potentials with various N. 
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Figure 5.13: The distribution of values of s in the ensembles for quartic (top) and cubic (bottom) 
random potentials for various values of N. The vertical axis represents the natural logarithm of 
the number of values in each bin. For the sake of clarity, the data for different values of N have 
been offset by constants, so only the slopes are meaningful. Purple diamonds correspond to = 2, 
green down-pointing triangles, = 4 ; red up-pointing triangles, N = 6 ] blue squares, N = 8 and 
black circles, N = 10. 
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Figure 5.14: The distribution of values of lowest barrier heights in the ensembles for quartic (top) 
and cubic (bottom) random potentials for various values of N. The vertical axis represents the 
natural logarithm of the number of values in each bin. For the sake of clarity, the data for different 
values of N have been offset by constants, so only the slopes are meaningful. Purple diamonds 
correspond to N = 2, green down-pointing triangles, = 4 ; red up-pointing triangles, N = 6 ] 
blue squares, N = 8 and black circles, = 10. 
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Figure 5.15: The distribution of values of the distance of the lowest saddle point to the false vacuum 
in the ensembles for quartic (top) and cubic (bottom) random potentials for various values of N. 
The vertical axis represents the natural logarithm of the number of values in each bin. For the sake 
of clarity, the data for different values of N have been offset by constants, so only the slopes are 
meaningful. Purple diamonds correspond to N = 2, green down-pointing triangles, = 4 ; red 
up-pointing triangles, N = 6 ] blue squares, N = 8 and black circles, N = 10. 
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Figure 5.16; The distribution of values of s (top) and lowest saddle point heights (bottom) in 
the ensembles for supersymmetric random potentials for various values of N. The vertical axis 
represents the natural logarithm of the number of values in each bin. For the sake of clarity, the 
data for different values of N have been offset by constants, so only the slopes are meaningful. 
Purple diamond correspond to = 1, green down-pointing triangles, N = 2 ■, red up-pointing 
triangles, N = 3 ] blue squares, N = A and black circles, N = 5. 
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Although the large number of possible vacua in the landscape of moduli of the underlying Calabi- 
Yau manifolds makes the string landscape a very intriguing place to start anthropic arguments for 
almost everything, it may contain the seeds of its demise within by destabilizing such an enormous 
number of vacua. Not only does the huge number of minima in high dimensional string theory 
landscapes (which can be the endpoints of tunneling events|^ offer the possibility of a significant 
enhancement of quantum tunneling processes, but we also have seen that the huge growth of 
the number of saddle points makes the lowest of the saddle points exponentially low and enhances 
thermal tunnelings. The main question of this section is to what extent the features of the landscape 
of random potentials that we have used in the previous sections provide an accurate insight into 
the properties of the string landscape. 

First we used a straight line that connects the false and true vacuum (to be more accurate, 
twice the contribution from a straight line that connects the false vacuum to a saddle points on 
the surrounding barrier) as an upper bound for the tension. But this very simple approximation 
may not be correct in the string landscape or in other high-dimensional spaces. The correct tun¬ 
neling trajectories may be very tricky and complicated. Explicit examples in the landscape are the 
conifunneling trajectories between monodromy-related flux vacua found in |79) (see also [801181| j. 
An additional, and potentially pivotal, complication is that the different flux vacua that are not 
monodromy related are generally minima of distinct potentials. Physically, such transitions invoke 
features not captured by our local field theoretic model, including for example the nucleation of 
branes to absorb changes in flux [371 Eg. These effects might significantly affect the tunneling 
action, and possibly mitigate the field theory instabilities we have identified. 

Second, we have assumed that as the stabilizing contributions to a given model are varied, the 
effective potentials around local minima will have expansions that are well modeled by random 


polynomials. In Eq. (5.2), if we keep to order D, there will be roughly random coefficients. 
But are there that many random free parameters to choose in a string landscape? In the landscape, 
the “random” parameters of the potential are the fluxes. There are far fewer fluxes in this system 
than the number of random coefficients we chose in our model. However, the associated minima 


^The density of vacua near the conifold point, Pconifoid, in a one-dimensional moduli space is described by 1/C{C + 
log r)^], where r is the distance from the conifold point [75]; applying this result near a generic point along the conifold 
locus an n-dimensional moduli space, shows he rapid growth in the number of vacua, f d"r Pconifoid with n. 
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of V will occur at different locations in the moduli space of a Calabi-Yau. As the location of 
these minima changes in the moduli space, it may help to pseudo-randomize the potential further. 
It seems reasonable to us that this will result in local expansions well modeled by the random 


potentials invoked in Sec. 5.1 but we do not have a firm argument. One important issue that we 
need to point out is that the decrease in barrier heights and tensions is merely a result of the fact 
that there are many positive and negative terms in the potential and that the number of coefficients 
grows exponentially with the number of fields. As long as these numbers are well scattered over 
positive and negative numbers, we expect the same results to hold. We probably do not need a 
perfect randomness as long as we have scattered the coefficients enough over positive and negative 
numbers. As mentioned, because the location of minimum can be anywhere in the landscape, that 
helps spreading these numbers, even if they are not random or freely chosen. 

Third, we used the results achieved from theories with canonically normalized kinetic terms 
to calculate the tunneling exponents. It is well known, however, that string vacua are densest in 
the vicinity of the conifold locus, where the classical moduli space metric suffers from a curvature 
singularity. In particular, near a generic point on the conifold locus we can choose local coordinates 
{Z^, Z'^,..., Z^) on the moduli space such that Z^ = 0 labels the conifold. Near Z^ = 0, the 
moduli space metric G behaves as: 


~ln(|ZY)- 


(5.25) 


The local form of the action then takes the form 

J ^/^[g^''Gij^^Z^^,Z^ -V{Z)], (5.26) 

where g is the space-time metric. In this expression the coordinates Z are the moduli space 
representation of the scalar fields (j) and V is their flux potential. At any non-singular point we can, 
of course, use a local change of field variables to absorb G into the Z, yielding a canonical kinetic 
term. But as we approach a conifold point, this change of variables corresponds to reducing the 
barrier heights in V (assuming V is continuous and is being expanded about a local minimum) and 
thus increases tunneling rates. In the regions of moduli space that are most densely populated with 
string vacua, we therefore expect the non-canonical kinetic terms to augment the destabilization 
we have found. 
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Fourth, since the flux potential V is derived from the superpotential W, and it is W that directly 
incorporates flux values, a more accurate representation of the landscape arises from randomly 
varying coefficients in a local expansion of W, and then using the result to calculate the random 
potentials. In principle, the relationships between the coefficients in V, which reflect its origin in 
IF, could alter our findings. For this reason we have repeated the calculations for a supersymmetric 
potential which is described in Appendix]^ The results, which were very similar to the random 
potentials, were presented in Sec. |5.2[ 

The examples we worked out in this chapter illustrate the challenges of investigating random 
potentials for large numbers of fields. In all of our studies, supersymmetric or not, computational 
considerations have forced us to only probe a limited range of values of N, the number of fields, 
and work with truncated potentials and use various approximations for tunneling rates. We are 


assuming that the pattern we have found, as evidenced in Figs. 5.3 5.11 will continue to hold as 
these constraints are relaxed. 


5.4 A multiverse explanation of the cosmological constant? 

Our results have important and immediate implications for attempts to use multiverse scenarios 
as an explanation for the smallness of the cosmological constant. These argument are based on 
simple counting. Even if the natural scale of the underlying theory for the cosmological constant 
is Planckian, one might naturally find vacua with A as low as if there are many more than 

10^^® vacua. However, a low cosmological constant is not the only feature of the observed universe. 
Our universe is several billion years old and whatever vacuum it is living in must be very long- 
lived. Therefore what we really need for an anthropic and natural explanation of the cosmological 
constant is an abundance (much more than of truly metastable and long-lived vacua. 

The generally accepted counting is that for a case of N fields, the number of vacua A/Vac is 
of order where g is some number of order unity, let’s say 10. However our results suggest 
that the chance of finding a stable vacuum decreases like where a is a model dependent 

number (in the cases we studied it was roughly 2.7 for quartic random potentials and 3.2 for 
supersymmetric potentials). Therefore, for very large N the chance of finding even one stable 
vacuum is exponentially small. If the number of vacua is g^, the number of metastable ones is 
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Figure 5.17: Parameter ranges allowing multiverse explanations of the cosmological constant. With 
a = 2.66 and g = 10, a sufficient number of metastable vacua is only possible for parameters in the 
region to the left of the solid line. This region is extended to the dashed line if g = 100. 
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A/'vac/(-Bmin) ~ 


(5.27) 


where f{B) given by Eq. (5.21). The requirement that this number be greater than 10^^^ is 


with 


N-^N^ > 120 

Ins' V Inc? 


= /3Bmin = XBm-,r, ~ lO’^^A 


Ctension / 


(5.28) 


(5.29) 


N must be a large number to satisfy Eq. (5.28), at least more than 120 A g = 10. The new feature 


we have found is that N cannot be too large, since the enhancement in the tunneling rate at large 
N will make the number of available stable vacua smaller rather larger. Therefore there can only 
be a window in N values, larger than let’s say 120, but small enough that the tunnelings leave us 
enough stable vacua to explain the cosmological constant. In fact there may be no such a window 


which satisfy this inequality, as illustrated in Eig. 5.17 where we used a = 2.66 from our results for 


non-supersymmetric potentials in Sec 5.2 The allowed region for b and N is on the left of the solid 
line for g = 10 and dashed line for g = 100. As we can see, there is no value of N which satisfy this 
criteria if 6 > 1.4 x 10“^ for = 10 and the range of allowed N is small until b gets much smaller 
than this threshold. Changing the value of g does not help much. For example for g = 100, the 
allowed values for b only increases by a factor of 4. We need to increase g exponentially to change 
the range of b significantly. 


In Eq. (5.28) a is model dependent. The effect of varying it is shown in Fig. 5.18 where we 


used a = 3.16 from our supersymmetric potentials in Sec. 5.2 Although the curve looks similar 
to the one of the previous case, the allowed values for b have shrunk by more than one order of 
magnitude. 


5.5 Summary 

In this chapter we studied the effect of increasing the number of fields on the stability of vacua in 
multi-field quantum theories. Our motivation arises from the landscape of string theory and whether 
or not it can furnish a natural framework to resolve the cosmological constant problem. We used 
random polynomials of quartic and cubic order with N fields {N < 10) as an approximation for the 
potential near its minimum. To make the computations tractable, we used some approximations 
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explained in Sec 5.1 and kept the number of fields relatively small. Even with these constraints, 
our data provide evidence that there is an exponentially large enhancement in the transition rates 
due to quantum and thermal fluctuations and that the rates increase so rapidly as a function of 
N that all but an exponentially small fraction of the generic vacua become unstable. Therefore 
the range of parameters that give a sufficiently large number of stable vacua to provide enough 
diversity to explain the cosmological constant problem is severely restricted. 

The assumptions we made were quite generic and potentially relevant to any model invoking 
anthropic explanation for the cosmological constant problem in which the required diversity of 
vacua is a consequence of involving a large number of fields. However, since our results are based 
on models of field theory, the impact of these considerations on the landscape of string theory and 
its ability to provide a natural explanation of cosmological constant problem needs further study. 
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Chapter 6 

Higher dimensional Einstein-Mcixwell 
landscapes 


Understanding the statistics, phenomenology, stability analysis and other features of the landscape 
of string theory is one of the most important and urgent problems in theoretical physics. Unfor¬ 
tunately, because of the complexity of this landscape, many of these features are not understood 
very well. As a way of understanding these features, different toy models have gotten attention and 
have helped us to uncover some of these features. Among these, the landscape of Einstein-Maxwell 
theory (for brevity we call it EM after this) is a very popular and interesting one. 

This is a model of electromagnetism in curved spacetime with a minimal coupling between 
electromagnetism and gravity. The magnetic flux can help to stabilize the compactification of some 
of the dimensions. For example, if we start from a six-dimensional spacetime and wrap magnetic 
fluxes around two of these dimensions, we end up with four non-compact dimensions. The remaining 
two dimensions have the geometry of a two-sphere. This compactification was first introduced in 
|83l [HI] and later on many aspects of this theory were understood. 

This model that we will describe later provides an excellent lab to test whether the field the¬ 
oretical approach we used in Chapter is applicable to potentials obtained from compactification 
because 

1. The compact extra dimensions are stabilized by fluxes wrapping around them. 


2. The kinetic terms are not canonical and even not diagonal. Nicely enough, there is a general 
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formula for making them diagonal and canonical. 

3. Although there are fluxes which corresponding to magnetic charges, we can go to a decom- 
pactified phase without nucleation of branes. 

4. There are inherently random numbers (the fluxes for different spheres). This resembles the 
situation in the string landscape. 

5. The effective potential obtained by compactification is exactly of cubic order. This makes the 
Taylor expansions we used in Chapter [^legitimate. 


The work we present in this Chapter is mainly based on a work in progress with Adam Brown 
and Alex Dahlen [85]. It will appear soon. We describe compactification of (4+2A^)-dimensional 
space in EM theory over C) (8> <8> ... <8> . In Sec. 6.1 we review some of the features of 

the 6-dimensional EM theory in which two of the dimensions are compactified over a two-sphere. 


This corresponds to = 1. In Sec. 6.2 we describe the generalization to higher dimensions and 
calculate the effective potential. We show some of the main features of the N = 2 case in Sec. 


6.4 and present some of the results and conjectures for higher dimensions in Sec. |6.5| Einally we 


summarize and discuss the future directions we will take in Sec. 16.6 


6.1 Review of the six-dimensional model 


In this section, we mainly follow the explanation and conventions in |86j . The action of six¬ 
dimensional EM is 

S = J ^ + A^ , (6.1) 

where ’^(e) and A denote the six-dimensional scalar curvature and cosmological constant. We 
use capital Latin letters to denote six-dimensional indices and Greek letters for four-dimensional 
indices. The equations of motion arising from this action are 


Rmn — -jgMNH 


1 

% 


Tmn, 




Qm 



0 


( 6 . 2 ) 
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where the energy-momentum tensor is 

Tmn = g^^FMpFNQ - -^Qmn - 9mnF . (6.3) 

We can compactify two of the dimensions over using a metric of the form 

ds^ = QMNdx^dx^ = Qf.udx^^dx'^ + ^ 4 ^ 

Here we defined the four-dimensional Planck mass as 


Ml = 


(6.5) 


The only ansatz for the electromagnetic equations consistent with the symmetries of the extra 
dimension and independent of the four-dimensional coordinates [83] is 

Tl 

=- —(cos6» ± 1) , (6.6) 

2g 

which is a monopole-like configuration in the extra-dimensional spheres. Here n is an integer and 
the signs correspond to the two patches necessary to describe a monopole configuration. Assuming 
that the two representations for the electromagnetic field are related by a single-valued gauge 
transformation quantizes n. The only nonzero components of Fmn are 


n 


Fg^ = -F^g = — sm ( 


(6.7) 


Using this ansatz for the electromagnetic sector, and the metric in Eq. (6.4), we can integrate out 


the two extra dimensions and get the four-dimensional effective action from Eq. (6.1) 


5 = / d'^xy/F^ -k + 17(V’)^ , 


( 6 . 8 ) 


where 

V(ih) - ^ ( r-^LIMp _ -2i,IMp , £3^^ -k>/Mp 

y2gml ^ Ml 

Eor any nonzero value of flux, this potential gets large and positive for large negative (small 
extra dimensions) and tends to zero for large positive V' (decompactifid phase). If there is no flux, 
the potential is unbounded from below and the extra-dimensional spheres collapse. The possible 
behaviors of the potential for nonzero fluxes is shown in Eig|6.1[ 
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V(iA) 



Figure 6.1: Different behaviors of the potential as we change the number of fluxes. For very large 
n, there is no stationary point (dotted purple line). By lowering n, the potential gets a single 
stationary point (dashed red line) and for smaller n (black line), it gets one maximum at a larger 
value of V' and one minimum at a smaller 'i/i. In the last case there is a stable compactified phase 
separated by a barrier from the decompactified phase. We are only interested in the last case. 


We are only interested in the case where there is a (classically) stable minimum. In this case, 
we can choose R in such a way that one of the stationary points appears at V' = 0. (Our convention 
is different from the one used in [86].) This choice of R is merely a shift in the values of V’ or, 
equivalently, a rescaling of R. Also, we define a dimensionless quantity h as = vrn^/ {2g^Mp) 
(Please notice that the electric charge in six-dimensional spacetime is no longer dimensionless). In 
terms of the rescaled R and n, the potential takes a very simple form. 


R^Vi^p) 

Ml 


= (2 - . 


( 6 . 10 ) 


Depending on n, the minima can be de Sitter, anti-de Sitter or Minkowski. For 0 < n < l/-v/2, 
this minimum is a de Sitter space. If n = l/\/2, the minimum is Minkowski. For larger n, the 
minimum is anti-de Sitter. But there are two possibilities. If l/\/2 < n < 1^2/3 there is a de 
Sitter minimum which is separated by a barrier with positive energy from the decompactification. 
For larger n’s there is no stationary point besides the anti-de Sitter minimum and this minimum is 
completely stable. Choosing a very large h creates a minimum with arbitrarily negative potential. 


These possibilities are shown in Fig. 6.2 This potential can tunnel to different values of h by 


creating a charged brane or by decompactifying. These instantons are explained in details in [86j . 
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R^WiiDIMp^ 



tpIMp 


Figure 6.2: Different possibilities for the minimum of the potential introduced in Eq. (6.10). As we 
change n from very small values, the minimum changes from being de Sitter (black dotted-dashed 
line), to Minkowski (red dashed), to an AdS with a positive maximum (blue line) to an AdS space 
which is never positive(dotted purple). These graphs correspond to h = 0.467, l/\/2,0.776 and 
0.826. 


6.2 Compactification and effective potential 


6.2.1 Calculation of effective potential 


In this section, we describe a compactification of (4 -|- 2V)-dimensional Einstein-Maxwell (EM) 
theory over (8) <8> ... <8> • As explained in the beginning of this chapter this model has 

^ -v 

N times 

very interesting features. Let’s start from the action of the EM theory in (4 -|- 2Al)-dimensional 
spacetime. As in the previous section, all the tilded quantities and also the capital Latin indices 
refer to the {N + 2)-dimensional quantities and the Greek indices denote the four-dimensional 
non-compact spacetime. The action is 


<S'em = J d - 


M, 


^ 1 

(4+^^(4+2n) ^ p.AB ^ ^ 

2 ' 4 


( 6 . 11 ) 


We compactify the spacetime over the product of N two-spheres and and write the metric as 


ds'^ = GMN{xK)dx^dx^ = e-^^^^^^^'>g^^dx^^dx’' + R^dQl^-^ . (6.12) 

Here is the length element of the i’th two-sphere. B and A are constant numbers and the 

■01’s only depend on the four-dimensional spacetime coordinates, i? is a length scale associated with 
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the size of the extra-dimensional spheres that keeps the dimensions correct. When we integrate 
out the extra dimensions, the action is in the Einstein frame only when A = B. Here we set them 
equal, and by a rescaling of the fields we set both to be equal to one. We are looking for solutions 
for the electromagnetic sector which are independent of the underlying four-dimensional spacetime. 
We want to choose these solutions in such a way that they stabilize the radius of the spheres. An 
ansatz which is compatible with the symmetries of the extra dimensions and independent of the 
four-dimensional coordinates is 


= ^sm0i , 

The other components are zero. The Lagrangian density of the electromagnetic tensor is 
FcdF^^ = 2 = 2 fI^^ = 2 ' 


(6.13) 


(6.14) 


Now we need to express the scalar curvature of the (4 -|- 2A^)-dimensional metric in Eq. (6.12) in 


terms of the associated scalar curvature of the four-dimensional space time. We have undertaken 
this calculation for the warped product of an arbitrary number of spaces in Appendix [C} The result 
for this special case is 


j^{A+2N) _ g'l' 


izA) + - -VT • VT - - ^ VV’i • Vt/; 


N 


i=l 




N 

E' 

1=1 


-Ipi 


(6.15) 


where T = covariant derivatives are calculated using the the connections 

obtained from the four-dimensional metric {g^u), so the inner product 

N \ / N 


VT-VT = V^( 


(6.16) 


^ 2 = 1 


We can readily integrate the action in Eq. (6.11) over the extra dimensions. We set the (4 -|- 2N)- 


dimensional Planck mass M(4+2iV) = 1- The effective actions of the different sectors are 


S'gr = - f ]^{A-kR?)^ 


^(4) + ^ |VT|2 - J ^+ I 2 E 




2=1 


2=1 


^It was argued in [8^ based on |84| that this is the only compatible answer for the case of the six-dimensional EM 
theory. We do not know whether or not this is the unique solution in the higher dimensional case. 
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Se-u = ld‘x{4wRYV^j22{i§2) 

S\ = [ d^x^y—g^AirR^)^Ae~^ . 




(6.17) 


Let’s define M = We can drop the second term in the first equation, because it is a 

total derivative. Putting everything together we get 

N , , N 


S = - f Sx-^gM^ !»<“> - 1 IV^r^ - ^ 1 |V*P + 4 E 

L i=l i=l 




-2 E (1^)4-*-’*-A.- 

2=1 ^ ^ 


(6.18) 


We have used the freedom to set the (4 + 2AA)-dimensional Planck mass to one, but still we can set 
RAK = 1. This is merely a shift in the ijjiS and is not setting a mass scale. This gives 

N 




- \ |V4/|2 - ^ i 


2 = 1 




(6.19) 


where 


Peff = f n? ) + e 


-'I' 


( 6 . 20 ) 


Here n* = If we define Oj = e the effective potential becomes a polynomial of order 

three. This verifies one of the claims we made in the beginning of this chapter . 


6.2.2 Bringing the kinetic terms into canonical form 


The kinetic terms for the scalar fields in Eq. (6.20) are not canonical. To diagonalize them, we 


need a field redefinition. There are many different ways to do this.We present two of them here. 
The first one is a symmetric field redefinition which treats all the fields on the same footing. The 
second one is not symmetric, but it is more useful for explicit calculations. 

Let’s start from the symmetric redefinition. Here rii = The kinetic terms of the scalar 

fields are not diagonal. To diagonalize them, we need a field redefinition. 

V’i = \/24i + , 

3 


( 6 . 21 ) 
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where 


A = -{N + l)±VWTl . 


To fix the notation, let’s choose the plus sign in Eq. (6.22). The inverse transformation is 

3 

In this notation, the kinetic terms are canonically normalized. 

|VT|V = . 

The effective potential in terms of the new fields is 


( 6 . 22 ) 


(6.23) 


(6.24) 


(6.25) 


where <I> = ^ (/>j . Although this field redefinition is very useful for the case that all the n^’s are 
the same, for other cases it is more convenient to use a field redefinitions which does not treat all 
the fields the same. 

V’l = + <^2 4- , 

V’2 = {N - 1)^2 , 

V’j = H-4 ^j-1 -{N - j + l)^j 


IpN = 4>1 + 4>2-\ -+ 4>N-1 — (t>N ■ 


(6.26) 


We can write this in a slicker way, 


= X] '^4 ~ ~ i>2 . 

j=i 


(6.27) 


The kinetic terms become 

|VTp + I= Y.^N + 1 - i){N + 2-i) 


N 


(6.28) 


i=l 


The inverse of this transformation is 


4>1 = + 'lp2 + ■ ■ ■ + i>n) , 
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4>2 = 


{N -1)N 


[•01 - (N - 1)02 + 03 + • • • + 0Af] , 


{N + l-j){N + 2-j) 


rr [01 - (IV + 1 - j)'lpj + 0J+1 + . . . + 0Ar] , 


0n = ^(01 -0n). 

Again we can write this in a closed form 

1 


{N + l-j)iN + 2-j) 


N 


tpl-{N +1- j)Tpj + ^ 0i 

i=j+l 


J >2 . 


(6.29) 


(6.30) 


This is diagonal. However, it is not canonically normalized and needs a trivial rescaling of the 
field.In order to go to a canonically normalized frame, we define 


0i — 


(iV + l-i)(A^ + 2 -i) 


The effective potential in terms of these fields is 

Here we present the field redefinition and effective potential for small values N. 


(6.31) 


(6.32) 


N=2 two-spheres 

The new fields are given by 


X 


1 


3 
2 

1 

V2 


3 

2 

1 

'V2 


01 

02 


and the effective potential Eq. (6.20) is 


Ves(X, 0 ) = 


N=3 two-spheres 

The new fields are given by 

f x\ 

01 

V 02 y 


1 

71 


( IX 

Vs \/3 Vs 
11-2 
V6 V6 V6 

VA 0 


v4 / 


01 ^ 
02 

V 03 y 


(6.33) 


(6.34) 


(6.35) 
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and the effective potential is 


VesiXA) = e 


-5xVfi 

-xiVs 




+ eA + 


1 

71 


( ^ ^ ^ 

V’2 

V’3 


2 

-3 


2 2 

1 1 

x/l2 x/l2 x/l2 ^/I2 
1 


0 ^ ^ 

^ s/6 V6 


V 7S 0 


0 


73 / V / 


(6.36) 


(6.37) 


N=4 two-spheres 

The new fields are given by 

/x\ 

(t>2 

\h J 

and the effective potential is 

Fefr(X, </>!, </.2, h) = e-37^^ I [(n?e-2^i + 

_g-^/72X I g<i>i^ ^-4>2/V^ g2<^2/73‘ g-03/76 g3</.3/76| g-2^/275X ^ 

This explicit form of potential is used in the calculations in the next few sections. 


6.3 Case of equal fluxes 


When all the fluxes are equal, it is easy to get some understanding of the potential in Eq. (6.25). 


Let’s assume that all the n(s take a common value n. We expect that at least some of the stationary 
points of the potential lie on the line where all the fields cjii take a common valne (f)- The effective 


potential Eq. (6.25) along this line is 

V^g: = iVn2e"V7lF(2+^)^ _ jsigs/2(i+N)4> . 


To make Eq. (6.39) look better, let’s define X = e 'I 24 >A+N/A) ^ 

Ves = X^ + Nn^X^+^ - NX^+^ 

The gradient vanishes at the location of critical points 

dV A + N 


dcpi 


A 


\/2X^ {n^X^(N + 2)-(l + N)X + l} . 


(6.39) 


(6.40) 


(6.41) 
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Equation Eq. (6.41) may admit at most two solutions, 

iV + 1 ± ^{N + iy - 4n‘^{N + 2) 


X = 


2n2(iV + 2) 


(6.42) 


If both roots are positive numbers, the larger one, which corresponds to smaller spheres, will be a 
minimum and the smaller one will be a maximum along the line </> = (/>. To understand whether 
these points are minima, maxima or saddle points, we need to calculate the other eigenvalues of the 
Hessian which correspond to fluctuations of potential in directions not parallel to <h. After some 
simplification, the Hessian matrix at the critical points is 


Hij = 2X 


N 


(4n^A^ - X){y^ + A(y + i) + 5, .) _ 2n^xV 


= 2X^[UAj + U2] , (6.43) 


where y = 1 + ^. Please notice that the terms which do not have indices on the right side are 
independent of i and j. The eigenvalues of this matrix are 2X^(Ui + NU 2 ) and an {N — l)-fold 
degenerate eigenvalue 2X^Ui. By looking at the behavior of the potential in Eq. (6.40) along 
the line cjii = 0, it is clear that the non-degenerate eigenvalue is negative at the smaller root and 
positive at the larger root. The other eigenvalues are all equal to Lf\ = 4n^A^ — X. The value of 
t/i at the larger root (the one with positive sign in Eq. ( 6.42[ )) is always positive. Therefore if this 
point exists, it is a minimum. 


N + 2^(1 + NY - 4n2(2 + N) ~ 

=-2VIV-^ ’ 

The value Ui at the smaller root (corresponding to bigger spheres) is 

N - 2v'(l + NY - 4n2(2 + N) ~ 

= -2TIV-^ ■ 

The sign of Ui can be easily expressed in terms of ric defined as 

ric = -\/2 -|- 3N . 


(6.44) 


(6.45) 


(6.46) 


If n < ric, Ui is negative and therefore this point corresponds to a maximum. For ric < n, Ui 
is positive and therefore this point is a saddle point with one negative mode and N — 1 positive 
modes. The vacua can be de Sitter, Minkowski and AdS. The Minkowski vacua happen at 


Vn 

nm = — 


(6.47) 


For n < riM, the vacuum will be an AdS space and for um < n, it is a de Sitter space. 
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There is no minimum for 


n > riniax 


A^ + 1 

2VN+ 2 ■ 


(6.48) 


We know that when the top of the barrier gets too flat, the CDL instantons cease to exist. The 
criteria for the existence of CDL bounces is |24l - l26] 


tj2 

■^top 


(6.49) 


When this criteria is met, the dominating bounce is the CDL one. However if the CDL does not 
exist, the Hawking-Moss bounce dominates. Therefore the tunneling mode has a transition from 
CDL to HM at V"{4)top) = 4iLtop ■ This happens at 


(1 +A^)\/6TiV 

ncDLHM = 2(4 + iV) 


(6.50) 


Interestingly enough, Uc, um and ncDLHM all coincide for N = 2. But this is only a coincidence. 
For all other N's this does not happen. These lines are shown in Fig |6.3[ 

Unfortunately, the case where the fluxes are not equal is not as easy as this. In the next section, 


we derive many features of the potentials for N = 2 and in Sec, 6.5 we derive some general results 
and also state some conjectures about the statistics and the Hubble parameter of the landscape of 
these vacua. 


6.4 Case N = 2 

In this section we study the features of the N = 2 case, which is the compactification of eight¬ 
dimensional Einstein-Maxwell theory on two two-spheres. 


6.4.1 Phase diagram of the N = 2 case 

Let’s assume that there are ni fluxes on the first two-sphere and n 2 on the second one. The effective 


potential for this case in canonically normalize fields is given in Eq. (6.34) . We want to know the 
number of minima, maxima and saddle points for different numbers of fluxes. A summary is shown 


in Fig. 6.4 The region above the top red line does not have any minimum. This resembles the 


= 1 case where increasing the flux causes the minimum to vanish. Along this line the minimum 
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n 



N 


Figure 6.3: The possible situations for the critical points with all the spheres of the same size for 


the effective potential of Eq. (6.25). There is no minimum (vacuum) above the black dot-dashed 
line. The blue solid line specifies the Minkowski vacuum, below it are the AdS and above the dS 
vacua. The transition between the CDL and HM bounces happens at the red dotted line. Above 
the purple dashed line, there is a saddle point and a minimum and under it a maximum and a 
minimum. 


and the saddle point annihilate each other. To find the equation that the red line satishes, let’s 


rewrite the potential in Eq. (6.20) as 




(6.51) 


where 

V = 1 + -p - 2n?e"2’^i - 2nie"^2 

At any stationary point of the potential 

^ dV dVV , 1 

~Wi~~ 


This leads to 


V = . 


(6.52) 


(6.53) 


(6.54) 


At the red line, the minimum merges with the saddle point. Therefore one of the eigenvalues of 
the Hessian should vanish. The Hessian matrix at the extremum is 

d‘^V _ y d f _yl I dV dv\ 

d-ipidil^j d'lpj \ V dt/ji J Vvdil^idiij dijji dif^j j 
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= V 


dV 


Vdipidil^j 


-1 = 


V 

V 


4n?e - e 


6,, - V 


}■ 


(6.55) 


Here we use a simple theorem about the determinant of matrices. For a matrix M defined by 
Mij = aiSij — b, the determinant of M vanishes if a~^ = b~^. Therefore the vanishing of one of 
the eigenvalue of the Hessian translates into 

4n^e~‘^Si — e~Si y ^ ^ 

i ^ 

This leads to 

y ^ = 1 . (6.57) 

4nJe-Si - 1 ^ ^ 

i * 

After we combine Eq. (6.54) and Eq. (|6.57|) and perform some manipulation, the equation for the 


red line in Eig. 6.4 is given by 


2358nfni + 58848nfn| + 65536n^n^ + 1125 (nf + nf) - 15420 {nfnl + nfn|) 
-2500 (nf + nf) + 22800 (nfnf + nfnf) - 67584 (nfnf + nfnf) = 0 . 


(6.58) 


The same equation also describes the lower part of the red line where a maximum and saddle point 
merge and disappear. On the left side of the vertical dashed blue line and also under the horizontal 


blue dashed line in Eig. 6.4, one of the saddle points hits the infinity and disappears. We give a 
proof for this in Sec )6.5[ where we prove a generalized version of this statement for arbitrary N. 

The vacua that we get from the minimum of these potentials can decay through decompactifi- 
cation. This decompactification can be carried either by CDL or HM bounces. But if the top of 
the barrier gets too flat, the CDL disappears, and using a similar technique, we find that the line 
where the transition between CDL-dominated decay to HM-dominated decay occurs along a curve 
[the criteria for this transition is V”{(j)top) = 

1286nfnf + 44256nfnf + 65536nfnf + 637 {n\ + nf) - 10012 (nfnf + nfn|) - 1764 (nf + nf) 

+17808 (nfnf + nfnf) - 59392 (nfn| + nfnf) = 0 . (6.59) 

This is the red dashed curve in Eig j6.5[ The Minkowski vacua occur along the green solid curve in 
Eig|6.5| which is described by 


1 1 

2 2 
nf nf 


= 4 . 


(6.60) 
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This is a special case of the result we present for arbitrary N in Sec 6.5 



Figure 6.4: The phase diagram of the number of minima (++), maxima(-) and saddle points(H—) 


for the N=2 case. The two curves satisfy Eqs. (6.58) and (6.59). There is no stationary point above 
the red line. On the left of the vertical dashed blue line (and similarly under the horizontal dashed 
blue line), a maximum and saddle point annihilate. The derivation of the curves and also the 
annihilations on the left of these lines are derived in the main text. 


Now we can give a full description of the phase diagram presented in Fig. |6.5[ In region 1 there 
is no stable minimum. Region 2 has a single dS vacuum and the transitions are dominated by the 
HM bounces. Region 3, the tiny area between the green solid line and the dashed red line, has a 
single dS vacuum and the transitions are dominated by CDL bounces. Region 4 has a AdS vacuum 
and a saddle point and the CDL bounces are dominating. Region 5 has a AdS vacuum, two saddle 
points and one maximum and the CDL bounces dominate. Region 7, 8 and 9 have a AdS minimum 
and the potential approaches zero from below for large 'ip’s and therefore there is no tunneling in 
this region. They are only different in the number of saddle points and maxima. Region 10 has a 
AdS vacuum, one maximum and two saddle points and the HM bounces are dominating. 
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Figure 6.5: The phase diagram for the N = 2 case. The green line given by Eq. (6.60) shows the 
location of the Minkowski vacua, above it is the dS and below it the AdS minima. There is no 


minimum above the top dotted-dashed blue line (region 1) which is given by Eq. (6.58). The red 


dashed line, which is given by Eq. (6.59), shows the boundary where the CDL bounces vanish and 


the HM bounces dominate . The most interesting part of this landscape is the region between the 
top red line and the green line where all the dS minima live in (regions 2 and 3). As proven in the 
text, its area, and therefore the number of dS vacua, is hnite. 
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6.4.2 Statistics of vacua 


It is not possible in general to find an analytic expression for the extrema of the potential Eq. (6.20) 


or, equivalently, Eq. (6.52). But we can infer some useful information about them. The area of the 


regions to the left or under the green line is infinite which means that there are an inhnite number 
of AdS vacua. Also, by making the n’s small enough, these minima can have arbitrarily negative 
cosmological constants. The story for the more interesting dS vacua is different. First we show 
that there are a finite number of them and then we find approximately the location of the minima 
and their cosmological constants for large values of fluxes. The dS vacua are bound between the 


two lines defined by Eqs. (|6.58|) and (|6.60|). For large values of ni, the tail of the green line is 

(6.61) 


1 


1 


“ 2 + 4n? ’ 


and the tail of the green blue dotted-dashed line is 

1 1 

Tloh — ~ “t“ - n 

2 Unj 

The area between these two lines is 


(6.62) 


dni [n2g - n2b] 


(6.63) 


which is finite. From Eq. (6.61) and Eq. (6.62) it is apparent that the right parameterization for 
the large values of n 2 is 

1 2 

ni = - + 6 , 


K 

n2 = — , 
e 


(6.64) 


where e can be arbitrarily small. In terms of this parameterization, the minimum is located at 


1 / 1 Vl - 12A:2 \ 5 

= 2 -b - -96 T 


9 


K2 


K2 I 






3a:2 


Trnin — 


2(1 + Vl - 12kA ( 24^2 - 1 - Vl - 12kA 


27A'4 


(6.65) 
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It is easy to see that K < . Now we can express I4iin in terms of ni and n 2 

2 ^1 + + (6 — 12ni)n2^ ^1 + (12 — 24ni)n2 + -\/l + (6 — 12ni)n2^ 


Idnin — 


27n| 


( 6 . 66 ) 


The eigenvalues of the Hessian in a canonically redefined field at the minimum are 


Ai — 


A 2 — 


Vi - i2K‘^ (^1 + Vi - i2if:2y 


9/74 

16 f 1 + Vl - 12kA 


9it'2 




(6.67) 


Therefore the potential energy at the minimum is of order e^. One of the eigenvalues is of order of 
e^. However, the other one is of order e^. Therefore one of the mass scales of excitations around 
the minimum is much smaller than the energy scale of the vacuum. 


6.5 Higher N'’s 

It is very difficult to completely map out the landscapes for N > 3. However we can infer some 
important and general features of them. 


6.5.1 Criteria for Minkowski vacua 


We can rewrite the effective potential in Eq. (6.20) as 


Ves = e-^^M , 


where 

At a Minkowski vacuum, H = 0 and dVjd'^i = 0. This leads to 


V 

dV 

dA 


0 , 

= 0 . 


( 6 . 68 ) 


(6.69) 


(6.70) 
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Form the second equation we can find all the V’i’s and plugging it back into F = 0 we get the 
criteria for Minkowski vacua. 


^ 1 


^ =4 
1=1 ^ 


(6.71) 


The dS vacua correspond to < 4 and AdS vacua to > 4. Equation Eq. (6.60) is 


N 1 


a special case of this result. 


6.5.2 Bounds on valnes of n,’s 


As we have seen in the previous sections, a stationary point does not exist for every choice of fluxes. 
Here we infer some bounds on the value of fluxes which is necessary for existence of a stationary 
point. It is more convenient to write the potential in Eq. (6.20) in terms of pi = 


Ves = Pi... pn fipi, ■■■, Pn) 


(6.72) 


where 


/(pi, ...,pn) = l + Y^ {n^iPi - pi) . 


The critical points must satisfy 


/ + 1 “ Pj) + ‘^^‘iP'i -Pi = ^ 

^Pi j=i 


N 


(6.73) 


(6.74) 


for every value of i. Therefore, at the critical point, pi{2nfpi — 1) should be independent of i. By 


completing the squares we can rewrite Eq. (6.74) as 

E (ft - = i + 


i=i 


4n2' 


Snf 


An?. 

3 ^ 


Therefore 


This leads to 


Mi 




1 


An^. 

3 


^ 1 8iV 
^2 — 2N + 1 


3 


(6.75) 


(6.76) 


(6.77) 


Because we cannot set the left-hand side of Eq. (6.75) to zero for all f’s, this is not the strictest 


bound we can set on n^’s. For large values of N, even this non-strict bound converges to the bound 
we found for the dS to AdS transition which means that the window in which we get a de Sitter 
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space gets smaller and smaller by going to higher dimensions. However we will see soon that the 
number of vacua diverges for N > 2. 

6.5.3 Location of the critical points 


We can rewrite Eq. (6.75) as 


1 


1 


4^2 ) sn? 

where A is defined as 


E 4^2 1 E 


n] pj 


2nl 


8nf 


1 




2„| 


Solving Pi in terms of A leads to 


Pi = 


4n? 


1 ± 


(6.78) 


(6.79) 


(6.80) 


If A is positive, only the positive sign applies to this problem. Now plugging Eq. (6.80) into 


Eq. (6.79), we get an equation for A, 

■^ = ^ ± \/l + . 


N + 2 


(6.81) 


Eor positive A, only positive signs apply here. However, for negative H’s, the situation is not that 
clear and both positive and negative signs are applicable. The possibility of plus or minus in each 
term of the sum makes this equation less useful, but we will still be able to use this equation later 
when we expand it in terms of small piS. 


6.5.4 Case of small fluxes 

We saw in Fig |6.4| that for small values of the n^’s, there always exists a stationary point which is a 
maximum(the diamond shape under the red line). In this section we show that this is indeed true 
for general N. 


Starting from Eq. (6.74) and expanding to 2nd order in Uj, we get 

1 


Pi = 


iV + 1 (iV + 1) 




2nj 


(iV + l)2 ’ 


(6.82) 
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and, using the definition of / in Eq. (6.73) 


f{Pli • ■ ■ ) /5Af)|criticalPoint — N 1 \ ^ y 


nj 


(IV+1)2 


+ o(n: 


(6.83) 


It is easy to show that 

= Pi... PN 2n‘fdij - ^ =-pi...pNiN + l)[dij + l + 0 {nj)] . 

upiUpj critical point L 

(6.84) 

The eigenvalues of this matrix to zeroth order in m are proportional to {N +1) and (iV + 1)2 which 
confirms that this extremum is indeed a maximum. 


6.5.5 Vanishing critical points because of hitting infinity 

In this section we find the criteria for changing the number of critical points due to the approach 
of some of the fields V’t to infinity for general N and find the location of these saddle points. It will 
be a straightforward calculation to determine whether they are minima, maxima or saddle points. 
If the r’th radion approaches infinity, then pi = will approach zero. Let’s separate the pi's 
that remain finite by labeling them as pi,..., pK and the ones that approach zero as px+i ■ ■ ■ Pn- 


From Eq. (6.80) we see that this behavior is possible if and only if A is an infinitesimally small 


negative number and we have to choose the positive root for i = 1... K and the negative root for 


i = K + 1,... N. Now expanding Eq. (6.81) to first order in A, we get 

K 


E 

i=l 


1 

4n? 


= 1 + (V + 1 - K)A = 1 - e , 


(6.85) 


where e is an infinitesimally small positive number. Therefore, for any subset ni... the number 


of stationary points of the potential changes at a hypersurface defined by Yld=i 
location of this critical point will be 




Pi = 



= 1. The 


( 6 . 86 ) 


The two blue dashed lines we found in Fig,6.4 (ni = i or n 2 = i) for N = 2 are special cases of 


this argument and we expect that a saddle point vanishes on the left (bottom) of these lines. Using 


the value of pi we found in Eq. (6.86) we can easily determine whether these points are minima, 
maxima or saddle points . 
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6.5.6 Divergence of the number of dS vacua for N > 3 

One of the most interesting features of this landscape is that not also is the possible number of AdS 
vacua infinite for N = 2, but also for larger N they accommodate an infinite number of dS vacua. 
In this section we prove this result for the case N = 3. It will be clear from this proof that the 
larger NA also have an infinite number of dS minima. Clearly, for any finite value for the fluxes 
ni, n 2 and ns there is a finite volume in parameter space. Therefore if we want to find an inhnite 


volume, these numbers should get very large. However, from Eq. (6.77), it is not possible for all of 
them to get large. Let’s assume that ni remains hnite and n 2 and ns get large. We are looking 


for dS minima. Using the argument in Sec,6.5.2 and Eqs. (6.77) and (6.71), ni should remain very 
close to 1/2. We use the following parameterization for ni and n 2 and ns: 


m = 2+6 


n2 = 


ns = 


K 

e 

L 


(6.87) 


where again e is a small positive number and K and L are positive numbers that should be chosen 
in such a way that ^ + ^ + ^<4. To the hrst non-zero approximation in e, the location of the 
minima is given by 


Pi — 2 + Ae^ , 

P2 = Be^ , 

P3 = Ce^ ■ 


In order to satisfy Eq. (6.74), A, B and C satisfy 


12+ A-B-C + B^K'^ + = 0 , 

A-2B-CP 3B‘^K^ + = 0 , 

A- B-2C + B'^K'^ + = 0 . 


( 6 . 88 ) 


(6.89) 


Unfortunately, it is not very easy to solve these equations analytically and analyze the minima 
as we could do in the case of two two-spheres. However, we can prove that the volume in the n 
space that creates dS vacua is infinite in contrast to the N = 2 case. In the case of two two-spheres, 
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we found a finite number of de Sitter vacua because the area between the graphs in Eqs. (6.61) 


and (6.62) was finite. However, if we rotate these two graphs around the rii axis, the integral in 


Eq. (6.63) diverges logarithmically. To make this argument more precise, let’s adopt a cylindrical 


coordinate system as follows: 

ni = ni, 


n2 = rcos9, = rsinO . 


(6.90) 


The two graphs that specify the boundary of the region containing the dS vacna are 


1 

= 2 


1 

ni = - 

2 


1 + 


1 + 


1 


1 


8 sin^ 9 cos^ 9 
G{9) 


(6.91) 


where the first one is a result of Eq. (6.71) and the second one comes from Eq. (6.89). Our numerical 
solutions show that 


F{9) = G{9) -^-— > 0. 

8 sin^ 9 cos^ 9 

Therefore the volume of phase space for the de Sitter vacua is 


oo r^/2 

rdrd9AH 


poo p 

Jr^ Jo 


poo 1 /*77/2 

( dr- I d9F{9) . 


(6.92) 


This integral is logarithmically divergent and therefore the number of dS vacua is infinite. A similar 
argument is correct for higher N and our conjecture is that this integral for different NA gets a 
modification of the form 


OO /*7r/2 


poo p 

Jr:, Jo 


^drdQN-2^-^ 


fOO 

dr 


^7r/2 

/ ‘‘0^(0) . 

Jo 


(6.93) 


This expression is more divergent for larger NA. Therefore in all of them we will have an infinite 
number of dS vacua and also an infinite number of AdS vacua with arbitrarily large negative 
cosmological constants. 


6.6 Summary and future directions 

In this chapter we studied compactifications of the (4 + 2A^)-dimensional Einstein-Maxwell theory 
over a product of a four-dimensional Lorentzian space and two-spheres. We calculated the four¬ 
dimensional effective potential. The landscape of vacua for = 2 is completely mapped out. We 
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showed that the number of dS vacua in this case is hnite and that there are an infinite number of 
AdS vacua. We studied the minima of this theory in the approximation that one of the fluxes gets 
very large. Although we could not completely study the landscape for higher N's, we still could 
derive some of its important features. Especially we showed that there are an infinite number of 
dS vacua for higher NA. 

There are many things left here that we will pursue in future works. One of them is the stability 
of these spheres for higher i modes. We also will check whether this model of compactiflcation 
becomes unstable for large N, as we could show for field theoretical models in Chapter 
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Chapter 7 

Conclusion 


The Universe can undergo large-scale phase transitions. When the hot and dense Universe got 
colder, there was a chance that it did not land on the true global minimum of the potential 
describing it. The way to move towards a more stable vacuum is a phase transition which is carried 
by nucleation of bubbles of the new phase. The subsequent expansion may complete the phase 
transition. In this thesis I explained aspects of this decay. 

In Chapter I laid the foundation for the subsequent chapters. I described the vacuum decay 
in a field theory. 

In Chapter]^ I explained the bubbles which have an 0(3) x 0(2) symmetry, which is different 
from the conventionally used 0(4) symmetry. I showed that these bubbles have higher actions and 
therefore are subdominant in the decay process. 

In Chapter I showed the results for the decay of a spatial vector field theory with different 
transverse and longitudinal speeds of sound. I showed that in some limits the flat walls get unstable 
and the bubbles develop kinks. 

In Chapter I studied the effects of large numbers of fields on the decay rate. I showed that 
the decay rates grow so quickly as to render the overwhelming majority of the vacua unstable, and 
therefore not good candidates for our observed Universe. 

In Chapter I presented a model which shares many of the features of the string landscape 
and makes a good lab for testing the assumptions made in Chapter We will investigate the 
consequences of large number of fields in this model in upcoming work. 
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Appendix A 

de Sitter space 


The d-dimensional de Sitter space dS^ can be thought as a hypersurface in a d-|-1-dimensional 
Minkowski space defined by m 


-X^ + Xi + ... + Xj=f 


(A.l) 


where I is the de Sitter radius. This space satisfies the Einstein equation with a cosmological 
constant related to i by 


A = 


{d-2){d-l) 

2£2 


Because de Sitter space has a horizon, it has a temperature [HHl |89] which is given by 


(A.2) 


(A.3) 


We now set i = 1. If we assume Xq is the time direction, then we can foliate the spacetime with 
constant time hypersurfaces. Each of these hypersurfaces has the geometry of a (d-l)-sphere of 
radius + Xq. This foliates the space with first contracting and then expanding spheres. This 
space has a SO{d, 1) symmetry group. A simple picture of this space is shown in Fig ]A.l| There 
are many different coordinates systems used for de Sitter space. Here we introduce the two which 
are needed in the thesis. The first one is the global coordinate system which describes the above 
embedding. Let’s devise a spherical coordinate system for the space-like directions 


UJ^ = cos 01 , 

= sin 01 cos 02 , 
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Figure A.l: A schematic view of the de Sitter space. Each of the points here is a (d-2)-sphere. 
The dashed line shows a (d-l)-sphere at the end of the contracting and beginning of the expanding 
phase. 


^ = sin 9i ... sin 9^-2 cos 9d-i , 

= sin 01 ... sin 9d-2 sin 0rf_i . (A.4) 

Here 0 < 9i < tt for 1 < i < d — 1 and 0 < 0rf_i < 27r. In terms of these variables, we can put a 
global coordinate system on the de Sitter space 

Xq = sinhr , 

Xi = iOi coshr . (A.5) 

The metric takes a very simple form 

+ cosh^ rdH^_i . (A.6) 

The Penrose diagram for de Sitter space is shown in Fig |A.2[ Not all the points are causally 
connected and for any point there is a horizon. The shaded area in this picture shows all the points 
that can send and receive signals from the north pole (the solid vertical line on the right). We can 
choose a coordinate system which only covers this shaded region 

Xq = y/l — sinh t , 
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r 



Figure A.2: Penrose diagram for the de Sitter space. The vertical solid lines are time-like points 
corresponding to the points on the north and south poles. As seen, it takes an infinite amount of 
time for a signal to travel between the poles and the north pole is only causally connected to the 
shaded region which is called the static patch. 

= ra;“ , a = l,...,d—l, 

Xd = Vl — cosht . (A.7) 

In this coordinate system the metric takes a very simple form 

ds^ = —(1 — r^)dt^ + ^^2 A • (A.8) 

This is a static metric and r = ±1 denotes the boundaries of the shaded area. Clearly the Euclidean 


version of the space given in Eq.(A.l) is a d-sphere. 
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Appendix B 

Examples of potentials with various 
orientation dependence of the tension. 


Here we present some analytic approaches and the numerical evaluation of the domain-wall tension 
given by the relaxation method [5lll56| . The first potential we study is a double well in one direction 
and a quadratic in the other direction. 


^(py) — + n 


(B.l) 


This is qualitatively similar to Eq. (4.7) when > 0. The differences are some 4th order terms 
involving 4>x, which is not very important when /3 > 0 stabilizes a trajectory near (px = 0. The two 
potentials can be roughly related by 


fx 

2 

2 


9 

m 

^ ■ 


(B.2) 


The two degenerate minima sit at (0, iy^). For the purely longitudinal (or transverse) wall 
oriented along the y (or x), we can solve the problem analytically and get the exact value of the 
tension. 
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Longitudinal wall (0 = 0): 


4>x{x,y) = 0 
cr = cr(0) = 


2V2fi^CL 


3A 


Transverse wall (0 = f): 


(l>x{x,y) = 0 

2^/2 yAcT 


a = cj( 7 r/ 2 ) = 


3A 


(B.3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B.8) 


For other orientations of the wall, we can evaluate the tension numerically. Before that, we can 
analyze two extreme cases. Using the method of rotating the potential as described in Sec 4.2, we 
have 


ct(0) = / dx + ye{(t>xAy) 


When jl 


such that 


f cA (P‘ 

= dx + -^4>y + y {4>x COS 9 + (j)y sin 0, c^y cos 6 - cj)x sin 0) 

oo, we effectively have a single field problem with 

7 ^x (f^y 

^ =-h = ’ 

cos 0 sm 0 

cos^ 6 + sin^ 0 -,2 


(B.9) 

(B.IO) 

(B.ll) 


a{9) = J dx ^ + y {(/>, 0) . 

Clearly, this gives us Eq. ( 4.19| ). 

The other extreme limit is /3 —)■ 0, in which the potential is flat in the (/>x direction. The two 
degenerate vacua approach two separated lines. Moving along these lines contributes nothing to 
the tension. As shown in Fig jB.l[ the path that minimizes the tension involves first moving along 
these lines to an appropriate angle cj), then connecting them through a straight line. The tension 


of this path is a function of both 0 and (j) through the orientation dependence in Eq. (4.19), and a 
simple projection of the length. 

C7(0) 1 


a{e,(j)) = 


Cl cos(0 — 4>) 


c\ cos^ (j) + sin^ (f) . 


(B.12) 
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Figure B.l: The cl and ct axes are the directions in which the field has purely longitudinal and 
transverse sound speeds. The two dots represent the two discrete vacua. In the limit /? —>• 0, the 
dashed lines through them are almost in the vacuum, too. The important portion of the domain- 
wall is the red (thick) path from one line to the other, which is free to pick the best orientation 
(p. 


Minimizing this with (p, we have 


(pm{d) = arccos 


cos 9 


Cy cos^ ^ sin^ 9 


(B.13) 


so the tension in this case should be 


a{9) = a[9,(pmi9)] . 


(B.14) 


We can apply the analysis in Sec, 4.2.2 and calculate the stability condition for the fiat longitu¬ 
dinal wall: 


1 (Pa 
a(0) d9‘^ 


= 1 - 


> -1 


(B.15) 


We can see that the wall becomes unstable as soon as cl > \/2ct. 

We next provide several plots with the numerical values on top of the three possible fits, 
Eqs. (4.18), (4.19) and (B.14). Fig|B.2| shows that the two extreme limits indeed fit very well 


with our analysis. Fig.(B.3) shows that with a more moderate choice of parameters, Eq. (4.18) is 


quite reliable, independent of the sound speed. 
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cr(6>) 


Ci = 15.8 Cj=5 A = 1 // =4 jS = 0.01 


(T{e) 


Cl = 15.8 ct=5 a = 1 // =4 jS = 500 




Figure B.2: The numerically calculated values of the tension for a double-well potential are shown 


in dots. The three analytical fits: Eq .(4.18) is the dashed line, Eq. (4.19) is the dot-dashed (blue) 


line, and Eq. (B.14) is the solid (red) line. We can see that for large (3, Eq. (4.19) is a good fit, and 


for small /3, Eq. (B.14) is a good fit 


Cl = ct=4 a = 1 fi=4 p = 5 

cr(6>) 



Cl = 10.2 cr =2 X= I ^=4 13 = 5 

cr(0) 



e 


Eigure B.3: The numerically calculated values of the tension for a double-well potential are shown 


in dots. Again the three analytical fits: Eq. (4.18) is the dashed line, Eq. (4.19) is the dot-dashed 


(blue) line, and Eq. (B.14) is the solid (red) line. The two figures use the same potential but 


different sound speed ratios. 
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Figure B.4: The more complicated potential introduced in Eq .(B.16), with q = 0.5 , ri = 


—2.5 , r 2 = 2.5 . We will use two values of S, 1.1 and 0.9, but that makes no visual difference. 


In the end, we provide a much more complicated potential as in Fig. (B.4). It has a general 


slope in the (j)x direction and two minima located at (0.001, ±2.498). So, trivially, the interpolation 
path will always involve both fields. 

2 ‘ 


V4>y) = ‘(1-5'exp 


-4{(l)x - sin( 


r2 - rr 


tanh^ 


hy - ri){(j)y - r 2 ) 


(B.16) 

We numerically evaluated the tension for various orientations and plotted it against the three 


analytical fits in Fig. (B.5). The overall shape can be quite different from any equation given in 


this paper. In particular, note that in the right portion of Fig. (B.5) the longitudinal domain-wall 
(actually an open set near 9 = 0) does not exist 0 


^This is for the same reason as described in [54]. The interpolation path breaks into two parts, connecting each 
vacuum individually with the —(px region. For vector helds, such runaway behavior also acquires an orientation 
dependence. 
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Figure B.5: From the potential given by Eq .(B.16), we again compare the numerical (j{0) with the 


three equations. In the left figure we have S' = 1.1. In the right figure we have S' = 1. For some 
orientations the domain-wall does not exist because the path runs away toward the —4>x direction. 
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Appendix C 

Curvature of warped product of 
metrics. 


In this appendix, we calculate the curvature tensor for a warped product of metrics. These results 
were used in Chapter to calculate the effective potential of the lower dimensional spacetimes 
after compactification. We do these calculations in two different parameterizations of the metric. 
Consider the warped product of manifolds AA x x .. which are respectively m, ni,n 2 ■ ■ ■ 
and Up dimensional. The coordinates of N4 are Xp and are the coordinates of 

C.l Parametrization 1 

The metric in the hrst parameterization is 

ds^ = GMNdx^dx^ = gf,uix)dxf"dx'' + (y(„)) , (C.l) 

where the index a runs from 1 to p. Latin indices from beginning of the alphabet specify which 
manifold Af is used. The nonzero components of the affine connection are 
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= U. - 2<5'' („v„, A + A 

(C.2) 

Ku = 

(C.3) 

p*a _ Y^a 

jUa jAa 

(C.4) 

pio — xia^J R 

^ Ma - ^ja 

(C.5) 


All covariant derivatives with Greek indices are with respect to the affine connection obtained form 
9)inu • The nonzero components of the Riemann tensor are 


^uajS ~ [a^ A + 2(5^ [Q,V^]AV,yA 

+ 25 !/[/ 3 Vq,]AV^A - 25^ (VA)^ , 


IT- = IT;. , - [fc (VBaf , 

Ja^af'a Ja^a^a I la\la ' ' 


Ja^a^a 

' (a) 

"^aja 


(C.6) 

(C.7) 


Kuu = [V^AV^Ra + VuBTV^A + ■ VRJ 


~ I'^a + ^iBaV^uA + VlyBaVfiA + XyBaX^Ba — Q^yXpBaX^A\ , 

iT;. , = -5Thf] e‘^^^+^»^XBa ■ VR, 

^bJalb 3 a khlh 


b • 


(a / b) 


(C.8) 

(C.9) 

(CAO) 


The nonzero components of the Ricci tensor are 

Rfiu = Rpy + (m — 2) V^V yA + V pAXyA — Qpy (VA) + g^yX'^A 

— na [XpXyBa + XpBaXyBa + XpBaXyA + XyBaXpA — gpyXA • VRa] 

[-("^ - 2)VA • XBa + V^R, + n,VR, • VR„] , 


(C.ll) 

(C.12) 


and the scalar curvature is 


R = + e-'^^-TT'^^ + 


2{m — 1)V^A — (m — l)(m — 2) (VA)^ 


- 2 naX‘^Ba - Ua {XBaf + 2 na{m - 2)VA • XBa - UaUbXBa ■ XBb 


The Laplacian operator on the scalar field (pi^^ Ua) is 


(C.13) 


X^cP = e^^Xf^^cP + e-^^-Xl)cP + e^^-Xp KR^ - (m - 2)A] . 


(a) 


(CA4) 
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C.2 Parametrization 2 

Another parametrization of the metric may be useful in other situations. Let’s assume the warped 
metric has the form 


ds'^ = GjMNdx^dx^ = g^u{x)dx>"dx’' + fa{x)hlf (?/(„)) 


The non-zero components of the affine connection are 


IVa VOi. ’ 

T^ia _ ■p^a 
ja^a ja^a ^ 

Mia “ 2 • 

The nonzero components of the Riemann tensor are 


oM _ oM 
^uafi ^uaf) ’ 


1 


-^iafcaL “ -^iafcaL + (^/a) 


OM = -h 

^ia^ja 2^iaja 

rU. = -sU 

Mai' 2 .?a 


^/a 'V^/aV./a - V^V./a 


7r/a“"V^/,V,/a - f-^V^V,fa 




(a / 6) 


The nonzero components of the Ricci tensor are 


R^u — Rfiu T 2 


1 


fa^yja • ^ufa " fa'^k^i'fi 


-u 


2^a ^ fiJa ^i/Ja Ja 


D _ P _ j!l ^(“1 

^*aia — -f^laia 2 Lja 


• V/fe - /-^ (V/a)2 + VVa 


and the scalar curvature is 

R = . VR - naf-^V^fa + {Vfaf 

The Laplacian operator for a scalar field (p{x, Ua) 


V^cP = V^^cP + f-^Vl^)CP + ^naf-^Vfa • V<P 


-u 


(C.15) 

(C.16) 

(C.17) 

(C.18) 

(C.19) 


(C.20) 

(C.21) 

(C.22) 

(C.23) 

(C.24) 

(C.25) 

(C.26) 


(C.27) 
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Appendix D 

Quartic supersymmetric potentials 


In this appendix, we explain the supersymmetric potentials used in Chapter to model the depen¬ 
dence of tunneling rates on the number of fields. We use chiral superfields ‘I'j, 

$(x, e) = A(x} + i9a^edmA{x) + ]eee9aA{x) + V29'iP{x) - -^99dmi^{x)a^9 + 99F{x) . (D.l) 

4 V 2 

The general form of Lagrangian for a set of chiral fields $ 1 ,..., is [90] 


C = 


eeee 


+ 




-|- h.c. 


ee 


where m, g and h are symmetric tensors. In terms of the component fields, this becomes 

1 


£ - + A-DA + F.F- + \^A.Fi 

T 9ijk{,AiAj F}^ Ipi'lpjAif;') p hijmn {EiAjAfnAji 3AiAj1prn'4^n) T ^iEi T h.c. 
Varying with respect to the auxiliary fields Ff^ and Fj^, we find 

Fk = -Xl-m*,A*-gL,A*A*-hT^,A*A*A*^, 

Ff. — Afc Tni}^Ai gijf^AiAj Hj^j^yij^AiAjA^ . 


(D.2) 


(D.3) 


(D.4) 


Plugging this back in Eq.(D.3) gives 


C 


+ A*nAi - ^miki’ii’k 
-g*^,^Pi^,Al-V{Ai,A*) , 


-m*,,iiii)k - gijkTpii^j^k 


(D.5) 
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where 


V = F*,Fk . 


(D.6) 


Using Eq.(D.4) gives 


V = Xk^k + {X%mikAi + Xkm*f^A*) 

+ fnikrrPjkAiA* + X^gijkAiAj + Xkg*jkA*A*j 

+ + '^*ikgabk^*i^aAb + Xkh*j^j^A*A*A*^ + X^hij^k^il^j^m 

+ Sijkdabk^i-^j-^a^b F f^ijmk'^^nk-^i-^j + hij^i^m^^AiAjAmA^ 

+ Slbk^ijmkl^aAlAiAjAm + gabkKjmk^aAbA* A* A*^ 

+ hijmkKhckA*A*jAl,AaAbAc . (D.7) 


To make the linear part of the potential vanish, we need 


At = 0 . 


(D.8) 


Keeping only therm to quartic order in Eq.(D.7) 


V = mikm*kAiA* + rriikglbkAiAlAl + m*,,gabkA*AaAb 

+ 9ijk9abkAiAjA^Af^ + h^j^yif^rrinkA^ Aj A^A^ + hijmk’kn-nkAiAjAjnA^ 


(D.9) 


We can choose a basis = SikMi^. This further simplifies the potential 


V = \Mk\Mk? + MkglhkAkAlAl + Mlg^bkAlA^A^, 

+ 9ijk9lbkAiAjA*Al + hT^^MkA*A*A*^Ak + hijmkM^AiAjAmAl . (D.IO) 

Let’s define = hijmkM^. Please notice there is no summation here on k and H is still 

symmetric in Similarly let’s define Gabk = dahkM^ where again there is no summation over 

k and the tensor is still symmetric in a, b. Let’s define everything in terms of real quantities as 
follows 


— (pi ? 


Hijmk — 
Gabc — 


Qabc 


Fijmk H“ '^LijYnk 
Pabc iQabc ; 
Rabc “1“ '^Sabc • 


(D.ll) 
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The potential in Eq.(D.lO) becomes 


V = \Mk\‘^(pl + \Mk\‘^Trl + 2Pijk [(l)i(j)j(j)k - T^iT^j4>k + 2(l)iTTjTTk] 

T ‘ 2 >Qijk [ T 4 ^i^j'^k ^'^i 4 ^j 4 ^k\ 

+ 0 i(pj(f>rn4^n [Rijk^mnk “ 1 “ -\- ‘IFij^nn ] + TTiTTjTTm'^n [SijkSmnk F Rijk^mnk ‘^P^ijmr}\ 

l^S-i-fYikSnjk “t“ ^^imk^njk ^ijk^mnk ^ijk^mnk “t“ ^^ijmn ^^ijnm\ 

-|- ‘l(f)i'Kj'KYn'^n ^RnjkSimk ‘^SnjkRimk F Lijfim ^Rijmn] 

+ ‘^'^i4^j4^m4^n ‘^RnjkSinik “1“ Lijnm — 3Lij 

mn ] . (D.12) 


This is the form we used for the calculations in Chapter We chose the real and imaginary parts 
of Qijk and hijki in the range [-1,1]. 
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